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FLOW OF VISCO-ELASTIC LIQUIDS IN TUBES 


| .. BY 
J. N. КАРОВ Амр Buasur (бок, Delhi « 
(Received—A pril 18, 1960) 


Abstract. In the present paper, we have obtained the condilions under which e rectihneer flow of а 
vísco-elastic fluid, which is isotropic in its state of rest and for which the stress components are expressible as 
polynomials in the gradients of velocity, gecon acceleration, ...... nth acceleration at the point coosidered, can 
be maintained by а uniform pressure gradient. We also show that the flow of such fimds in tubes can be 
identified with those of certain non-Newtonian liquids. In addition we have obtained expressions for the 
dissipation of energy for certain simple flows and use these to derive certain restrictions on the functions of 
scalar invariants occurring in the expression for the stress tensor from the fact that the dissipation function 
should be essentially non-negative. 


Introduction. The relation between stress and rate of deformation tensors for an 
isotropic flow in which the stress at any point is assumed to depend only on the velocity 
gradients at that point was considered by Reiner (1945) and Rivlin (1948). If the flow is 
incompressible and the stress components are expressible as polynomials in the velocity 
gradient only, then the stress tensor can be expressed in the form. 


ty = Ody + jdady — poy, (1) 
where dy is the rate-of-deformation tensor defined in terms of velocity vector by 
dg =з (003 ttj), dy = 0 (2) 


p is an arbitrary hydrostatic pressure and comma denotes differentiation with respect to 
one of the co-ordinates. Also 0 and у are polynomials in invariants I and 111 defined by 
П = ду Ш = дуда (8) 

More recently, Rivlin and Erickson (1955) have shown that if we assume that in & 
visco-elastic fluid which is isotropic in its state of rest, the stress components are expre- 


ssible as polynomials in the gradients of velocity, second acceleration,............ (n— 1) 
acceleration at the point considered, then the stress matrix T = ||ty|| may be expressed 
ав a matrix polynomials in n kinematic matrices Ај, 4,,......... „А, and whose coefficients 


are expressible as polynomials in traces of products formed out of these matrices. The 
kinematic matrices are defined by the equations. 


(r) 
A, = lláy || А (г = 1,2,...... yn) (4) 
(1) 
Ag = vij t 51 = 24% (5) 


(т) (г) (г) 
Agtt) = S4, С) а у Дыт Ammi Tr = 1,2.. n 1] (6) 
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In a later paper, Rivlin (1955) considered certain simple types of steady-state 
laminar flows viz. rectilinear laminar flow, torsional flow between two parallel plane 
discks, helical flow in the annular region between the two coaxial cylinders, for these 
mere general visco-elastic fluids Та all these cases, A, = 0 for r > 2 and 


T = —pl+a,A,+0,4,+%,4,7+0,4,7+0,(4,A,+A,A,) 
ta (4,4, + A,A,7) c (А, А, + A, A, E s (АА,  A724,7) (7) 


where I is the unit matrix and a,,a, ......... ‚® are polynomials in the ten scalar invariants 
tr. Ay, tr. Ду, tr. Ау, tr. Ag, tr. Аг, tr. Ал, tr. А.А, tr. АРА, tr. A,A? and tr. А,2А,2, 


Ib is obvious that the fluids characterised by (1) are particular cases of those 
characterised by (7) (i.e. when о, = 0, а, = as = а, = а, = а, = 0). 


Rivlin and Green (1956) have considered the flow of non-Newtonian fluids charac- 
terised by (1) through tubes and in particular show that in general, unless the tube is 
circular, a rectilinear flow cannot be maintained by a uniform pressuie gradient and 
secondary flows can arise in the absence of body forces. This conclusion was earlier 
reached by Erickson (1956) and later examined by Stone (1956). 


In the present paper, we show that a rectilinear flow of a liquid characterised by (7) 
can be maintained by a uniform pressure gradient under the same conditions as obtained 
by Rivlin and Green (1956) for the less general flow characterised by (1). In fact, we 
prove a more general result that the rectilinear flows in tubes of the visco-elastic flows 
characterised by (7) can be identified with those of certain non-Newtonian fluids charac- 
terised by (1), 


We also obtain expression for the dissipation of energy and obtain certain restric- 
tions on a's from the fat that this should not be negative. 


2. The Basio Equations, Let the rectilinear flow be given by 


v, —0, 0 = 0, v, = 2,9) (8) 
so that 
0 0 f 2f, Зы 0 
А, = 0 0 f, and 4, — 2] 2f, 0 '9) 
в fy 0 0 0 0 
and А, =0 [r2 2] (10) 


From (9) we easily obtain the matrices A,*, 4,2, А,А,, A,A,, A,A,”, A,*A,, 4,°4,, 
A,A,’, Ar A’, А, А,“ and their traces which are given by 
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- tr. А, = 0, tr. А, = fg? + fy’), tr. Ay? = Ufe fy), 
tr. А,“ = Alfo + fy)’, tr, A,A,? = 0, tr. A,A, = 0, 


tr. A,A, = 0, tr. 4:34, == lfa? + fA)? (11) 
ir. АА = 4(fa7 +), tr. A, = 
tr. А, = B(fo + fy’)? 
Also the components of the stress tensor are given by 
А = —р+ (203 + aa)fa? + (а, + ae) fe (ја + југ) + Вавјаз(јаз + fy?)® 
а = >p (224 +0) + (а, 0 ју (ја + fy?) + Bo f (fa? + f) 
tss = -р+а, (е + fy?) 
у (12) 
tia = (2a, + а) ау (а, + авјјвју(ја“ +.) + Ва fuf (fu) +f)? 
tis = а] + даба (ја + ју )4arfolfa® + fy"). 
у = tify + да ји (ја + fy?) + Anu (fa +f) 
From (11) it appears that ал, изу аз......... ‚ аз are all function of fg? + fy’. 
Again using (8) and (12) the equations of motion give 
oX +22 = 6 = Зе (29; + а.) а tla, а ја (fat + fy’) + Bashy d (fs? +}, x } 
uL +4(94 + аварий + fy") + Вазы Гу (fa? + fy?) ®t 
> Op ди: 3472): 2 212 
сҮ *By = Bg =з + аа) Јајо + 4 (a4 аа) fel y (fe +747) FBasfyfo fo + f)? (13) 
P 5, 0n а) + (а, + af (јаз + fy?) + Вар" (fu? + |") (14) 
0+2 = fa Datel fa? + fy?) + да ја ја f) 
y leafy fur + fet) + df Y (18) 


3. Identification with the motion of non-Newtonian fluids. Let us make the 
substitutions 


2, +29,( + fy?) + Aa. (fo? tf) == 0 (16) 

(Qa, + as + (а, + ав) (fo? + fy”) + Bag (fo? + 57 = ү, (17) 

where 0 and у are functions of [м +, since all the e's are functions or the same 
quantity. 
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With these substitutions, equations (18), (15) become 


хъдр 905,49. 18 
e +5 Я Velta fef y) (18) 
id." MED RE mE ) (19) 
+3, 5 (8/5) + = (6f,) (20) 


But these are precisely the equaticns of motion obtamed by Rivlin and Green for 
the motion of non-Newtonian fluids characterised by the coefficients 9 and № of viscosity 
and cross-viscosity respectively. Thus we find that the motion of visco-olastic flows 
characterised by (7) is identical with that of non-Newtonian flows for wl ich 0 and у are 
given by (16) and (17), In other words to write the equations of motion for the visco- 
elastic liquids, we write these for non-Newtonian liquids and replace «, and аз by 

€4 + Зав (Де + fy’) + Фа (ры + ће) 
and (дал оз) + Alay + s) (fa? + fy?) + 8% (о? + fy)” 
respectively. In writing the expressions for the components of stress, the same substitu- 
tions can be used except that in the expression for bss the original expression for оз is to 
be retained. 

à. Reotilinear motion under uniform pressure gradient. From the above identi- 
fication and from the results of Green and Erickson, we deduce that a purely rectilinear 
flow can be maintained by a constant pressure gradient if 

() p= kG i.e, in particular if 
Dals = Da бшщ} (21) 
| 95 €. 

This condition is obviously satisfied if all a'8 are constant. 
or (ii) fis a function of a linear combination of æ and у 
or (ii) fis a function of (z* + у") 

Discussion of other simple flows. ‘Rivlin (1955) has discussed certain simple flows 
of visco-elastic fluids characterised by (7). It can be shown by arguments similar to 
those used above that all those flows can be deduced from those of ordinary non-Newto- 
nian flows by certain substitutions. We give below these substitutions as well as the 
modifications in the components of the stress tensor. 

(a) Rectilinear laminar flow: Let the velocity components be given by 

ор = kv, v2 =0, v3 = 0, (22) 

then the substitutions are 
a, by a; + аи + 4o, k* (28) 
га by (да, Фа) + (а, + ak? + Ва, (24) 
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In calculating £,,, however, the original value of a, 1s to be retained. 


(b) Torsional flow of a cylindrical mass о] a fluid, Let the velocity components 
be given by : à E 


U,— 0, ve = yrs, v,—0, ~- (25) 

then the substitutions are 
ay by a, + 22,75? + 49,7404 (26) 
аз by (да, +a) + 4(&, + a)r’ + Ват (27) 


Here only for tg the original value are to be retained. 
(c) Helical flow in an annular space | 
Неге v, — 0, ve=wr, т; = и, (28) 
and the substitutions are 
a, by а + да ! rw? + u^?) + 4a, (ul? т?) (29) 
аз by (21, + аз) + (ag + а) (720/2 +u?) + Ва (u'? + rw’) (80) 


In calculating the stress components 1,,, tro, trs we use these substitutions, while in 
calculating the other components we retain the original values. 


6. Dissipation of energy. For finding this, we calculate the expression 


4 


in each ease we get 
(a) for rectilinear motion between planes: 
D = К (а, +2a,k? + 4a,h*) (81) 
(b) for torsional flow between two parallel discks 
D = ут (а, + 805027? + Дарт“ (82) 
(с) for helical flow in the annular region between two coaxial cylinders. 
р = (гь 4 w?) a, — да (w^ + иЗ) + атша + и] (88) 
(d) for flow in a tube | 
D = (foi + fy?) [ay + дама! + fy?) + dns + 17) ] (84) 


Since the dissipation cannot be negative, we get that in the four cases discussed 
above, the coefficients a,, «а, «, should be such function of the invariants that 


@,+Qagk? + 4a k* > 0; а + дот" + Дач“ > 0; 
a, + да (u? + 74073) + 4a и"? + та") > 0, (85) 
а: + За (јан + fy?) + Хал (ја + ју)" == 0 
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It is interesting to note that the expressions which we require to be non-negative are 
those corresponding to the coefficient of viscosity in the substitutions obtained above. 
Also the similarity in the above expressions is interesting. 


DELHI UNIVERSITY 
DELHI 
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ON GENERALISED HANKEL-TRANSFORM-—III 


Ву 
Вам Komar, Roorkee, U.P. 


(Recerted—April 18. 1980) 


1. Agarwal (1950) iniroduced the generalised Hankel-transform 


e 


s) = @* | up iri ват јада a) 


and obtained certain inversion formulae (1950, 1953) and the fundamental properties 
(1951) for it. Stankovic (1955) obtained the general form of functions self-reciprecal 
under this transform. I have in a series of papers studied the recurrence relations 
(1954, 1955), integral representations (1957), infinite series expansions (1957), connections 
with other transforms (1956, 1957), convergence theorem (1959), etc. of the transform 


g(y) = || (xy), (с) (ade, (2) 


0 


where J а) is the Bessel-Maitland function defined by’ 
л < (==) 
Ј 
a(n) = 2 атаа y ncmo (8) 


In view of the interesting nature of this generalised Hankel-transform, I give 
herewith certain rules for it. Some of these rules are similar to those of the Hankel- 
transform while the analogues of others do not exist in the Hankel-transform theory. І 
have also given some examples to illustrate their use. 


In what follows, we sball call g(y) the generalised Hankel-transform of f(x) or 


J; „-transform of f(x) when they are connected by the relation (2). 
2. If the J, -transform of f(x) is g(y), then we have the following rules :— 
и . 1 
(1) Jy,,-transform of Каз) 18 59910). 


(2) J, -transform of xf(z) 
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= (L/py) [A — ид (у; А-у) -9(9; Ар 1, v)], 
(L/Ay)[g(y; Хо Ту 1) tag (ys А + psv +2) 1, "— 
ЈОДА Dgly ; Atl iv +1) ^ ng(y ; ААТ, у +2)]. 


(8) J,,,-transform of “ f(a) 


= (yl [g(y; à—1,u,v-1) + ng(y 5; Хи ЊУ] 
= (yl) lA- ugly; Х—рању— 2) – g(y ; X n1," 2)], 
= y[(A+1 gly; Atl, 71) -ug(y ; Xl вр, ]. 


The rules (2) and (8) are obvious transformations of the following result due to 


Wright 11988) : — 
огфр+р = 98–1 + (1—0)95, 

g Ф zi 

Fus = #0629) = È ратая" 

ni 

and hence | Ф,(ж) = Ф(и,1+А; —а), p>o0. 

А H 
(4) Ja, transform орг") is y^"g(y; А,д,у + m). 


в А т " 
(5) Ji,transformo f x™f(x) is (Pu ps [y79(y; Х—тииу)], 





provided m is any positive integer including zero. 

о For, ду Атыра) = f (аа Ге», 

t у 0 
then, by virtue of the relation (Wright, 1988) 
(ddz)eg = Cp: 

° д 
f tinyay), 
0 


E fast; Атон} 
whence, оп generalisation, we get the required result. 


" o 
© ЈА, -transform of tfla) ia y f agla; oed, 
ra 


provided 0 < в < L with an additional condition (А + д) > —4 in case р = 1. 


wf wards = у f еә | летучая 
Li 0 y 


For 
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= J Guy) 7 fede. 
0 


(7) J, transform of 1 е9) 


с 


Hu bon] т" "18 4, -a(y] (2n)*(g(n; А — уд + идр ћуја 
0 
Since f(y) = (ду)-"е- гр, (а (у), 
= (пј ај рт 1478, 


= 9(p) 
and therefore 
р-^-'(р-*) = (2 југа лат P 
ра te 
= (тј) #247) ay, plat), 
provided R(A-- фри р) > —1. 
Therefore (Kumar, 1958) 


/ ЈА (20) (95) 498 6789), (а|(ду)%)йу 
0 


= (Ән ибаа). 
The required result can now be obtained by proceeding as in (0). 
(8) Ji transform of абоза f poetas mala) (Да 
ig pl\c)g(y, А+с, и, | 


| А -= j MT 0 -— /; ed 
For | етае | осона) oet jan 
© € | 
= yf £c ovre) | хоч) 7 Yn фи) et, (zy)dad£, 
0 0 


= „го f (вул оков 


0 
= АГ(е)9 (у; A оу). 
9-~2057P—I 
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The z-integral has been evaluated by term integration, a process easily justifiable. 
n А & „% 
(9) Ja,,-transform ој z= f || Mis / | ё" (баба, dE.. dëm 
о 'o во 
is g(y; A—my,u,v— т), where m is any positive integer including 0. For, И we write 


Р, (£m) = || s / РЕН || : || " p-"fOdeat, dt, un db 


oo 


then / су а f №. (Са ја С dz, 
0 


0 


v f 762 | ла, 
0 fn 


=y f Лун) F (Emde my 
0 


By repeating the above process, the required result is obtained. 


3. We shall now find out the rules of finding out the J; -transform of the function 


involving the differential co efficients of f(x). 


(1) ЈА transform of ў (а) = у[9(у; X iu в, у») -vgly; А,д›—1], 


provided aJ (ау) б) tends tu zero both ав z — 0 and as 2 > со, 


This can easily be obtained on integrating by parts and then making use of 3(A) - 


a 
The J,,,-transform of higher derivatives can be obtained by the repeated application 


of the above result. For example, в 


"(уз Ашу) = (у; AF pav) уд (у; Xov DE 
= gly; А+дии у) —259(у; A u,uv — 1) 
Tw -19(y; A,u,v—2)]. 


в ал“ 
(2) Jy,,-transform of (2) forf(x)} = y"g(y; A+ та иу), 


provided wd (у) (=; Ё ү" fz*/(x)} vanishes both at zero and at oo. 
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< а d 
For (уа (y) (a= Jta 
А dz 


=y f «чле, 


by virtue of the relation 2(B) and the conditions of the rule, Hence, on generalisation, 
we obtain the required result, 


(8) Jhy-transform of a + Y df 


T 
= y'[g(y ; Хари v) –уд(у з Ау 1), 


provided (ау)! (а) and =» (ey) (2) vanish both at 0 and at со, 
SDN 7 x Уз : - 
9 в 2 = в А 
For | yep PF an = -yf erre) ать еу) dos 
a А 


0 T 


Therefore 


= 


dz sdr 


с Е ИГ. 

J етеу E + Уа = у | ст араа, 

0 | } 0 | Фе 
= y'[g(y 5 А+, иу) —vg(2 ; At pv —1)]. 


(4) The rules given in sections 2 and 3 can be employed to find out J, ,-transforms 


of certain functions and to evaluate certain integrals. Ав an illustration, we give below 
few examples, 


Example, 1, For v = 1, 2(6) gives 


[ Ji(ry)n7 tene Hedy = пћу“+ i- (ау. 
“a 


Therefore, by rule 6, we have ~ 


Ф < 
= - e B г ia 
f iere tentes = кї | maa, 
0 ^ 


= te 
= hy J ruetntt, 
y 


. = Эту (ау). 
а А 
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i ia 
Hence у, -transform of 2-7 2126—4169 ig m yd (ау!) 


€ + 
ant ЈА (ayt), 
a 


and Ледовое = 
0 


provided 0 < p zz; 1 with an additional condition R(A)>—3/2. 
Example 9. Since 


i A 
DATED 
0 


then, by rule (5), we have 


2 А " d^ [1,* 1 
] а Ј ,lay)de = (-) ds [но], 


where т ів any positive integer including zero and 0 < д S 1 with an additional condition 
Ки) > 2m +8/2. 
Ezample 3. From 2(C), by the help of the rule 5, we have 


| euer D, sto] Qn) 
0 


=( ) ( ) Th |» у - ( y | 
= mov = (п) iv- J ме + | 
2 к Att pomtip а 


Example, 4. If we take f(x) = 1, then the Ji, transform of 


ln Го) 
E) CESET 
ón making use the following result due to Gupta (1048) :— 


ГИ +v) 1 


S д 
vJ dz = = . =, 
] (ey). (zy)dz TOFA p-m) y 


valid when 0 < u < 1 and R(v)>=1 with an additional condition R(À=2w)>-—$ in case 
=L. 


1 m 2 
Also af ni à- EE f (a) PI, нзрт) der, 
0 


n I(L+$v) 1 
DESETIN 
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Hence, by rule 7, we have 
= 
-yt 1276+1 
велен (удал = Г). ETE 
[^ ьа = ech T 
valid for R(v) > 4, 
Example 5, With the help of the integral (Kumar, 1955) 


с 


У 2 
f ерене? „аде "Неа буда 
6 


в 
= ис (а 1) ареала), 


where m, n and с are any positive integers including zero, and the rule 8, 16 can be 


easily shown that 
œ a Ф 
| artimi? (ay) f Ee-ee tHenl2gtdgde 
= 2-»-imiT(o)rty-in737,4 plyt). 
This result ean be easily verified by changing the order of integration and then 
evaluating the £-integral. 


DEPARTMENT OF МАТНЕМАТ1ОВ, 
URIVERSITY OF ROORKEE 
ROORKEB 
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SOME FUNDAMENTAL CONSIDERATIONS ABOUT 
MASS IN GENERAL AND ABOUT 
NEGATIVE MASSES 


Bv 
M. P. Verma, Varanasi, U.P. 


(Communicated by Prof, V. V. Narlilar— Recewed—May 12, 1960) 


Abstract, Ivery particle can be assigned three numbers representing its three masses 012 , the inertial, 
the active gravitational and the passive gravitational. In S:etion (1) we consider qualitatively the motion of 
two particles due to their mutual attraction or repuleion when their masses can be positive or negative and 
indicate twenty different possibilities In Section (2) the motion of test particles and of light pulses in 
Schwarzschild's external field and 1n the external field of a matter dipole has been considered The new 
possibilities regarding the nature of mass considered here lead to unexpected consequences regarding the 
velocities of test particles and light pulses in а gravitational field. It has been shown that the velocity of a 
light pulse exceeds the fundamental velocity c (=1 in natural units) and 13 not bounded above. The test 
particles, while they always wove with velocities less than the velocities of light at the points of their 
instantaneous lozation, acquire velocity maxima which are also unbounded above. 


1. In classical physics the concept of mass is introduced in two different contexts, 
firstly in the second law of motion 





P = mf, (1.1) 
and secondly in the law of gravitation 
ray T у (1.2) 


While the mass m appearing in equation (1.1) measures the capacity of the body to 
resist motion due to a force acting on ıt, the masses m, and т, appearing in equation 
(1.2) measure the capacities of the two bodies, separated by a distance 7, to attract or get 
attracted. We thus arrive ab three distinct notions of mass which can explicitly be 
stated as follows. 


(A) The inertial mass of a body is a measure of tts capacity to resist motion when a 
force acts on it 


(B) The active gravitational mass of a body is a measure of its capacity to attract 
another body, 


(C) The passive gravitational mass of a body is a measure of its capacity to respond 
to the attraction of another body. à 


If we distinguish the three masses associated with a body by giving suffixes 1, a, 


16 M. P. VERMA 


and p, the equations (1.1) and (1.2) become 


Р = тј, (1.8) 

Fig = у 2191р, | (1,44) 
p 

jy me, (1.48) 


where F'y measures the force of attraction on the body j due b` the body i. 


The accelerations experienced by the two bodies due to their mutual attraction wil 





be given by 
ћ= Pa = Y maa P (1.5) 
Mi T тү 
T. = Е,, = X m,a a (1.6) 
ma T Ma 


We observe that (1.5) and (1.6) involve only the ratios of the passive gravitational to the 
inertial masses of the bcdies. 


The units having been suitably chosen Newton’s third law of motion implies the 
equality of the active and passive gravitational masses of a body which immediately 
follows from equations (1.4А) and (1.4B). The equality of the inertial and the gravita- 
tional masses of а body is axiomatic. Thus in classica! physics the three masses of a body 
are represented by the same number, 


General relativity does not make any distinction between the active and passive 
gravitational masses. ‘The equality of the inertial and the gravitational masses is 
postulated in the principle of equivalence which plays a fundamental role in the formula- 
tion of general relativity. Experimentally the ratio is found io differ from unity by а 
term of the order of 1075. The gravitational mass appears for the first time as a constant 
of integration in Schwarzschild’s exterior solution. This mass is essentially the active 
gravitational mass, The inertial and the passive gravilational masses do not appear 


anywhere explicitly. However, we can still give meaning to these terms. In the 
operational definition of the field viz., 


Field intensity = Lt, Foree experienced (л) 
mass > 0 mass, 


the mass appearing is necessarily the passive ртауіЊа лопа! mass. The snertial mass can 
be given the sume meaning as in equation (1,1). 


The force appearing in equation (1.7) thus involves the inertia! mass so that we get 
on the right hand side of equation (1.7) only a ratio of the inertial to the passive gravita-- 
tional mass. Now the principle of equivalence demands that the ratio of the inertial to 
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the passive gravitational miss should be the same for all bodies, and that in a way 
explains why these two masses do not explicitly appear in the equations of motion in 
general relativity. 

Recently the possibility of negative masses has been suggested and discussed by 
Bondi (1957). А number of possibilities arise.for the accelerations of two bodies under 
their mutual attraction or repulsion. We give below tables indicating the different 
possibilities on the basis of equations (1.5) and (1.6) which could be taken as a guide fo; 
the corresponding relativistic problems. The ratio of inertial to passive gravitational 
massis taken to be +1 or —1, ‘Towards’ appearing in column ‘Body Г’ will mean 
towards body 2 and so on. 


ТАВЉЕ I 
mi/my = +1 
MASS ACCELERATION 
BODY 1 BODY 2 BODY 1 BODY 2 
8. Inertial Passive Active Inertial Passive Active 
No. gravita- gravita- gravita- gravita- 
tional tional tional tional 
1. + + + + + + Towards Towards 
2. + + + + + = Away Towards 
8. + + + – — t Towards Towards 
4. + + + — — — Away Towards 
5. + + = + + - Away Away 
6. + + Е = B + Towards Away 
7. + + – – – — Away Away 
8, - – + — - + Towards Towards 
9. - - + — Е – Away Towards 
10. = – – – – – Away Away 
TABLE П 
ттр = —1 
MASS ACCELERATION 
BODY 1 BODY 2 BODY 1 BODY 2 
S. Inertial Passive Active Inertial Passive 
No. gravita- gravita- gravita- 
tional tional tional 
1 + - + + — T Away Away 
2. + — + — + + Away Away 
8. + - + + — — Towards Away 
4, + – + — T — Towards Away 
б. - + + — + + Away Away 
6. — + + + - — Towards Away 
7, — + + — 4 — Towards Away 
8. + - — + - — Towards Towards 
9. + — — - T — Towards Towards 
10. ЁР + — = + – To vards Towards 


8—9057P— 1 
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We observe that for the ratio mim, = +1, а body 1s always ascelerated towards а 
positive active mass and away from a negative active mass, as observed by Bondi but for 
the ratio m;/m, = —1, а body is always accelerated towards a negative active mass and 
away from a positive active mass. | 


For each ratio тутр only two of the three masses of a body can be-given any of 
the two signs independently, Hence there аге in all 2x 2 = 4 possible combinations of 
signs for the three masses of a body. ‘Therefore, the number of possible combinations of 
signs for the two bodies will be4x4 = 16. Tke number of distinct cases is reduced to 
10 on account of symmetry. 


a. We now investigate the situations arising out of importing a negative active 


gravitational mass within the frame-work of general relativity. In Schwarzschild’s 
exterior solution if we replace m by – т we get 


=t 
ditus -(1+ =) di пав -Pasféde + (1+ 2" Уди, (2.1) 
т 
Tbe metric is obviously regular everywhere except at the point r — 0. For the radial 


motion of a test particle in this field, following the method indicated by Singh and 
Pandey in a recent communication, we obtain from the equations of the geodesics 


у _ _ 8m 2m N41 т 2m, 

T + (1+) +i (1+ 2") (2.2) 
where a dot overhead indicates а differentiation with respect to time, The first integral 
of this equation gives : 4 

ji = Zm (да e 2m V a 

r= ер) G Z 81-4), (2.8) 


whore т = 0 forr = а. Equation (2.3) shows that forr to be real r has always to be 
greater than a, 


16 is easy to see that т reaches a maximum for 


6ma 





^ Am—a’ (2.4) 
and this maximum velocity is given by 
Е 2 2m 
Tmax = 8.8 1+ а Y (2.0) 
which is greater than unity for 
ac 4m 


SOME FUNDAMENTAL CONSIDERATIONS ABOUT MASS IN GENERAL ЕТО. 19 


‘and does not have an upper bound. Thus if a test particle statis from a point 


_ z> 4m 
=“ < g8i-3 


it will aquire a velocity grealer than unity and the maximum velocity will tend to infinity 
as a tends to zero. 

Equation (2.4) also shows that the test particle attains a velocity-maximum only if 
it starts from a position т = a «24m. For the body at the origin to have a mass equal in 


magnitude to that of the sun this will require the numerioal value of its volume density 
to be greater than 2.1015 gms/e.c. which is unusually high. 


For the radial propagation of light we have to put 
ds = 0 with 0 = 9 = 0, 


‚ во that from equation (2.1) we geb the radial velocity r of light $o be 


#=1+®®, (2.6) 


which is greater than one for all values of т and is unbounded above. However, equations 
(2.8) and (2.6) show that the radial velocity of light pulses is greater than the radial velo- 
city of a test particle for all r. 

Similar results can be obtained for Bondi’s general relativistic model of an 
uniformly accelerated matter dipole, The metric for this model in a co-moving frame of 
reference is 





da? = 679 dt? —e-**?ec(dr? + dg?) — re 2402, (2.7) 

3 _ OoOo m = та у 

where Ф = &log([r? + $?]* +8} ОВ (AJ. (2.8) 
с = 41001 +42(7 + 27) ~ 1j, (2.9) 


for pointa external to the material distributions. m, and m, are the masses of the two 
spheres of radii a, and а, placed wilh their centres on the z-axis at s = h, and s = h; 
respectively, Ћ, and №, (h> ћу) being positive, For а, << h., а, << (h,—h,) and 
в << (h,—h,) we have 


DET 
TR, = DA , | 
| (2.10) . 
та = и . | | ` 


The singularities of the metric are located on the s.axis abs = h,, я = h, гой 
à < 0, the last one being purely a coordinate singularity. 
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For this metric the equations of the geodesics are 


р Suy. (us e 


4o 2,09 ЈЕ «(er - dr da 241 
+e 20021 27520, (2.11) 
4:0 1 – 98) ar drd6 „Әә ds dé : 
== 49 = Е .12 
ds? t (3 Эг/ ds ав Se ds ds 9 ote 


die (8%) (27) - (46) Va n eno 
ds? * Os де ЭЕ ds н Pa Os ds 


3 
= S 2) (@- enar ds ds 
4s -2, Ө? dt = р 
S OsX da ix Or er/ds ds 0, ant) 
аң 89 dr di 299 dz dt 
m wp QUT Se =0 2. 
р ds! “Orde da Os de de ~ O ‚н 
If initially js = = 0, 
» | 436 
we get from equation (2.12) У = 0. 


Hence motion is restricted to the plane 0 = а, 


Now equation (2.14) can be written as 


0192-7520 2. 
de? dade ý (449) 
which gives as its first integral 
dt 
SE = Де-®, - (2.16; 
With the help of equation (2.15) we can reduce equations (2.11) and (2.12) to 
OR es 3799) . = p 4 -2,09 (22-2) • = 
T ib Fr (1? — 23) + eto a. t? Sr a) 0, (2.17) 
and — 229 + (2 -emye —32)4 gie ut + «(8 - oe ra = 0. (2.18) 
Now for r — O it ean easily be seen that 
Po, 


or Or 
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Therefore, if we put т = 0 when г = a as an initial condition equation (2.17) will give 


T = 0. 
Equation (2.18) then gives the equation of motion of а test particle along the a-axis 


which can be simplified and integrated to give an expression for a, However, 2* > is. de- 
termined more easily from the metric (2.7) with ths help of equation (2.16) as 


а Aa? _ 42a | [4.2 Е (a — 2)(aa — h,h,) 
Е expl SERRA exp oe cess Wee BER St 210 





(ha фи h,)° р 


where we have pub д = 0 for # = a and а = LUN 


Let a be greater than h,. Then for 2 to be real ғ has always to ba less than a. А 
test particle, therefore, fells from a position я = a > h, towards the dipo:e and reaches 


its maximum velocity а given by 


- 2 - { Qa \ 
Bmax = — 8 exp} -- cM. _ 2.90 
8t 5 (z (a—h ( 


at a point в such that 


z «(a — a) (az — hyh,) 4 
р = п __} = 8, 
а { (a—h,)(a—-h,)(a—h,)(e- al Ы 
For a,2>h 
s 
z >} 
- 2а 
or 8 > т 


Equation (2.20) clearly shows that Bmax does not have an upper boundt 


For light the velocity 2’ is found to be given by 


" = 3s oxp{ - = TES CES сл) (2.21) 


+ ltisan inherent property of the uniformly accelerated metrio ds? = 24di?—daI—dyi—dal of which 
(2.7) is в transform for mj=m,=0 that test particles and light pulses do not have an upper bound for their 
velocities in the space-time described by it. 
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which is also obviously unbounded above. In this case again equations (2.19) and (2.21) 
show that the particle velocity is less than the velocity of light at any point along the 
axis of s. 

The equations of the geodesics ignore entirely the considerations of mass of a test 
particle. The treatment given above, remains valid when for the test particle the ratio 
mim, = +1. The case тијтр = —1 and also the case when та 0 can be included 
only if the problem is solved as a two body problem so as to take account of the masses 
of both the bodies. This requires that the problem of n bodies be solved in general 
relativity by deriving the equations of motion from the field equations themselves in the 
general case of three types of finite masses associated with each body. 

The above investigation shows that the introduction of negative masses within the 
frame-work of general relativity contradicts the fundamental hypothesis regarding the 
velocity of light. It needs a careful examination as to how the concept of negative mass 
can be accommodated п a relativity theory. 

The author's thanks are due to Prof. V. V. Narlikar, Head of the Department of 
Mathematies, Banaras Hindu University for having suggested the problem. His thanks 
are also due to Dr. K. P. Singh for helpful discussions. 
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1. Introduction. The integral equation 


to) =] eftt, Rip) >0, (1.1) 
0 


represents the classical Laplace Transform and the functions фір) and f/t), related by (1.1), 

are said to be operationally related to, each other. 9(p) is called the image and f(t), the 
original. Symbolically we write 

olp) = f(é) or f(t) = e(p). (1,2) 

Many generalizations of the aforesaid transform have been given from time to time 

by different mathematicians like Meijer (1940, 1941), Boas (1942), Varma (1947, 1951) 

and others. In the present note we take up a new generalization on the lines of Meijer 


and Varma, We give some theorems analogous to the well-known Goldstein’s theorem 
(Goldstein, 1982) and also establish the Uniqueness theorem and an Inversion Formula. 


2. We consider the integral relation 
e(p) =p | етн, рр) 73-й, Rip) > 0, (2.1) 
0 


provided of course, the integral on the right exists. Here also we call Ф(р), the image 
and f(t), the original and symbolically we write 


k k 
— — А 
plp) A f(t) or f(t) А e(p) (2.9) 
<— ~ 
т m 
It is evident that (2.1) reduces to the Meijer transform if we take 
А = К, (2.8) 
and to the generalization due to Varma if we assume 
А = -—m, (2.4) 
Further if we take A=k=> + т, (2.5) 


we fall back upon the original Laplaee transform (1.1). 


3. The Uniqueness Theorem. In case of the Laplace transform it has been 
established by Lerch that the function 9(p) is the image of a unique function f(t), Before 
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establishing a similar theorem in connection with the present generalization we prove the 
following lemma. 


wo 


Lemma, If / e РУ, a (pt)teF(t)dt = 0, (8.1) 
then | F(t)=0, (8.2) 
provided that 

(i) F(t) is continuous т t > 0, 

(i) F(t) = 0(), Б(р фт +8)2> — 8, when t is small and 

(i) R- m)». 

Proof. We multiply (8.1) by p'"-1(p—5)*7"-9? and integrate with respect to p 
between the limits (z,o), а 2 0 


Then we obtain J yn-ip аваар | а Три (ptjt^F(t)dt = 0, 
s 0 


oro f erat | e- tW p, mipi) p™ Кр — в) "ар = 0, (8.5) 
0 . 


ihe change in the order of integrations being justifiable by absolute convergence of the 
integrals under the conditions imposed. 


Now evaluating the p-integral, we obtain 


T'k—m-1)Zk*n-i || eie F()dt = 0, 
0 


Г ачта = 0. (8.4) 
0 
Henes by Lerch’s theorem F(t)=0, (Lerch, 1908). 
Thus having established the above lemma, the proof of the Uniqueness Theorem, 
stated below, follows as a direct consequence. 
Theorem 1. Let f,(t) and falt) be continuous in t > 0 


k 


and h(t) X ep), 


—— 


k 


— 


and also fit) A ep), (8.6: 
PM 


(3.5) 


т 
Then fat) = f(t). (8.7) 
4. In connection with the Laplace transform, Goldstein (1982) has proved the 


theorem : 
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"у Фр) f, нл) 
and т) + falt), (4.2) 
then ] sort = jw (4.8) 


provided the changes in the orders of integrations involved аге justifiable,” 
An analogue of the above theorem in the present case may be given as follows. 
Theorem 2, Let /,(1) and f,(t) be continuous in 12» 0, 


k 


— 


and ep) А |), (4.4) 
i . 
m 
k 


— 


and ep) А falt). ME (4.5) 


< 


т 
Then f «i03 Јл (4.9) 
0 0 


provided the integrals converge and 
(i) fit) = v(t), RO, £m-Ay»-1, when + 18 small, 
and (i) fat) = O(t9), R(8,+m—A)>—1, when t is small, 


The proof 1в quite simple, being similar to that of the Goldstein’s theorem, and 
hence we omit. Е 


Below we give two more theorems analogous to the above theorem but slightly more 
' general. 


Theorem 3. Let f,(t) and f,() be continuous in t > 0, 
) zu k 


T в.в) м fill) | | (4.1) 
pn - 


cc Ё 


and 7o cep. | 7 (4.8) 


© 


Thon [ stor oeil | e cat opa, - (4.9) 


0 0 
provided the integrals converge and 


(hf) = 065), RQ, +т—А,)>-—1 | 


and () falt) = 005), Rat mA) –1 
4—9057P—1 


when t is small, 
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Proof, We have 


T= | ој без 
J 


= f деста | commutat; hf liy, 
0 0 


Бу (47). 


(4.10) 


Now changing the order of integrations, which is justifiable under the conditions 


imposed upon f,(t) and f,(t), we get 


I = | f oyrsdy | e WOW, a s (yt)(yl) Hf (Edt, 
0 0 


= [ бу, 
0 


by (4.8). Hence the theorem 
Evidently theorem 2. follows from theorem 8., if we take A,=A,. 
Theorem 4. Let f,(1) and f,(t) be continuous int >0, 
k, 


— 


апа Ф,(р) А, f(t), 


and Pal p) А falt). 


та 


kı К, 


Е di мај i 
Then | 1.0%) т A, (р) M || fa(at)e, (t) 
0 << <—"o 
a m, Ms 
provided that (р) existe and - 


G) ће) = 008), R(5, x m, 1, —1 | 
(ii) falt) = O(t), E(5, xm, —Х > —1 
and (Hi) 9, (#)t4~1 and o,(t)ih-1 belong 10 L(0,oo). 
Proof. Let us consider the integral 
vp) e pf «cimus tontos) 44 | је 
0 0 


=> | «02 | ccm реро) 7f tas, 
0 0 


(4.11) 


(4.19) 


(4.18) 


(4.14) 


when t 18 small, 


(4.15) 
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the change in the order of integrations being justifiable under the conditions imposed 
upon f,(t) and e,(1). 
Hence, in view of (4.12), we obtain 


I i dt | 
Фр) = [ a(t), (р) т | (4.16) 


о 
Similarly we obtain 


s" (p) = pf ови ара рвача | falto lt) 
0 0 


= / $: (И Фа(р/%) E | (4.17) 
0 


Now writing p/t for t in (4 17) we obtain (4,16), Hence we conclude 
Ф(р) = e"(p) = e(p)(say). 
Evidently theorem 2, foilows from theorem 4, by making use of the Uniqueness 
theorem (theorem 1.) proved above, if we take 


ky = k Ay = Ху тү = ms. 


5. Examples. (1) Let - ft) = t, Важт-^) 2 —1, (5.1) 
and falt) = #«+2-16-80 (5.2) 
Then, according to thecrem 8, 9,(p) = Py(atm—A, +1) pi (5.8) 
Га - А-А, +1) ' 
and Palp) = 2Bxpt- sk (2(8p)9). (5.4) 
Hence / Фа (0f, (D0«7-1dt = Г (аот — A, 1)g*7 7271, (5.5) 
0 
an J sif atit = Dem А ege, (5.6) 


0 
which is true according (о theorem 8. 
(2. Let f(t) = tt, К(а тА) > —1 and fa(t) = t^e7?', R(u xm Ау) >—1, (5.7) 


Da(otm-—A + 1) 





Then according to theorem 8, — €,(p) = Tía-k-A +1) p", (5.8) 
_ icd; pt+m—-Ayt1m—k; 
«nd Фа(р) = Pius cA E) Е a +1 ( | 3i) 
Ee-k-A, t1) (B+ реттом p-k-A*1; BrP 
R(p+28)>0. (5.9) 
Hence Í (ПИ 


0 
t Ty'A + ага + ВУЧА — В). 
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—Tepim-ASD fe UT mh 8 ) di 


Ifu-k-Mt*1 4 д—К—А+1; ё+ї/ (prie 


1 Tm-AÀ,tl,gm-k; 
2 Г, (ит -А +1) Я! п-ов" з t di ; 
Ги — k - ,- 1) J Е ЕГ фе -А-А+ 


-I(-m- + (а m —A, T 1). (а -A,—À4 а), т) 
Dp Во Ху I)e Mns Е p-k—A,1; 
R(u--À,—A54—a) > 0. 
_ Vyletm—-Ay + ЮГА, —А—9), 











(5.10) 
I(x—k—2A, + Вк 
5 cat T,(o+m—aA a СЕ 

tf (01s S1gi = + 1 Btieti- —a-1 di 

and | СЛАНА i + e" Pts d 
_ Гаја тА 1) Ги  À,— Ха—а) (5 11 

Га-Е-А, +1) ВЕНА А-а. ? | 
which js again true according to theorem 8. . 

(B). Let f. (t) = t* and БО = ttime Pt (am, А) >l. (5.12) 
Then according to theorem 4,, - 9,(p) = = eee (5,18) 
and Фа(р) = 28%pi -Ak am, (2(Bp)#). (5.14) 


Hence | кад = [ порно 


0 


=? f è kym, 809 7d, 


0 


= Гат А, явь, (а 7, +1)>]R(m,)| (6.15) 


and r= f 09,04 
0 


oo 


2 Гк(а ет A + 1) аль, af в ВІЗ [Ak о 2dt 
0 


Г(а—К,—А +1) 


\ =] 


Гат, —2A, +1) 


E 1 ВА, 8-1 ~tte—h,—A, 
пен пере 0 ] Жеш 


l.(a-m,—A ОГ –Е, А a les 
= ағыт A p n „Ва-Ња- а -1. (620) 


` 
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Now we obtain 


Ф(р) =p Г +, (ра) (ре) ~“ (а) аз 


0 
Я Ј immu uus pr) dde 
0 


= Га (а i, —А, + IT g(a т, —-À4 +1) 
Та А Паз Mie 


(5.17) 
and s" (p) = p / ө WIV aus (pe) (ра) ->ti (0) 48 
о 


_ Г (а фт),--А, + DI — Ка -А +1) 
D(a - А, —А, 4-1)g 771p i 





со 
| ву арт, (ра)2 => +%42, 
0 


= Г, (ат. А, +10 Га фт, А, +1) 
TG-E-A +1 MER С, 


=9'(p), by (5.17). 
'8. An Inversion Formula. We have 





(5.18) 


ы 


р)=р [ e VV, Lu s (pt) (pt) ^ tf(t)dt, (6.1) 
0 
Let us assume that Ф(в), the Mellin transform of op), existe. Then 


Фа) = f wp «visis opt ^ tdt, ame tir, 
0 


0 
= | t if(t)dt | o7 „(рр + “зар, (6,2) 
0 0 


the change in the order of integrations being justifiable if we assume that £-«-1f(1) belongs 
to L(0, co) and В(в+т-^)> -1. 


Hence evaluating the p-integral, we obtain 


1 Ty(s+m—-A+1) wget 
Е Та E LEE 


-a-ig — l(8-k—A-1) 
ot f o ке P QW xa (0.8) 


Now applying the Mellin’s inversion formula, we get 


qo Г Rak A1) gy 


9n J Гявжт-л+1) 





4[f(t+0)  f(£—0)] = в) 8, s=o ir, (6.4) 


o> R(k- A—1). 
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Hence the theorem : 


Theorem 5. Let f(t) be bounded іп #>0, and of bounded parietem ín the 
neighbourhood of the point t, : 


o 


and ot) = р | от мо (ђе irat 
0 


Then А1000) +f(t-O)] = 4 I. Иа s = он, 


o> R(k+A-1), 
where Ф (в) is the Mellin transform of 9(p), 
provided that 2 f((t7 77! and q(t)tv7! belong to L(0, ee), 
and R(s -m — A) - 1. | 


Further if we take f(t) to be continuous in t > 0, we obtain the inversion formula т 
the form 


Too 
201 пре =) у 


provided the other conditions of the theorem are satisfied, 


We may verify the inversion formula by considering an example here. 


We know that 


Г.(иЕт-Л+И pee | ie - 176 -A- 
PETAT ip ue? = | eX" ТУ ут (pt)(pt) ^ t^ 
5 I(u—k—A-c1) А 


p-k-At+l; 
x F( а t (6.8) 
pim-A-i; 


Then according to the inversion formula, established above, we have 
Г =й) ОИ 2) 
Iuxm-A-1) uim-A-c1; 
pof pa-keig t) 
== ea tds, 6.7 
an A Гата) а дш Tor 


R(g) > а > R(k+a-1), 
t Ра; ty; -ЭЕ Fal; £t y, 8–7; –0). - 
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which may be shown to be true by evaluating the integral by swinging round the contour 
во as to enclose the poles on the right of the contour). 


DEPARTMENT OP MATHEMATIOS, 
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BOOK REVIEW 


S. Mitra and 6. К. Dutt :—A Tezibook of the Integral Calculus, №. Heffer 
& Sons Ltd., Cambridge, Great Britain, 1953, xiv + 89; pp. 80s. 


The book is supposed to cover the B. Sc. Syllabus of the Indian Universities. The 
authors wished to deal with definite integrals first but the position in India, as they 
remark, prevented them from doing so, We do not know whether the authors studied 
in the Calcutta University, Even as early as the year 1916, the only text book followed 
here was Integral Calculus for Beginners by J. Edwards (MaeMillan & Co. Ltd ) which 
begins with the Integral as the limit of a sum. | 


In establishing the uniquaness of the Indefinite Integrals (р. 2) the authors have 
referred to § 58 of their text of the Differential Calculus, where the discussion given is 
not sufficient. 


The authors want to develop Integral Calculus in a logical manner, unfortunately 
they assume the existence of area without proceeding to explain tha difficulties involved 
in defining area enclosed by a closed curve (p. 72). 


In defining surface integrals, no mention is made of the surface of Schwarz 
(p. 268). 


The authors have used Fourier series in partial differential equations before 
treating Fourier series. 


The way the development of Fourier Series has been done is too difficult to be 
followed by an average student, 


However, there ig no major flaw in the treatment of the subject. 
Lastly the price of the book is beyond the reach of the most of the Indian students, 
The printing and the binding of the book are good, 


В. Dorr, 


CALCUTTA MATHEMATICAL SOCIETY. 


Report of the Council on the state of affairs of the Society for the year 1960, 
placed in the Annual General Meeting of the Society. 


The Council of the Caleutta Mathematical Society has the pleasure to submit the 
following report on the state of affairs of the Socisty for the year 1960 as required by the 
provisions of Rule 25 of the Society’s Constitution. 

The Council. The council of the Society for the year 1960 consisting of officers 
end other members elected at the last Annual General meeting with Assistant Secretary 
and Editorial Secretary was constituted as follows: 

President 


8. N. Bose, 


Vice-Presidents ` 
R. N. Sen, | | У, У. Narlikar, 


P. C. Mahalanobis, B. № Prasad, 
N. R. Sen, 


Treasurer 

5. О, Ghosh, 
Secretary 

P. P. Chattarji, 


Editorial Secretary 
P. К. Ghosh. 


Other Members of the Council 


. P. K. Bose, Saktikanta Chakraborty, (Asst. Secretary.) 
A. C, Banerji, . B. C. Chatterji, 
Sanat Kumar Basu, А. С. Choudhury, 
Ram Behari, S. С. Das Gupta, 
М, б. Bhabde, B. R. Seth, 
M. C, Chaki, B. S, Ray, 


General. In this year, in addition to usual activities, the Society has undertaken 
some new ventures viz. the renovation of the library of the Society, the organisation of 
Symposia on recent advances in different branches of Mathematics, the 
orgsnisation of Sjmposia on teaching of Mathematics in Schools to devise ways and 
means for improving the standard of teaching of mathematics in India, the arrangement 
of publication ‘of the short report of the progress of the researches in Mathematics in 
India in the last fifty years: The Society could not undertake a few more important 
works viz. the arrangement of publication of the Proceedings of the Symposia on Fluid 
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Mechanics, the organisation of the Symposia on teaching of Mathematics in the Graduate 
and Post-Graduate level, the increase in the number of subscribed journals in the Society. 
A short report of the activities of the Council in different heads is given below: 


Symposia. (1) Symposium on Teaching of Mathematics in Schools, On the 28th., 
20th. and 30th. of March, 1960 a symposium on Teaching on Mathematics in the schools was 
held. Prof. 5 K. Chakraborty presided over the General Symposium. Principal В. K.Basu, 
Prof, D. №. Ray and Proi. 8. К. Chakraborty presided over discussions for Higher Secondary 
stage, Middle stage and Primary stage respectively, Many eminent Mathematicians and 
distinguished teachers from Colleges and Schools participated in the symposium. A 
short report of the proceedings of the Symposium will be published soon, 

(П). Symposium on Teaching of Fluid Mechanics т Post Graduate and Research 
level. On the 28rd. of September, 1960 a Symposium on teaching of Fluid Mechanics 
in Post-Graduate and Research level was held. Prof. N. В. Sen presided over the 
Symposium, Prof. N, L. Ghosh, Dr. J. N. Kapur and Dr. A. K. Roy made important 
suggestions. ‘The report is being made ready for publication, 

(Il) Symposium on Fluid Mechanics. On the 22nd., 28rd.and 24th. of September 
1900 a Symposium on Fluid Mechanics was held. Prof. М. В. Sen presided over the 
Symposium, Ж; 


Speaker Topic 
N. R. Sen bes On the study of Fluid Dynamics. 
N. L. Ghosh ive On some recent investigations regarding the 
Drag and lift on aerofoils. 
J. N. Kapur ө « Some aspects on superposability in Fluid 
Dynamies. 
А. K. Ray m Study of some aspects of the Flow 


characteristics in shock wave boundary layer 
interaction problems. 


8. D. Nigam & К, 5. Sinha Some problems in the theory of rotating 
Fluids. 

М. К. Jain T On extremal point collocation method for 
Fluid flow problems. 

R. 5. Nanda се: Viscous flows with suction or injection. 

R. N. Bhattacharya Bed ейесь on wave resistance of а ship 


moving in a circle. 


A. R, Sen zur Stationary phase and missi g waves. 

J. М. Kapur & В, К. Jain Parallel flow in an annular channel in 
hydromagnetics. 

К. М. Ghosh о The early period decay in case of any general 


type Turbulence with specific hypothesis 
borrowed from normal distributions, 
‘Effort is being made to publish the proceedings, 
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Pablications. During ths year under review the Society has published three issues 
of the Bulletin viz. Vol, öl, No. 4 ani Vol 82 Nos. 1 & 2. Тһе Council takes the 
opportunity to record here its indebtedness to the authorities of the Calcutta University 
for printing the Bulletin free of charge and to the Officers and members of the Staff 
of the University Press for their valued services, It may, however, be mentioned that 


the Soeiety bears the cost of paper and blocks and the matter of Bulletin upto May, 
1960 was composed by the Society’s compositors, 


The Council 18 glad to announce that the printing of the short report of the progress 
of researches ш Mathematics in India during tue last fifty year is made possible by 
benevolent patronage of the Government of India and the Government of West Bengal. 

The Council 13 glad to announce that tho printing of the Proceedings of the 
Sympsoium on Teaching of Mathematics in Schools is in progress. The report on 


the Teaching ot Fluid Mechanics in Post-Graduate and Research level is also bemg 
made ready. 


The Council is glad to announce that the printing of the first part of the Golden 
Jubilee Commemoration Volume has been completed and will soon come out of the press. 
The printing of the second part is in progress. 


Exchange of Publication. The transmission of the Society’s publications to various 
Countries of the World has been carried on regularly during the year. The Society has 
received one hundred and thirtytwo journals in exchange this year. 


Library. The Society has purchased only five books and subscribed sixteen 
journals this year. Some useful journals lying unbcund for several years have been bound. 
The Council feels that the Society should purchase more books and journals if the library 
is to serve the need of the active research workers to the desired extent, which, however, 
is 006 possible for shortage of funds. The rack has been expanded and a book case 
ћав been purchased for better upkeeping of Journals. 


Meetings Daring 1960. The Council has mot nine times this year and there have 
been six Gen-ral meetings of members in which sixteen original papers have been 
discussed and nine original papers have been taken as read. 


Membership, The Council desires to report that during the year under review 
twenty-one new names have been added to the list of Society’s members. Two gentlemen 
have become life membere of the Society. 


Delegations to Learned Societies. In reply to invitations from the Indian Science 
Congress Association the Council nominated the following members to represent the 
Society at their 48th, Session at Roorkee. 


(1) Principal 8 К, Basu, M.A., Ph.D., Presidency College, Calcutta-12. 
(2 Dr. M. C. Chaki, M.Sc., D.Phil., Calcutta University, Calcutta, 
(8) Prot. М. L. Ghosh, М.82., D.Phil., Presidency College, Calcutta, 
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(4) Prof, М. Ray, D.8c., Agra College, Agra. 

(B) Prof. B. B. Sen, D.Sc., Jadavpur University, Jadavpur. 

(Ó Prof, В. В. Seth, D.Sc., Indian Institute of Technology, Kharagpur. 

At the 20th, Session of the Indian Mathematical Society at Chandigarh, the Council 
nominated the following members to represent the Bocieby. 

(1) Dr, A. C. Choudhury, M.8c., D.Phil , Calcutta University, Caloutta, 

(2) Dr. G. Bandyopadhyay, M.Sc., D.Phil., Indian Institute of Technology, 

Kharagpur. 

The Council nominated the following members to represent the Society at the 

Silver Jubilee Celebrations of the National Institute of Sciences of India at Delhi. 


(1) Prof, 8. М. Bose, Е.В 5 , National Professor. 
(2) Principal, S. K. Basu, Ph D., Presidency College, Caleutta-12. 
(8) Prof. Ram Behari, D 8c , Delhi University, Delhi. 
(4) Prof. P. L- Bhatnagar, D.Sc., Indian Institute of Sc., Bangalore. 
(5) Prof, В. К. Chakraborty, D.Sc. B.E. College, Howrah. 
(6) Prof. B. М. Prasad, D.Sc., Allahabad University, Allahabad. 
(7) Prot, B. B. Sen, D.Sc., Jadavpur University, Caloutta-82, 
(8) Prof. М. В. беп, D.8c., Ph.D., Calcutta. 
(9) Prof. R. N. Sen, Ph D., Calcutta University, Calcutta. 
(10) Prof. А. C. Banerji, Ex-Vice-Chancellor, Allahabad University. 


Finance, The Annual Accounts of the Society have been presented to the Council 
in the standardized form by the auditors Bri P. L. Ganguli and Dr. М. Mitra. The Council 
offers them its sincere thanks for their honorary services. The Council also takes this 
opportunity io make a few comments on the audited statement of accounts durin, the 
year under review. A glance at the receipts and disbursements reveals some improvement 
in the financial position of the Society, thanks to the grant given by the Government of 
India which was not forthcoming for the last two years.. 

[t is, however, obvious that unless the Society receives more grants-in-aid from 
the Government and donations from the members and the public, it will be difficult 
to increage the activities of the Society. 

The Society received the follow ing grants for the year under review :— 


(i) Government of India—Hs 4,500,- 
(i) Government of W. Bengal—Rs, 2,000/- 
(iii) National Institute of Sciences of India—Rs. 1,000/- 
(iv) Govt of West Bengal (Sympos'um on Teaching)—Rs. 1,009/- 

The Council offers its grateful thanks to the Government of India, the Government 
of West Bengal and the National Institute of Sc.ences of India for these grants, 

Obituary. This year the Society Icet two of 165 esteemed members Dr, Н. M. 
Sengupta, Reader, Department of Pure Mathematics, Calcutta University and Dr. N, G. 
Shabde, Chairman, §.8.C.E. Board, М.Р. The Council takes the opportunity 
to offer condolence to the bereaved families, 
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CALCUTTA MATHEMATICAL SOCIETY 


93, Acharya Prafulla Chandra Road, 
Caleutta-9 (India) 
Dated, February 27, 1961. 


É NOTICE 


The Calcutta Mathematical Society has entered into a reciprocity agreement with 
the American Mathematical Society. The conditions of the agreement are given 
below : 


Under such an agreement, any member of the Calcutta Mathematical Society 
may join the American Matbematical Society by submitting an application for membership 
on the usual form. He would continue to pay the usual dues to the Calcutta 
Mathemeties! Society but would pay half dues only to the American Mathematical 
Society. At present, the full rate is $ 14,00 and reciprocating members pay annual 
dues of $ 7.00. Members of the Society under one of the reciprocity agreements enjoy 
all the usual privileges of membership such as receiving the Notices, the Proceedings, 
and the Bulletin and presenting papers at meetings of the Society, purchasing the 
publications of the Society and certain other publications at reduced rates and subseri- 
bing to Mathematical Reviews at the members’ rate of $ 16.00 per year. They are also 
permitted to substitute the Transactions for the Proceedings upon payment of a premium 
to dues of $ 10.00 in 4-volume years. They are also allowed to substitute Mathematical 
Reviews for the Proceedings by paying a premium of $ 10.00 to the dues, 

On the other hand, under such a reciprocity agreement, any member of the 
American Mathematical Society would be sdmitted to membership зо the Calcutta 
Mathematica! Society at his reqest and would pay half dues to the Calcutta Mathematical 
Society. He would then enjoy all the usual privileges of membership of the Calcutta 
Mathematical Society. 

There is one exception to the above schedule of dues. Ifa member of our Society 
became a member of the American Mathematical Society through such a reciprocity 
agreement and then visited the U.S.A., he would pay full dues to the American 
Mathematica! Society for the time spent inthe U.S.A. A similar arrangement would 
take effect for a member of the American Mathematical Society visiting our country. 

Members of our Society who are desirous of bemg members of the American 
Mathematical Society under reciprocity agreement should write to 

The Executive Director, 
American Mathematical Society 
190, Hope Street, 

Providence 6, Bhode Island, 
U.8.A. 


Application forms for membership of the American Mathematical Society may be 


had form the office of our Society. | 
P, P. СНАТТАВЛ, 
Becretary, 


Calcutta Mathematical Society. 


All correspondence with the Society, subscriptions to the Bulletin, admission fees 
and annual: contributions of members are to be sent to the Secretary, Caloutta 
Mathematical Society, 98, Acharya Prafulla Chandra Road, Caloutta-9, 


Papers intended for publication in the Bulletin of the Society, and all Editorial 


Correspondence, should be addressed to the Editorial Secretary, Caloutta Mathematical 
Society, 


The publications of the Caloutta Mathematical Society’ may be purchased direct 
from the Society’s office, or from its agente—Messrs. Bowes & Bowes, Booksellers and 
Publishers, 1, Trinity Street, Cambridge, England. 


NOTICE TO AUTHORS 


The manuscript of each paper communicated for publication in this Journal should 
be legibly written (preferably type-written) on one side of the paper and should be accom- 
panied by a short abstract of the paper. 


References to literature in the text should be given, whenever possible, in chrono- 
logical order, only the names of the authors and years of publication, in brackets, being 
given. They should be cited in full at the end of paper, the author’s name following form, 
vis., name or names of authors; year of publication; name of the journal! (abbreviation) ; 


number of volume; and lastly, the page number, The following would bea useful 
illustration : 


e 


Wilson, B, M. (1922), Proc, Lond, Math. Soo. (2), 21, 285. 


Authors of papers printed in the Bulletin are entitled to receive, free of cost, 50 
separate copies of their communications, They can, however, by previous notice to the 
Secretary, ascertain whether it will be possible to obtain more nopies even on payment 
of the usual charges. 
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ON LATTICOIDS 


By 


А. C. Сноорноку, Calcutta 
(Recerved—June 1, 1960) 


{. Introduction. An upper semilattice is an additive commutative semigroup in 
which every element is idempotent, It can be therefore expected that many properties 
of an abelian group can easily be generalised to properties of an upper semilattice. One 
such: property of an abelian grcup has already’ been generalised (Choudhury). It 
has been proved there that the endomappings of an upper semilaitice form an upper 
m-semilattice. This corresponds to the wellknown result that the endomappings of 
an rdditive abelian group form a ring This facb in the theory of an abelian group 
produces the notion of a module and all the results associated with it. It is natural 
therefore to expect that the m-semilattice of the endomappings of U will generate the 
corresponding notions of a latticoid and many beautiful resulta will follow from this 
notion, The object of this paper is to build up the theory ofa latticoid. An upper 
semilattice U which admite a scalar multiplication with the elements of an m-semilattice 
Q has been defined to be а lattieo:d. The condition that а subset may be a sublatticoid 
hus been found and a dependence has been developed. Further Lhe isotone mappings 
between a right lattieoid U and а two sided Jatticoid V have been proved to form a left 
Ја соја and has been denoted вв Iso(U,V). A commutative m-semilattice Q can be 
considered ав a two sided latticoid and consequently U* = Iso(U,Q) 18 а left Jatticoid 
when U isa night latticoid. U* may be called the dual of U. The connection between 
U and U* has been determined, Further 16 has been shown that sn isotone mapping f 
between two latticoids U,V over О induces ап isotone mapping ff between U* and V*. 
Lastly these results have used to obtain a representation theory for an upper 
m-semilattice. | 

9. Lattioold, Ар upper semilnttice U which admits the elements of m-semilattice 
Q as а right scalar multiplier ів culled a right latticoid 1: the scalar multiplier satisfes the 


laws :— Е 
(i) aae U Вбр ав О, ава, бео. 


(ii) аја+ђ)=аа +ађ, (2 +Вја=аф + Ba 
(ut) (aa;b = alab) 
(iv) al = а when 1 is the topmost element of О, 
If the scalar multiplication is on the left 1,8. if аа is defined, then the system will be 
called a left lattico'd. 
One can easily construct an example of a latticoid. For, from the elements of Q, 
consiruct elements. . 


and impo:e on the elements а the rules :— 
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(i) «a = (аа, а,а,......... „ала) 
(ü) а+В (ај, Dy, а, + b,,......... „ан + bn) 
when а = (4;,03....... ,ün) 
B= (b,,b,,......... bak 


Then the set U of elements о form a right latticoid. For, as a, b; are the elemenl: 
of an m-semilattice, the sum а; + b; are associative, commutative and idempotent. Hence 
U is a semilattice. The laws (1) to (iv) can be trivially verified. Hence U is a right 
latticoid, 

Let U be aright latticoid over the m semilattice О. Then the mappings D,:+—>2a 
is an endomorphism of U. For Ра +В) = (2+ 8)0 = aa-- ga = Dala) + БВ). 

D, will be called the dilation of the semilattice U. The condition (аајђ = alab) can 
be written as Da(2) = БьОо(«). Hence DDa = Da. This gives a rule of multiplication 
for the dilation D,. 

Similarly the mapping #,:а Э аа where а isa fixed el ment of О, is а homomor- 
phism of Q into U. For Ё, (а) = аа and t,la +b) = «(a-- b) = aa ab = 1.(а) + t. (b). 

Sublatticoids of a latticoid U over Q are those subsets which are again latticoids 
over Q. Then one can prove the theorem: 

A subset U,, is a sublatticoid if and only if 

(i) а+880,, when а,Вє U, 
(ii) аа e U,, when аа U,,a6Q, 

Proof. If U, is a latticoid, then these conditions must hold by the definition. 
Conversely, lst these conditions hold. Then if a,&j4 € U, 2+8=86+а and (+В) +1 = 
at+(B+y) since z,8,Yy є U. Also аво eU, ata = a. Hence U, 15 а semilattice. Further 
as a € U, and a,b в Q, 

а(а +b) = aa+ab 

Ав ae U,, a,b e Q, the condition (ii) gives that «a,ab e U, and tl e condition (1) gives 
that «a об є U,. 

As a,b в Q and О is an m-semilattice a+b є О and hence a(a +b) = за+ 26 i18 true 
in U, Та a similiar manner, other laws can be proved. Thus О, 1з a latticoid. Hence 
the theorem is proved. 


Theorem :—Sublatticoids of U form a lattice with U, N О, as the lattice product 
and U,+U, as the lattice sum, 


Proof. In the Boolean algebra of subsets of U, determine a topology X — X where 


X is any веб and X is the smallest sublatticod containing X. This smallest latticoid 
exists аз the meet of two sudlatticoids is а sublatticoid. The closed elements form a 


lattice in which the product and the sum are U, (1 U, and U,+U, = U, U U,. 
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The elements of U, + О, are 2+8 where a e 0,660, 


Proof. Asa,8 U, U О, +в 0, +0, = О, ЏИ. The seb of. all clements «+8 
where ae U,, Be U, form а sublatticoid. 

For, (a, +81) + (а, +82) = (a, +.) + (8, B5! = аз + Ва where a, =a, +a, and Bs =$; + Ва. 
As U,,U, are sublatticoids and «ла є U,, 8,836 Ол a, a, € U,, 8, 0,6 Оз. Further, as 
U,,U, are subl«tticoids and a eU,, Be Us, а’ = aa e U,, В = Bae U,. Hence 


(« +вја = са +ва = «^g 
where 4 є U,, B'e U,. Thus the set of elements «+В forms a s iblatticoid. But by 
definition U,+U, = 0, U С, is the smallest sublatticoid containing U, and U,, U, + О, is 
identical with the sublatticoid formed by the elements а +В. 
The representation of ihe elements of U, + О, is not unique, unless U, N U, = 0. 


For, if U, П О, + 0, thereis ап element ye U, N Ua. Tet ze U,, Вв U, be any 
two elements and put а“ = а+у, В = 8 y. Then “+В = («+ +(8+ у) = «+В+у = 
а +8' a 48, Thus the representation is not unique. 


The sum U,+U, when U, U,=0 is called a direct edt and is denoted as 
U, ФЦ.. 

3. Dependence. In a right Jatticoid U, an element а is called dependent on the 
elements 


1,03; че Lg 


if a = ай Фада + апад, n being finite. Let I be a setofandices. Let aielbea 
seb of elements of U. An clement а is said to be finitely dependent ор S if а is dependent 
on a finite number of elements of 9. 


Theorem. The clement a is dependent on azta... ‚ая tf ай = аза) Наза:+...... + алан 
and а is а reversible element of О, 7 


Proof. Ав а ів reversible, a^! exist and aa^! —1 where 1 is the top element of Q 


Then if «a = aa, + 2362 t ...... +andy, one gets by multiplying with а“7, (аа‘а-? = 
(ла ја“! + faga,)a~?...... (па ја“ 7. 


But (ваја ! = ala a7?) = al = а by axiom IV. 
Putting aja? = by, 

а, == ab. Бара +...... +anbn where 5,6 О. 
Hence a is dependent оп a,,a,...... n. 
The elements dependent on a,,25,...... an form a sublatticoid. 
Proof. Let z,«' be two elements dependent on a,,2,,...... jn. 
Then а = ауйу+@@,+......... + алал 


a! = 2121 Када -+...... + andr’ 
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Hece ate’ = ala, а,’ as (а 4 a,)...... tan(@p tan’). Also if a is any element 
of Q. 


аа = (2,a,)a+ (ха ја +... + (ааа. 
= a(a,a) + а (ада) +... + а, (ада). 


Hence both а +а/ and aa whereue Q are dependent оп «,,2, ‚з. Thus the ploposi-_ 
tion bold. 


If а is not dependent on 2,,25,...,95, then а is called independent of z,,x,, „ха. Thus 
if о is independent of а;,а,,.. бу, 


а а,а, + аза +... + ада for all аџ,а,,.. ‚ав, 
The set of independent elc ments of a 1а соја їв called а basis of the ћи вала, 


The smallest latticoxd contaming a given set S of elements of О is called a latticoid 
generated by 8. | 


Theorem. Lei S be a subset of а Та соја U. In order that an element ae U may 


belong to the aublatticoid V, generated by S, it is necessary and sufficient that a ia finitely 
dependent on 8. . 


Proof. Let « be dependent on the subset 8 of elements. 


Then, Ту the definition of dependence, 


€ = аа, ада, +...... +anQn, о 6 8. ETET 
As a; e S, the set A; of all аза form a sublatticoid of U. б 
For (via)b = одађ) e А; and aatob = ua +b) e Aj. 
The get, алб t. ads... auda € Ay+A,+...tAne V. 
Hence ae V. 


Thus the condition is sufficient. 
The condition is also necessery, 


For, let V be generated by S, Then if а, %,--.%, are elements of 8,...2,4, аза, + 
sys +4,d,¢V. These elements fcrm a sublattiioid of V. As by difinition, V is the 
smallest eublatticoid generated by S, the elements of V are expressible as у 

2418, Ноа. +... Hanin 
where а, 2,,..., ал із а fiaitg number of elemen's of 8. Hence the elements of V 


arc finitely dependent on S. 


8, Isotone mappings. A homomorphism ora linear mapping of a latticoid U 
into another Jatticoid V, both of which are over the same m-semilattice Q, is a 
homomorphism of the uj per semilattice U into the upper semilattice V such that 


Каа) = Каја for any а є U ав О. Thus the linear mapping is such that 


f(a+B) = f(a) +/(8) 
Каа) = (ада. 
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This is a special case of the isotone mapping of the latticoid. Ап isotone mapping 
f of the latticoid U into the latticoid V 15 а mapping f which keeps the inclusion invariant 
and is such that /(ха) = f(a)a. Similarly one can define antitone mapping of the Ја соја. 
The set of isotone mappings of а latticoid U into another latticoid V will be denoted as 
Іво (U,V). The set of all linear mappings UV will be denoted as Hom (U, V). Then 
Hom (U, V) S Іво (U, У) 


Theorem. Iso (U,V) of two right latticoids U,V, is an upper semilattice, 
Hom (0, V) is а sub-semilattice of Іво (U, V). 
Proof. Define the sum of two isotone mappings f, g, by the rule 
| (|+а)2 = f(x) -g(») where ze U. 
The sum {+g во define is isotone ; for, let а, В be two elements of U such that а €x В. 
Then (| + да = fa +ga 
(/+9)В = 18 +98 
As f, g, are isotone, [(2) Є f(8) and g(a) © 4(8) Hence ff) + g(x) са f(8) +9(8). 
(f +9)а = (+98. | 
Also (f+g)faa) = Цаа) + д(аа) = Каја + g(a'a = [f(a) + (а) Ја 


= [(f+g)z]a. 
Hence [++ «Iso (U, У). 
It follows from the definition of the sum, that 
ї+9 = 9+], ftg)t+h =ft+(gth) [+ =f. 5 


Thus Іво (0, V) is an upper semilattice, . As linear maps are isofone, Hom (И, F)isa 
subsemilattice of Iso (U, V). 


The mapping I, defined by the condition that 


la=leV for allaeU 
is the topmost element of Iso (0, У). 
For, (1 +h) = Та + ва = 1 +ћа = 1 = la 
Непсе 1+ћ = 11.е- А61 
, Theorem. so (U, V) is an upper m-semilatlice wilh the identity mapping i PM 
_ the identity. : 
| Proof. In the above theorem, Iso (U, V) has been proved to be an upper semilattice. 
Now define the produet fg by the rule 
(јдја = figa). 


This product is associative, 


For [ffg)h]a = f(g(ho)) = [f(gh)]a · 
It is also right distributive. , 
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For [(f+g)h]a = (1 + (Ва) 
= (fg) when ha = а є U 
= ја + ga' 
= fha + ghz = (fh gh)a 
If iis the identity mapping, 
Е ia ag 
Then gia = gz. Hence gi = gq. 
Similarly ig = g. 


A two sided latticoid is one which is both a left Jatticoid as well as a right labticoid 


and in which (аа)Ь = a(ab) hold. 


Theorem. If Uisa right latticoid and V is a two sided latlicoid over an upper 


m.semilatlice Q, then Iso (U, V) can be made a left latticoid. 


Proof. As V is a right latticoid, Iso (U, V) is an upper semilattice. 

Now defiine af, a є О, f є Iso (U, V) by the condition that (af)a = alfa) 

This is possible, as fe€V which is a two sided Ја соја. 

Then ај е1во (U, V). 

For (af)(sb) = a[f(ab)] = af (fa)b] (by the rule a(rb) = (ar)b in V) = [(ађајђ 
Also lt c&B,a,BeU. Then = f(a) S КВ) in V as f is isotone, 


So af(a) Є af(B), since the left scalar multiplication is isotone in V. 


Thus (af) са (af)B. 
Hence ај e Iso (U, V). 
It is easy to verify that alf +g) = ај+ад 
(a+b)f =af+bf 
if =f= fi. 


Thus in Iso (О, V), one can define a left scalar multiplication making it a left latticoid. 
In the above theorem, Q is a right distributive m-semilattice. Suppose that Q is 


both sided distiibutive, In this case, Q is a two sided latticoid provided the product is 
considered as the scalar products. For, 


(i) abeQ, when a,b¢Q 
(ii) a(b+c) = ab +ао, (b+c)a = ђа + са 
(Hi) (ab)c = a(ac) 


(iv) al=a 


Iso (U, QY is a left latticoid which will be denoted as U* i.e. U* = Iso (0,0). 


U* will be called the dual of U* 
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ра 


If U is a right latticoid, then U* has been proved to be a left 1а ома. Га a similar 


manner, 16 сап be proved that if U is a left latticoid, thon О" із n right litio. 
then U* is a right latticoid. Consequently if U іза right latucotd, U** = (U*)* за 
right latticoid U** will be called the bidual of U. 


We shall now study the relation between U and U**, Let з ba a fixs! elamsat cf 


U and f e U* = Iso(U, Q). Then f(x) є О. Denote by Pa; the mipping 


f>fla 
maps U* into Q. Then e.) = 10). 

Theorem. The mapping 9, : U*4Q is isotone on the left laltcoid U*. 
Proof. $.(а]) = (ара. = alfa) = ae.(f). 
Also 156 fag. ‘Then ftg =g. 
Hence Sf + 9) = s. (9) 
But Paf +g) = (f + да = ја tga = Falf) + v.(g). 
Thus $.(9} = Balf) + 9,(9) 
Hence e,(f) € Ф.(9) 

and 9. ів isotone on U*, 
Writing f(a) аз ја, one notes that 9,/f) = fa and 

hence је g implies ја са ga. 


Also | o © В implies ја & f8. 
Thus one gets that the product fa is isotone with each factor. 


As 9, 18 isotone on U*, Ф, is an element of U** The mapping: a3, is a ^natucl 


m .pping of О into U**, 


Theorem. The natural mapping у: 0-0 *% ів isotone. 


Proof. Let «,8 be two elements of U such that а < В. Suppose that f is ап isoton 


mapping of U, so that fe U. 


Then 16) e ГВ). 

As $. (f) = fle), Pelf) = FR), 
Pa (f) Є plf). 

Ав this is true for every feU*, $. & Ф. 


As Yla) = Фа, this gives thet (о) с (В). 
Also elf) = Каа) = Каја = e.(f)a = [v.a]f 
Hence | $a = 9,0 1.6. (аа) = у(аја, 


Thus the natural mapp по У is isotone, 
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Next let U,V be two latticoids over Q and f:U—V be an isotone mapping iê. 
fe Iso(U, V). 


7 v 
If ve V*, v maps У iuto Q and hence vf: U—>V—Q is an isotone mapping of U irto 
О and is an element of U*, Then the map 


ft. v-»vf 
is а шар V*—U* and will be denoted as |. It will be called the transpose of f. 
Theorem, The transpose ft is isotone i.e. ft e Iso (V*,U*) 


Proof. Let о, Sv, where v,,v, e V*. 


then us 01,0, f are all isotone, v,f S v,f. 
But as јој jtv) = vf. 
Heuce jv) € f'(v,). 
Further, as ає О, xe U, then 

[(va)j]2 = (ва) (уә) = (v(/3))a) = (ој) х а 

= (vf;aja. 

Непсе fiva) = за. 
Thus f! is isotone and fte Іво‘ У*, 0%). 


The mapping у: /->{ із а mapping of Iso (U,V) into Iso (V*, U *). 
Theorem. Ths mapping vy is isotone 1,8. у є Iso(I30(U,V), Iso,V,U1) 
Proof. As p:f—>/', (f) = F. 

Now let fg. Then vf < vg. i.e, f'(v) S gv) for all ve V*. 


Thus fag i.e. of) e plg) 

This proves thut у 1s isotone, 

Now let aeQ,veV*, 

Then (uft)v == v(af). 

If ae U, (о(ај))а = v(af(x)) = a 5] (з)) = a((vf)a). 
Hence vlaf) = a(vf). 

Thus (ао = af'(v) for all v e V*, 


Hence y(af) = ap(f) i.e. у admits the scalar multiplication, 
Hence У is isotone. 
This imply that y є Iso(Iso(U,V:, Iso(U*, У“), 


One ean evidently gen лаве this result. 


/ o 
Let U-V-—W ђе two isotone maps. Then f'g 9f їз а map Іво V,Wj—-Iso(U,W 
This map f° 18 isotone i,e, f^ в Isof1so(V,]V,, Iso(U,W)). 
Further the тар y fj 
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is isotone and hence y e Iso (Iso (U,V), Iso (Iso (V,W), Iso (U,W))). 


7 4 
Theorem. If U—YV—W ате isotone, then 
(off = fg’. 


Proof, Let we W* i.e, let ш be an isotone map У—>»а. 


As vov yw 
gw = v is an isotone map ГО. 
Also by definition, 
јој 
g*:w-wg 
(д) :w —w(gf) 
Now (gf)hw = wlgf) 


= (wg)f, as w g,f are mappings and hence the produc: is associative 


= g'twif. 
As g'(u) = wg =, 
80 aes (gfw = vf = fv) 
= fi(g'(w) 
= fig'(w). 


As this is true for every №, 
: (gf = е. 


All the above considerations have been made for the i otone mappings whieh form 
Iso(U,V). But these can also be made for linear mappings UV. 


5. Representation. Let E be а given upper m зе асе and V be a given upper 
вешт а ев, 


A homomorphism {:R->H(U) where E(U) is the upper m-semilattice of isotone 
mappings of О, 's called а representation- of В. 


Theorem, А representation јој R determines а left latticoid over R and conversely. 


Proof, Let f be a given mapping R>E'U). Then’byf,ra where а ів an element 
of E(U). Now define a scalar multiplication ra by the cule та. = ala) = f(r)s where r e E 
and «в U. Asa e E(U), a 18 an endomapping of U and consequently a(a) є U. Thus rz is в 
scalar multiplication of the elements of U with those of R. 


2—S057P— 2 
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Ава is an endomorphism of the upper semilattice О, 
ala, Фаза) = 4'2,) t ala) 

Непсе та, ta) = Ta, + та; 

Also, by the definition of the sum of tw» endomorphisiis 
(а, +а,)х = ааа + а 

Hence (пр +r a = 7,2 + 742 

Thus the scalar produci satisfies both distributive laws. 

Also, by tha definition of the product of two en Jomorphisms 
(аја ја = a,fa,3) 

Hence (тата) = т, (таа) 

Consequently U becomes а latticoid over R. 


Conversely, let U bs а lattieoid «ver an upper m-cemilattice В. ‘Then i£ ris n рун 
element of В, define a mapping a: U—U by the rule 


а :д— та, where ze U 

Then a(x) = та. 

Ав T(3 24) = Га + Тул, ala, + 2,) = a(2,) +а(а,). 

Thus a is an endomorphism of U. 

Also as (т +73) = та + rax, (а, a4)(3) = а (о) +а (а). 

Consequently the sum of the mappings а,.а, is given by thislaw. Lastly, as 
(пута = 7, (гаа), (8,24) = а (а,(2)). This gives the product of two mappings. 

With these definitions of sum and product, these endomorphisms form an upper 
m вета се H(U) The mapping f: RE (0, is а homomorphism. For, let f :r-»a є E(U;. 
Then f(r) = a and /(т)х = a'a) = r7 where x є U. Consequently /(г, +72)1 = (r, таја = 
тух yx = fira + (а) = (/(т,)х + га). 

f(r,rj)a = (гута = т, (т) = fir) f r)a) = (тј a. 
Thus f(r, * rj) = f(r.) fira); (гта) = flr fin) 
Consequently f is а hom omorpuisin of the upper m-semilatties В int» the upper 


m-semilattce E(U). Hence f is a 1epresentation of В through U. 
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NOTE ON THE STRESSES DUE TO A CENTRE OF 
DILATATION IN AN INFINITE SLAB 
WITH FIXED PLANE FACES 


Bv 


Summi Kana ВноммисЕ, Jadavpur 
(Ccmmunicated by Dr, Закн Ката Chakravarli— Received - Мау 19, 1962) 


Introduction. The effect of a centre of dilatation in a semi-infinite elastic solid 
with its plane face fixed has been determined by the author (1960) in a previous paper. 
The object of the present note is to find the stresses and displacements due to a centre 
of dilatation inside an infinite, isotropic elastic slab of finite and uniform thickness when 
parallel plane faces cf the slab are kept fixed. 


Solution. We take one of the plane faces of the slab as the plane z = 0 with the 
ступ on this face, Тһе axis of z ав drawn into the body. If there be a centre of 
dilatution at ође point (0, 0, о), we have the displacements given by (ef. 
Bhowmick, 1960) 








EDS ae 1 _ Belete) ] 
ET m zi RS R^ @—40)k,* | 
1 , 1 _ Galste) | | 
- =Р [- = + = = i 1 
d vfu = КАЙ BAR (1) 
ss | 
pai eu 92 " Geet | 
Ré Ку 8—АојВу (8—4c)H, Г | 
where Rj? = 13 +y?4+(z-c)’, 
l R? = +y? + (а + 0), 
P is a constant depending on the strength of the nucleus, and с is Poisson's ratio. 
It is evident that и = v = w = 0 whenz = 0, Writing г for (z? + y?)*, 
1 a 
we have BS f e~- J (ar)da. 2 >c 
1 
У (2) 


——— ag € 


i-[f eto (ar)da, 


0 
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Ј,(ат) being a Bessel’s function of the first kind and zero 


order. We can 
write u,v, W ав 


= 
| 


P| / e-e- Гу (ur) da — / в-««+е) [7 (a7)]da 2 | x67 «0*9 [Jo (ar)]da |, 
Д ex A On В ёт 


e 
|| 


© © © 
2l f e «0-09 [J (ar)]da— ] в-*ө+о [Jun Jd —2 | ase 7t D LT (us) da, 
0 ду 0 ду | ey 


0 





o со 
w= P[ J — ger) (ar)da + / as (nda 5 2. | взет faridah 
– дат 
0 0 0 


If weta Ње other face of the slab as given by а = 2c, we have on this face 


[шине = f sS tendis, | 
0 


гер = [ГОУ 2.6145, 
" / = 


| 

| (8) 
BE 

| 


Шә / Ria) lard, 
0 
where Qla) = P(e7*^ —8e-9-* — 4ace- 8), 
l R(a) = — Р[а6-* +а6- 8: + 4o! (8 — 4o)7 1079-7]. 
We are to find another system of displacements (U, V, W) which satisfy the 
equations of equilibrium, vanish at а = 0, and nullify the displacements (3) on я = 2c. 


In the absence of body forces, equations of équilibrium are (cf. Love, 1027, p. 133) 


р , 1 9 98 э) 
i up VWs = CS. 9.20): 
: УУ, 1—20 \-2’ ду да 
Let A =-2(1—20) E F 
9: 
such that | УЕ = 0, 
and УШУ, =v “2, $, „је 
Hence we can take U = gib. +49, " 
Ox e 
F 
yz E + Фа, ' 
OF 
l W 2s = 
| 8 5; + $, 


where 9,, Фа. Ф; are harmonic functions, 


NOTE ON TEE STRESSES DUE IO A CENTRE OF DILATION IN AN ETC. 
Let Ф, = ia а) sinh “217, (от) аз, 
: ад 
9, = | A(x)sinh р) lda, 
“о 


a 


фз = [ a Bia) sinh аз J (ата. 











70 
y . 
Since A О OY UB, ыч уд 

дт ду ағ Әз 

we shall have НЕЕ. f a?(B cosh az — А sinh а2) (от) 

, Oz 4с – 3 o 

whence ju cd J alB sinh az — А cosh az)d „(ат)“. 

4с —8 Д 


The above value of F gives 








Pac f da (wah gee A ай ЧӨ ЧЕТ f Aabe о enla 
4с —8 А : @х А Ox 

У = [ az(B sinh ag — А cosh on) [J,(at) da + f^ sinh 228.07, (от) ]48, 
4c – 8 я oy ^ ду 





х © 

РЯ sf a?g B cosh «z — А sinh а2)7 olar)da + || aB sinh sz Ј „(2т)42. 
== 

0 о 


On # = 2с we must have и+0 = 0, 9+7 = 0, wt W = 0. 


Therefore "5 (B sinh дас — А cosh 256) + A sinh 22c = - Q(z) 
с 





2 
za 5 (B cosh 220 — A sinh 2ac) + «В sinh 2ас = -- R(a). 
== А 


From tho above equations we geb 


_ (8— 4e) [ R(a)2c sinh 2ас + (а) {(8 — 4c) sinh 22¢ — ас cosh 2«c] ] 


Ala) = 
(а) [ (8—4c)* sinh? 2ac — 4a*c*] 





(8 —40-, [Q(a)2ex sinh 2ac — Е'аја " (8 — 40) sinh дас + Зао cosh2scl] | 


В(а) = 
(а) [(8 — 40" sinh? дас — 4а*с* | 
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Values of stress components can be calculated from the components of displacement 


ut+U,v+VandwtW. If #2, and Bey be the normal stresses deduced from (u,v,w) and 


(U,V,W) respectively, we obtain 
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Ja, oL 2EP [5 zx 
5776 (I-ej8-4e) В RE 


Where E? = 23 +y? tc, 


= E 
[224 ],=о= + (1+5 


аву aide, 
pcm 


Pu'ting о = .25 and c = 1, it is found that at т = 0 on the face 2 = 0. 


[22 Jomo = [ғ2, Жу = 2 28ЕР approximately. 
7=0 r9 
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SOME FORMULAE CONNECTING SELF-RECIPROCAL 
FUNCTIONS 


By 
8. Мавоор, Banaras, U P. 


(Communicated by Dr. Вгу Mohan - Recested- March 15, 1960) 


1. В Mohan (1984) gave certain theorems connecting different e’asses of fune'i па 
self-reciprocal in the Hanks! transf tm Iu thes note an attempt is being made 
to generalise tho-e theorems. In the end a few new resul's are obtained therefrom. 


Following Hardy and Titchmarsh we sha!l вау that a function is R,, if 16 is self- 
reciprocal for J, transforms and itis —J?,, if it is skew-reciprecal for J, transforms. Also, 
for Ry and R_, we shall write В, and Re respectively. 


We shall make use of the following known result (Hardy and Ti:elimarsh, 1980). 


А. necessary and sufficient condition that а function f'x) should be В, is that it 
should be of the form 


cio 


© = f see imita as, (1.1) 


where 0<с<1, and y(s) = (1—8). (1.2) 


9. One of the theorems of B. Mohan 15 
If f(z) тв Ra, the function 


UE ef 
giz = = ] o(iog = урау 
z y 
0 

is R, provided that | 


kd io 
Q(z, = RS / ГД + 4 +48)Г(@+{4уУ—{4в)'в)е®43 (x > 0) 
2лї p | (2.1) 
=0, (z < 0) 


where А is any positive integer and (в) satisfies (1.2) 
We generalise this theorem in the following form. 
Theorem I. Zf f(x) is of the form 


t+ ia 


Ка) = v j Гра 48) Г + Ду — 18) ГО + ФА + фа) Х (ата | i (2.2) 
" cio 


where X(s) satisfies (1,2), 
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then the function g(x) = i 3 (log = Vay 
. х 
0 


18 Hy, provided that Qiz) 18 given by (2.1). 
By (2.2) we have 


otir 


g(z) = А J Q(log = M | Nds. + фи + 48) (2+ dv — БЕДА +4ь)Х у tds, 


where A(s) = X(1— 8). PME provided the invcrsion of the order of integrati n 


is justified, we have 


2 


eio р. 
йа) = 2 | | октава та оганда | о(ов = уч, 
Е u t-te д : 
| МАЈЕ = 
= = | 2 Г-З + PDEA аге е (аа f. Q(u во- Оди, 
d 0 


8-15 


Now, using а form of Mellin’s Inversion Form Иа (Hardy, 1221) from (2.1) we have 
Ј 9s = Td arena +89), 
0 
where w(s) = (1—8). Changing s into 1—8, we get: 


/ c6-D7Q(zldr = Г(#+{4д-—4в)Г(} + {у + {з)ш(в). 


Hence 

etio 
gle) = эт D4 фа + do) d» - 99) HAA + Га dv de) 
x Г + фи —{43)у(8)ш'в)ж”*4в, 

ctio 

ш 2i T(3 +AA +48) (в) 7 da, 
Qnt ТТА 
where у(а) = Ги 4e) (3-- 4» —49)TQ t dv ага + iv — #8,Х (8) (3) 


As (8) satisfies tha equation (8) = y(1— в), it follows from (1.1) that g(z) ав R}. 


8. Corollary 1. When А = у, the above theorem reduces to the proposition 
established by B. Mohan referred to above. | 


Corollary 2. When и = v, the above corollary reduces further to the ИН 
simpler form : 
If f(x) is R,, the function 


glz) = УК = JfaMy 
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is R, provided that 


kt+in 


Q(z) = = | ез (в) ав (= > 0) | 
= 0, (z <0) 
where А(8) = A(1— в). 


Here we shall give a few interesting examples of this theorem. 


(i) Let 


мо = [ (Иов (3a) (да + Bm) pu pem а? || 


по Ру + тог (а В+ {а*+(в—{ е" * у a-Btl 'a'cr(s-4) 





m=O 


where Rib) > 0, | Ig | (8 -}), and — |b|« V [a* + (s —4)*] -- (s - 4). 
Then, by a formula given by Watson (1944), we get 


Ola) = ot Ја-в ад ба) (020) | 


xy-e-B у г 
= 0, (2420) | 
x x 
4 а log = Ју (bio =) 
Hence ga) | пы) AN 
0 


A—a-B 
(=) 


Аз а particular case, when А = ®+В8 we get 


a 


=z 
Я AL =) ( =) 
= = | - _ lg = Wag- = 
gu) а) yi а (а og y «+81 blog j f(y)dy. 





; т n (8—$)? b? с? 
(it, Let A(s) (bc) J «| Qbe | 


where R(s) > (6 +10 ic) and R(n) >-$. Then, with the help of another formula given 
by Wat-on (1944), we get 


Q(x) = еї, (а), (са) esos y 
- 0, @<0) j 
te т 
1 1 x „© 
Hence g(z) = | wil’ log РУДЕ log Y Voy. 


(üi) Let 


ма) = ОФА + 9) ФА + Ia ФА) ibe 
Мит TED (a  3)a7 237 *(s — је А 





=H) р 47-а 
Fen ens СИ 


ats 
8—2057P—2 
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where а < 8—4 and а +А із an even integer. Then, bya formula given by Titchmarsh 
(1927), we get 
Q(z) = a*-1ei Ki(az) (x > 0) | 


= 0, : (т <0) 


where K,(z) is a Bessel's function of imaginary argument of the second kind, 
Hence gix) = 1 HE (лав = y (а log £y. 
Putting а = 1 we get, as a particular case, 
1 [1 
ate) =] ураа = ade 


where A is an odd integer. 


(iv) Let 
we 25- Hb + фа + p) EYE bn P) p e [^ т = | 
| nl(n-g)(s—3)'*7** 3,n48 (8232 


where b and в— are both positive and n-- p-- & is an even positive integer. 


Then, by a formula given by-B. Moban (1941), we get 
Q(z) = 22-804 H, (bx) (ж > 0) 


у i0 0) 


where Н (с) is а Struve's function of order v. 
sf ant) "(tos =) 
(2) = gu = = dy. 
Hence g(a) at) yi DEC H, bles f(y)dy 
Putting p = 8 ав a particular case we get 
_ 1 е x 
g(a) = l H,( blog = titu, 


where n is zero or an even positive integer. 


2 D(oIT (0-3) on-2 n es OP, 0, 0+3. -an] 
8 ? 


(v) Let Мв) = лї (s — 4)?» ps Ва ' (3—4)? А 


where о is a positive integer 


(п) foralls—4 and b when p+2 <q 





(b) for | e< Eal | when р+і = 9 
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Then, by a formula given by Sinha (1944), we get 


Qiz) = xe-lett РРА Up – а | (x > 0) 


Ву --. Ва , 


| 
= 0, (2 < 0) J 


i 2р-1 | СР a 
Hence 4(2) = / («== aryl M id -(% log) оа. 
Е v Ва J 


4. One more theorem due to B. Mohan states that: 
If ра) is Ra, the funct on 


MES 
Я = > 1 
а = f eee ody 
is R,, provided that 
1 Кеа - \ 
O(n) = a J Фвр + ћи + 88) (d+ dv + 85)o(s)e7*da (x > 0) | 
Ол, TESA } (4.1) 
= 0, (x « 0) 


where ofa) = o(1—8) 
‘The generalisation of the above yields the following. 
Theorem II. If f(x) is a function of the form 


~ 


С too 
Ка) = zi оГ аи 35) P 3v 9T ФА — 3) (в) вав, 
f cio A 
where 0 «c < 1 and (s) = (1 в), then the function 


Ma 


ga = f (ов Кум, 
0 


is В, provided that Q(x) is given by (4.1). 
The proof of this is exactly similar to that of theorem I proved above. 


Corollary 4. When А = v the theorem 18 reducible to the one due lo B. Mohan 
given above. 


1 thank Dr. Brij Mohan for his help and guivance in the preparation of this note. 
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BENDING OF CERTAIN CLAMPED THIN 
CURVILINEAR ELASTIC PLATES 
NORMALLY LOADED OVER 
A CONCENTRIC CIRCLE 


By 


В. 8. Dnaurwarn, Ludhiana 
(Communicated by Dr. D. N. Мита, Received—July 11, 1960) 


In this paper complex variable methods are used to derive exact solutions in closed 
forms for the small deflexions of certain thin elastic plates subject to normal lcading over 
a concentric circle, the loading considered being symmetrical with respect to the centre 
of the circle. The isotropic plates considered are bounded by certain types of curvilinear 
edges along which the platos are rigidly clamped. The circular plate and a plate bounded 
by an inverse of an ellipse with respect to its centre are included as special cases. 


1. Introduction, Problems concerning the tiansverse flexure of thin clamped or 
simply supported elastic plates having various shapes and subject to various distributions 
of normal loading have been considered by several investigators, In a series of papers 
Bassali (1956a, 1956b) and Bassali and Dawoud (1956) discussed the bending of thin 
circular plates subject to different types of normal leading distributed over the area of an 
eccentric circle; methods of complex function theory being used to obtain the solution, 
Sen (1942) used curvilinear co-ordinates to solve the problems of bending of c'amped and 
uniformly loaded thin elastic plates bounded by quartic curves having the forms of the 
cardioid, the loop of а lemniscate of Bernou'li, the inverse of an ellipse with respect to 
its centre and inverse of an ellipse with respect to its focus (the elliptic limacon). The 
solution rppropriate to a concentrated load at any point of a clamped elliptic plate has 
been obtained in series form by Sengupta (1948, 1949). A simple method of deter- 
mining the deflexion of certain types of elastic plates, subject to transverse isolated 
loads ut certain points was given by Sen (1984) The same method was applied by Das 
(1950) to establish closed expressions for the deflexion at any point of (а) а clamped 
plate bounded by the inverse of an ellipse with respect to its сепіге and loaded by а 
concentrated force at its centre, (b) a clamped plate bounded by an elliptic limacon and 
subject to an isolated load at the focus. The solution for a clamped regular polygonal 
plate with n sides and subject to uniform loading distributed over the entire plate has 
been given ın a closed form by Stevenson (1943). When the same plate is under an 
isolated load at any point or loaded over a concentric circular area the solutions have 
been obtained by Dawoud (1950) or Bassali (1958) respectively. The solutions for 
cerlain clamped thin elastic slabs with curvilinear boundaries, cubject to transyerse 
concentrated forces or couples applied at arbitrary or specified points have been obtained 
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by Bassali (1959). Mukhelishvili’s method of dealing with biharmonic equation has been 
applied by Gray (1952) to the analysis of clamped thin elastic slabs under concentrated 
loads. 

In this paper complex variable methods are u od to find closed expressions for the 
complex potentials and deflexion nt any poiat cf а curvilinear plat, subject to normal 
loading overa concentric circular area. The load taken is symmetrical with respect to 
the centre of the circle. The plate, taken in the #-рјапе can be mapped on the region 
inside the unit circle in the ¢-plane by the mapping function 

в = об /(1+т +1), (c>0,|m|<1/p), 
where c, m and p are real constants, p being а positive integer. The tentative method of 
solution adopted here is to assume suitable forms for the complex potentials and we then 
show that they fit the conditions given, provided certam constants involved take on 
certain definite values, This method has been extensively used by Stevenson (1942, 
1043) and Bassali (1958, 1950), 

9. Formulation. Consider a clamped thin isotropic plate bounded by a closed 
curse Cy. Whun the middle plane of a normally loaded plate is assumed to be in the 
complex plane 2 = 2+1) its deflection w normal to the plane of the plate is governed by 
the bibarmenic equation 

4 — 
учо = 16 To Ple 0/0, | (2.1) 
where p(z, з) is the transverse load intensity and D = 2ER'/[B(1—v*) 18 the flexural 
rigidity in which E and v are Ycung's modulus and Poisson's ratio respectively, for the 
material of the plate, 2h being the thickness of the plate and bars used to denote conju- 
gate complex quantities. 


The general solution of (2.1) may be expressed in terms of {wo analytic functions 
(aj, V (21 and a particular integral W(z, 2) corresponding to the load intensity p(z, 2): 
w = вф(а) + 29(2) + х(2) + X(z) + W(e,a), (2.2) 


Х(а) = f (зав. 


W(a,z) may be taken ав zero if z lies in an unloaded region. 
Boundary conditions that the boundary C, is clamped may be tuken ав | 


w = 0, = = 0 along Cp, . (9.8) 


where д/дљ denotes differentiation along the outward drawn normal. It is easily seen 
that these conditions are equivalent to 


ш = 0, 90 = 0 along бу. (2.4) 


BENDING OF CERTAIN CLAMPED THIN CURVILINEAR ELASTIC PLATES ETC, 63 


Using (2.2), conditions (2.4) can be written as 


29(z) + #9(2) + Х(2) + х (2) + W(a,2) = 0 along Cp (2 5) 
zy (в) + 9(z) + (2) + ow = 0 along Ср (2.6) 


"where accents denote differentiation with respect to the argument. 


3. Clamped thin curvilinear plate normally loaded over a concentric circle. 
The transformation 





а= обје а, b= el, = (8.1) 


о:>0, |m|zalp, 


where c, m are real constants and p a positive integer, maps tho interior of a closed curve 
C, in the e-plane on the interior of a unit circle у in а compl-x ¢-piane. ‘The parametric 
equations of the closed curve C, are given by 





x cos 6 + m cos рб Ya sin 9 — т sin рб __ 
с e 


1 + та + дт cos в (2 +1)6' 


1 + т? + т cos (p+1)9’ (8.2) 


It is easily seen that Ср is the inverse with respect to the circle | z | = c of а closed 
curve Г, defined by 
z = с(сов 6 + т cos p6), у = c(sin 9— т sin рб). ` (8.3) 


The curve C, corresponding to р = 1 is the inverse of the ellipse z*/a* + y*/b* = 1 
with respect to a concentric circle of radius d, where in this case 


0 — 
g Те та D (0 << m x 1), 
_ 2d? _a—b 
du ru on - ve 


The resulting figure in tbe z-plane is called the lemniscate of Booth and is shown 
in Fig. lonp 70. When m is almost equal to unity it differs little from that produced by 
two touching circles of equal radii [Muskhelishvili (1949), p. 180; Sokolnikoff (1948), р 
176, 183]. The limiting case in which b > 0 and m 1 is усті consideration, In this case 
the plate in the s-plane is infinite and is clamped along the two semi infinite lines 
xz >> 4с, у = 0 andjrzx—te, у = 0 represented respectively by AB-(or AD) and А'В (or 
A'D) in Fig. 1, these lines being the only boundaries of the plate, 


When т = 1/p(p> 1) the curve Cp is the inverse with respect to а concentric 
circle of а star shaped bypocycloid Г, with p+1 cusps [Muskhelishvili (1949), p. 178]. 
When m= -2/р(р+1) the closed curve Гу approximates to a regular rectilinear polygon 
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with p+1 sides and the corresponding curve C, is therefore approximately bounded by 
p+1 circular ares each of which passes through the centre and two consecutive vertices 
of a regular rectilinear polygon with p+1 вез. The curves С; and Г, corresponding to 
p= 2,m=—},p = 8, т= апар = 7, m=—yy are shown in Figs. 2, 8 and 4 respec- 
tively. 

Let the plate be normally loaded over the area of a concentric circle C of radius 


с 


= Tal 


and clamped along the boundary. Let 1 and 2 refer to the two regions inside 


C and between C and C, respectively, Let the transverse load intensity p on the plate 
be given by 


т 
— 2-28 


p-pr"*- Ро(#2) 2, n > 2 over 1 
Pa = 0 over 2 (8.5) 


For n = 2 we have uniform loading over 1, 


The particular integrals W, and W, may be taken as 
УУ (2,2) = - po (22) * [n*n^ D, (n^ = +2), 


W.(a,2; = 0. (3.6) 
Let the complex potentials in the regions 1 and 2 be denoted by 9,(2) = e(t), 


Х, (#) = X,(t) and palz) = 9,(), %(z) = Х,(6) respectively. Since W.(z,z) = 0, the boun- 
dary conditions (2.5) and (2.6) are transformed to 


ш(с)о (с) + (с), (с)Х + Хо) +Х (о) = 0 (8.7) 


©(©)Ф,/ lo) + w'(m)e,(v) +Х„ (о) = 0 (3.8) 
where o = сі ів the value of t on ү. 
It has been shown by Dawoud (1954) that when the two regions are continuously 


and differently loaded, the necesrary and sufficient conditions for the physical continuity 
across C can be taken as 








1D, == D Ow, _ Ow, cw, _ Ow, 

1 2» g Cz , да: 825 , 
фи, _ Fw, ‚ 
Oro, бб, EC (ea 
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Now usirg (2 2) and (8 6) we find that the continuity condition (2 9; are satisfied if 





= > —— 41 87a" 
[22(2) + 29(2) +x(2)+xlz)] = ——„ (8.10) 
2 пп 
e DRESS 1 2 
[++] = 2% (8.11) 
2 пп z 
Ас, а. ed 1 
[p (e) +r] = 2 (8.12) 
leto] = 2 (8.18) 
3 2 
where x = Р./16я0, (8.14) 
P, = ?sp,a"[n. (8.15) 


Lt may be noted that Р, 13 tie otal load ^n the plate. 


As tentative solutions for tha.» continuity conditions we assu ос that 


[9:29] = (4 Blog 2 (8.16) 
2 


[х(2)] = E+F log А (8.17) 


where А, B, Е and F are real constants. Substituvng from (8.16? and (8.17) in (3.12) 
and (8 13) and equating the coefficients of te same pow rs of а on both sides we get 


Asi -1), Bes (3.18) 


Hence we have [ofa] = (2 —1+ log £) (8.19) 
Y 2 


Substituting frum (3.19) and (3 17) in (3 10) anl (3 11) and equating the coefficients 


of the same powers of 2 on both sides we get 


Е = ла + 1), Р = xa^", (8.20) 
n! 
where a? = naj. (8,21) 
А 1 1 4 (B 22) 
Hence we have [xa] = ва s +1 +106 = o 
2 


4 -2057P—2 
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Let us assume tentatively that 


9,02) = 946) = az (a, + ap! — log S), (8,23) 


Xs) = xt) = af f, e bui) - a logt], (8.24) 


$ 
where @, by, ар and by are real constants, 
Now from :8.19), (3.22), (8,28) and (3.24) we easily find that 


$1(2) = $.) = aa, etit t log 5 ), (327) 
5 
Х, (3) = ХИ = [trs eui +a? log |], (3.26) 
5 
where ao = а, +(1/п) -1, 


by = by *a*(0/n) +1{/е%, | 
(3.27) 
а' = ay t m((1/n; — 1), 
b^, = by. ) 
It is clear from (8.28; — (8.26) and (8.1) that 9,, x, and Ф, x, are regular in ther 
respective regions. 
Now we have to determine the real constants аб, b, ар and b, from the boundary 
conditions (8.7) and (3.8). Substituting from (3.1), (8.28) and (8.24) in (8 7) we find 
да, taplo? t! + 07771) + (b, + byo? * )(1-- mo ^97) + (by + bpa? 7) (1 + mst) = 9. 
Equating the coefficients of various powers of o to zeic we obtain 
а, +6, + mbp = 0 | 


(8.28) 
ay t b, x mb, = 0 


Similarly substituting from (3 1), '8.25) and (8.20) in (8 8) and simplifying we get 
[aot (р + 2)ay0?3 ! (1+ то?) 71 m(p + 1) (a, Каро" +) ја + (а) + ayo" 97) 1 тре?) 
ер 1)(bp— mb, (1 + mo^P-!)g? 
LE 72 , 
= [| + (14-2m3) 5 је“ +m eate [1 + (24 m?) a |= +m? а gm 
й с с ct 
which 15 satisfied if 
a 
с 





да, m'p-1)(b,— mb,) — mpa, = 1 + (1 + 2т7) 


2 


атра, ~ (р + 1)(byp— mb.) —(p+2)a, = "тј + (2 + m!) 21 ] (3.29) 


| 
| 
| 
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The five equations given by (8.28) and (8.20) are satisfied for the following values 


of aa. by, ay and by: 
a, = [L—m? + + (р + 19m? тА, 
bo =—[L+k{l+(2p+ И? Г, 


— аа 





(3.80) 
ар = km, 
b, = m|1-k(1-m*)]]/A, 
a’? 
where k= ст, А = 211 +рт?). (8.81) 
Now from (3 27) and (8.80) we have 
а’, = : — [14 (2b + 1)m* + + 2(p 1 1,m? — ту ГА, 
, k 3 
b's = —-—[I-k(1—m)]/A, 
n (3.82) 


| А 
‘p= m itk), 
а р „(1 ) 


ы ee .—--——— 


| 


b’p = m[1—k(L—m7)]/A. 

The bcundüry ecnditicns along Ср aud the continuity conditions across C are thus 
satisfied and the complex potentials are now completely determined by (3,28) – (3.26) 
and (3.80) – (8 82). 


The deflecticn at any point of the regi ns 1 and 2 is given by the following 
expressions : 





, ,b! , mk |, 7 т Ду 
w, = гај {а 0 Yo + тб’ ро"? + ( we a) cos (p+ yep (1° +a’) pn x] (8.83) 
w, = 2| nda bo, + тђро"Р + mi(1 - a)? cos (p+ 119} — (r? + а") log e}, (8.84) 
о 


where т? = 28. 


The deflection w, at the centre of the plate is given by 


w, = 2X,(0) = С а? log J 
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Substituting the values of z and b’, from (8.14) and (3.32) we bave 


5 Р.с? [ (i =] 3.85) 
w= 16:2 2k xw leg P ри? ^ . ( 


The bending mcmentg, twisting moments M,, Ms, M p and shearing forces. Q,, Qo cll 





рег unit Jengih, at any pomt cf the plate may be computed ! y using the standard 
formulae. We know that for unlcaded region of the plate 


М,+ M, = —8D'1* v)[v'.(e) + v^, (2)], (3.36) 


and since alcng the clamped bcundary vM, = Mg, we have 





[Ма ђе, = -ep| vx) , Фа ] (8.87) 
wi "uve 
Now rubstituting the values of ẹ,(¢) and e(t) from (8.28) and (8.1) ard using (3.38) 
and (8.14) we get 


P 
M =< = ву ар са фл ама ee 3 | У 
[Mel 2r[1+m*p’+2mp’ cos (p + 1)6] о. 
+ тір +8 —2a,— p'k(1 + m*)] cos (р + 1)6 + m*k cos (2p 4-2)6] (8 88) 
where р'=р+2. 


For m — 0, the previous results reduced to those already found by Вавва!! and 
Dawoud (1956) for a clamped circular plate normally and symmetrically loaded over а 
concentric circle. 


The solution corresponding to a concentrated transverse load P, acting at the centre 
of the plate, can be derived Гот the foregoing results by making the radius а of the 
loaded area approach zero. Now taking а = a’ = = 0, we get the following expressions 
for the complex potentials and deflexion for this case: 





1-—m? 
Z =з] 
ole) = doe] үст, ~2 log 
јаје — — 29. Тот 





2(1 + т?р) 1+ mute 


1 - 9? —m?(1— gf) 


ш = =" 


which agrees with the solution derived by Bassali (1959), 


In the limiting case in which p = 1 and m tends to unily the plate becomes infinite 
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and is bounded by and clamped a:ong the semi-infinite lines AB and A/D of 199.1. In 
this case taking m = p = 1 we get the following expressions for w, and w,: 








P, [ | 1 4k-n’ 1 
у E {5—-1+К+ 2 (з 
Pı 87р : n 4n/g? T 4 
+( EA LM, вов 20! - (r* 4.4") log е]. (3.40) 
n'o n 
P 5 1+4К © 
© = 80 L {к- д + Іо? 
«(1-3 ). сов 26} —(т°+ а?) log e]. (8.41) 
e 


In case the same plate is acted upon by a concentrated force at the centre of the 
plate, the expression for the deflection may be derived from (3.41) by taking a’ = k = 0. 
In this case we get 


w = a "[e-2) -log e], (3.42) 


which is a particular case of relation (12) of Dean (1953). 

When the plate is of the shape of a lemniscate of Bootn and ів acted upon by a 
concentrated force Р, at the centre, the expression for diflection can be obtained by 
taking p = 1 in the last relation (8.89) and hence putting p = 1 and substituting tbe value 
of а. from (8,14) in (8.89) we get, 





P, r? [ а 1 2 2 | 
= ——|1—т*—— + — (1 + m?) log о |. (8.43) 
lOgD 1+? ш e ен ) loge 
Now putting mig = ett) (0 — m « 1), 


in (8.1) we see that the transformation (8.1) takes the form 
, я = q вес (uti) 
where 9 = c[2m*, и = 0, тїр = e-*, (8.44) 
Now setting m = 678 (0 << 8 < oo) we have 
о = ef" с = 996-8, 
It is now evident that straight line v = В in (u,v) plane corresponds to the boundary 
y of the unit circle in the -plane and hence to the boundary C, of the lemniscate of Booth 
in the z-plane. 


With the above substitutions (8.48) reduces to 





2 от 2(8 —v) + (sinh 2v — sinh 28) /совћ 28 ] 


8.4 
87р cosh 2v + cos 2и (3.45) 


which coincides with the so'ulion derived by Das (1950) using the method of Sen (1934), 
on noting the difference in notation. 


TO - 
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p=2, m a-- 


Fig. 2 


LENDING OF CERTAIN OLAMPED THIN CURVILNER ELASTIC PLATES ETC. 71 





R. S. DHALIWAL 


Acknowledgement. Author, in conclusion, wih his sincere thanks to Dr, D. N 
Mitra, cf Indian Institute of Technology, Kharagpur, for his kind advice in the preparation 


of this paper. 


GURU NANAK RNGINPERING COLLRGE, 


LUDHIANA* 


References 


Вавзаћ, W.A. (1°66a), Proc. Camb, Phil Soc., 88, 784. 
Basaali, W A (1956b), — ————— ——, ва, 743. 
Basel, W.A, aod Dawoud, В.Н (1956), —— —— -, 82, 584. 
Bassali, W.A (1958). Вил Faculty of Sci , Al xandria, 2, 3. 
Bascal), W А. (1959), Proc. Camb Phil. Soc , 88, 121. 
Dae, 8 C. (1950), Bull Cal. Math Soc., 42, 89. 
Dawoud, К К (1920), Ph D Thesis, London 
Dawoud, ВН (1954), Proc Math Phys. Soc., Egypt., 6, 51 
Dean, W.R. (1953), Proc Camb, Phil, Soc., 49, 319. 
Саву, С.А M. (1952), J. Appl. Mech., 19, 490. 
Muskhelishvili, N.I. (1949), Some basic problems of the 
methemattcal theory of elasticity, 3rd Ed., M. вост. 
Sen, В. (1912). Phil, Mag. (т), 88, 294. 
8 n, В, (1984), Indian J. Phys. Maths., 8, 17. 
Sengupta, H. M. (1948), Bull. Cal Math, Soc , 40, 17. 
Ber.gupta, Н.М. (3940) —— ——————— ~, M, 168. 
Sokoln Кой, I.S. (1946), Mathematical theory of elasticity, 1st Ed., London. 
Stevenson, A.C (1948), Phil. Mag. 47), 34, 105. 
Stevenzon, А.С (1912), ———— ———, 83, 639. 


* Presently, Departmen: of Mathematics, In hin Institute of Technology, Eharugpur. 


ON THE MEANS OF AN ENTIRE FUNCTION AND 
ITS DERIVATIVES 
Br 
SATYA NARAIN Srivastava, Lucknow, 
‘Communicated by Dr. S. К. Bose —Received May 5, 1960) 


1. Let f(z) be an entire function of order о and lower order A and let 


2x 
barf) = af | fire) | ад (1.1) 
2r А 
2x 
prf) = | |f) (тет) | 846 (1.2) 
У 21 А 
1 qur 
ma (тј) = ff кв“) | *32^dxd6 (1.8) 
m AA 
T 2x 
must) = —L Ј f tfo) ачаа, (L4) 
miktl A 


where f) (2) is the n^ derivative of f(z). It is known that [ Rahman (1956), p. 193] 


` Dun 89 log log вв (т) _ 8 
тю Inf log hi А“ 


lim sup log log тз„(т,]) _ 0 
2B OR ат) NE ЈИ 
r+ inf рт À 
It is also known that [Srivastava (1959), p 280] 
log [4 нт} (2) || 
вир ихт,) _е 
те inf “Joop OS 


In this paper we would extend this resull to 0 «2 5 < 1 and to maur). We have also 
ohtained a few properties of msz(r). The results are given in the form of theorems. 


2 Theorem í. For every entire function (e), other than a polynomial, of order o 


and lower order А, 
og (rese И 
sup таз (т.ј) g 
rre inf log r Es 





where k > 5>1 
6 —2057 P—2 
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The proof of this theorem is based on the following two lemmas : 


For every entire function f(z), other than а polynomial, 


Lemma 1. 
8 
maj (rf) > mae n [Em , 
rlogr 
for all т> т, = (аа k 2821 


Proof, We have 


mul ft!) = Ls f | | fO (aote, |^84240 








f (xot?) — Гр — zee") 
-aff в Sew as a*dzd6 - 

1 (е) | - f(z x — ree) | р 
> dim | a ie = EON sagas (2.2) 


By Minkowski’s mequality n 1989, p. 384] 


ГГ | f(ze*) | -] f(z— z—xec%)) s] = [( (Г If(ze*; tae) - (Рива) 


0 


Hence, 
mga (rf) 22 Hina m LAT en vo) (Ге s) zi da] 
0 “0 
= im а 1 (ef iem nue) - - (2 | PEPEL 2°) ) Pa jJ Yaz] 


Again, using Minkowski's inequality, we obtain 


maj (7/09) = im png Tall (if Г Ком) |“ мав}! - | [| | (= ~ zee) Bat Чад | 
| 





1 
> lim abb Hm, л} - fan- +] в) тте) 


+: 


> lim 


«o0 НЕ єт 


[scit ima ir ] 


log m T 
og ть (т, then zince g(r) is а positive indefinitely increasing 


Now, take g(r) — 
log т 


function of rfor T> ту = то(ј), in facs log mai -a(r,f) із а convex function of log 7, and > 


we have 
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5 


i-1 

$ 

ma, (г, > Е {ты б mti | 
5 


= [Hem] emet jb]. 


5 
e. matr, [O) > та. m] | 





1 , 
tog maalt, 3) is а convex function of log r can be prcved in the following manner: 


1 k+l 
d(log то (и) _ Ser 2H | fte”) NOS t IE Кеке“) Рао dd 
dilog ry ay } 





Tak Hi 


1 kal | E 
E | f (re) Pde — нг а | и ІД нац 


ma) 7 





= [ 24 -#+1 | 


maar) 


15 


which increases with т for r2» ту, since ду!) 1s a convex function of nar). [Rabman 


(1956), pp. 194, 195]. 


2 \ 
Hence 215 тилт) > 0 for r> т, 
d(log ту? 


Lemma 2 If f(z), other than a polynomial, із regular for|z | < r and if 2 = ге", 


О<х< т апі Е > о 1, then 





5 T 

(1+) [1 (r, 2 - 

mai) = ps [a] 15 ph- 91 f ue pele es 
0 


Proof. By Cauchy's theorem 
3 1 T у AN 
Ifa) 15 S- l 1/0700) |40, 0<а<т, |zi-a. 
Also, by Holder's inc quality, [Titehmai: h, 1989, р. 882] 


2x 5 
| а) lx T !f(re ae) 
0 


тт, 


Ш 


= (2r) (т 2} 2л 


T | те“) "ад 


16 SATYA NARAIN SRIVASTAVA 


rè 
7 араа 


Now, атр те 2 +7 en) and intcgrating the above inequality with respect 


to 6 between the limits 0 and 27, we get 





. + ё 4% 
1 log x log us(z,f) log x 
Pe delere ne Y ные 
wef) = (Qn) i" mp log usta, р) g log = lcg pala f) . 
i 
But (zE) —0 as coo, since log us(z,f) is п steadily increasing convex 
log рна, ў) 


function of log 2. Hence for апу є > 0, we can find an д, = 2,(¢,f), su^h that for zm, 
(2 log x je 
lcg pa(z,f) 


Henze, pte, f) < мә canere] жа 





1 
(бај 


Integrating the above inequality after multiplying by z* with respect to x between 
the limits o and r, we get 


J entere © a zu (ее Лу pa(z + ze da 
0 


(ду)! zlogz 





3 г 
8 ЕР f jzk- 
* (3&)-1l jog т Д uz + 26)“ +, 





a 
5 t 
Hence, msx(r,f) - (146) {в p'r f) 


2 [i 
(zc + е) А 
(2т)5-1 log r | LEA pela + cela Ida, 


Proof of Theorem 1, We have from Lemma 1, 


T 

(rf) 15 

lo [4 Mi MUI Е. ] 
sup B ms, -s(r,f) 





e 
fü infuLÍLIILIU ul : 
és m log 7 А 
also from lemma 2. 
1 
log| 7| max rf 0) bl 

. sup їпь,к—ь{Т,]) 0 
lim inf my 
TRO log T 


and hence the result, 


ON THE MEANS OF AN ENTIRE FUNCTION AND ITS DERIVATIVES 17 


Corollary. If we take out the least value of 1/2 from the integral n (2,2) and 
proceed in the same manncr as in lemma 1 we obtain 


. AS 
Ee] G2) 


mair, f) > mun] Sr log T 


Аво if f(z, is an integral function of order о > 1, then 


mai lr f 0)) = d | f D (ze) [ct cde 
л E 


т теча 


< m ГГ: Met e= a два xd 


2z 


> А те | | (сей) "а а где, 8>1, 


0 0 


- nee 


f! a 





< O(r(et:-1), 





on using the result [Srivastav, 1957, p 3681, | t 


Hence, msx(r f") < ee (2.4) 


where A is independent of r. 
- From (2 3) and (2.4) follows 


ты} 


Toc ' ]ogr 





8. Theorem 2. For every entire function f(z), other than a polynomial, 


1 
[4 mlr f®) Y] 
Р patr, f) 
lim sup Шо Qe 5l 
TO logr | 


where, т tends to infinity through values outside an exceptional set of at most finite 


measure. 
The proof of this theorem is based on the following two lemmas: 


Lemma 1. For every entire function f(z), outside an é6xceptional set of at most 


finite measure [Valiron, 1928, p. 108], 
ут) d sc 
air) zm uf ( 2 i-e “ү, 
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where v(r) denotes the rank of the. maximum term in f(a), for | 2 | = rand 0 8 < 12. 


2x 
Proof. Wehave plr, f О) = af | | О (rest) 226 
7 


() 
„| | ЕС [ue не 
For ordinary valucs of т, we have [Valiron, 1928, p. 103] 


(1 (z) 


(=) 








= [1 + маву 75] 20) ‚ }Һ\<К, [а[= т. 
Hence, mlr fO) = af | + AG) (В) Ires у | ftre?) ра 
vr) Y 1 - Vor, 30 | 
= (20) team yf лег pao, 
0 


5 L 
is pi f90) = ptr ph P Ji- sq y 


Lemma 2. For0<8<1 and for every ров ов value of є, 


ur, f, << Aplr fret, 
where A is independent of r, 


Proof. We have 


pale f) = =f | f f (re^) fao ae 

















| f(re*^ Рад 





TERT 


< Hom f | (ева) Ва 
m 0 


= Aps(r,f)r(o*«- 0t 








' КӨ 

on using the result, | {O (те®) < O(riot«-1), 
f(re*?) 

Hence, plr, f) < Aplr, fret, 


Proof of Theorem 2. From lemma 1, PNE, 


log|r HE / B ] 
lim вир palt,f) slim sup log vtr) Ут) _ 
т+Ф logr т og г 
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Also from lemma 2, follows 


ьар wel) | 
пр J. 


hm sup 





ro logr uet 
(1) $ 
eps | 
Hence, lim sup Bal?,J) Si 
ro log r 


4. Let f(z) = 5 ала" be an integral function of order e and lower order A, M(r the 
тей 


maximum modulus of f(z) for |z| = r. Further it 1s known [Rahman (1956), p. 1921 
that if fiz) is regular in |e |< В, andit 2 = re”, Ox r«H,870, 


А (R3 — 7; (R^ – те, | f(Re®) do | 
Па) Р < i R?—2Rreos(6) - 9 4 r* —2Hrcos(0) - o 4- r*" (4.1) 


Frcm (1.1) and (4.1) it is evident that 


B ss MD = AT aq). 4.2) 


Theorem 3. If v(r) denotes the rank of the mazimum term, y(r), in the Taylor 


expansion of f(z) = Бач, | | = rendo <А Сесе then 


hm log man(7) Tis, i) (1- ~, 
To КӘР vis) kg t vir) lcg v = 


; log та ут) ( 1 1 ) 
1 © > 8 € , 
у г. v(r) log v(r) 





and if f(s) be of regular growth, then 


lim log nua) _ 
To v(r)logr 


lim leg marl) _ 
To v(r) log v(r) 


Proof. For functions of finite order wo know that [Srivastav (1956), p. 80-81] 





lim sup CEMO „ lim sup IEAM) <(!- х); (4.8) 
To (1) leg? rte v(r)log r e 
: log M(r) : log u(r) ( l1 1 ) 
] oN! = Jim oe ENE Le. 4.4 
ns e v(r) log v(r) To вор у(!) log у(т) А e den 
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Now, using the inequality (4,2) for | 2| = <r, k >0, we get 
pelz) = {M(zxi}}, 


or ds | mates d моја 
< жн {М (т)} pf ons 
_ 9 
= gei e 
л 9 
ie О, тт) s; МОЈУ. 


Since vir) is positive und unbounded function of г, therefore, taking limit: 
dividing by v(r)log r, we get 


lim up Eml) clm вир 2198 МИ) — ша sup Blog a(t) | 


To Ур v(r)log r To vw(r)logr тэ v(r)log r 
| li log ms(r) À 
Hence from (4.8), im вир SENT = 8(1– = ). 
. ‚ TO v(r) log T о 
è log me, x(T "Hoi 6 log M(r) c 6 log alr) 
Sgain; m sup v(r) log v L = а y(r) log vir) ея вор у(т) log у(т, ` 


TTherefcre, from (4.4), 


log m x (7) ( 1 1 ) 
lim Do ИКЕА — 2275 
е: вир vizi log vtr) ô z: 


B. Theorem 4. Let f(z) be an integral function, other than a polynomial, ој finite 
order e, n(r) denote the number of zeros of f(a) т [а | rand f(0) #0, Further, 





li 
O if ent UO = 
. ¿log та (т) gê 
him зар 98 aar) i 5 
then ипар = = eri (5.1) 
ЈЕ а а 
on | lim inf femal) ~, 8 5>0. (5.3) 


т+® rlogr 2' 
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Proof, (i) We have 


log p(z) = log 2- f” ijao, рав | 
0 


and since | f(ze) |? is а positive continuous function, we get from [Titehmarsh, p. 811] 


log us(z = log | Исе“) [40 
= ЈЕ tUn(üdtic-8log|f()| zr, (5.3) 
0 
on using Jenson's formula. 
Henoe, if lim int ?U = В, we have for апу s >> 0 and r2 ту, n(r) > (8—:)r", there- 


fore, log us/z) > 8(8— "1 +8 lag | f(0) |, for z > то. 
Since log ius(z)z*] is a positive function, 


У (из(т)а ах. > +f log (раја ја > d t- | r22) 


+8 log | #{0) |. (1-2 ) (са т - 1)-(1-2:)). 


Since log{ I z)orda} f log {ua(xja*{dx 
0 0 








therefore, log 2 Геге) > P шу + 8 log | f(0) ЈЕ = n) 


e(y+ 1) 
+8 (лов тез log ю)- (1-* 


А log тад (ер) > y +00) 


Proceeding to limits, we get 


linn ing log metr f) > 88 | 
To TP о(0 + 1) 





an 


ре; => 1, we have for any => 0 and т> м, n(r) > (L—s)rlog r, 
og 


(ii) Again, if lim 
ro 
and from (5.8), we get 
log m(x) > 20-9 f log tdt + 8 log Í f(0) { 
0 


= 8(1—s)(x log 2-2) +8 log | f(0) 
-or log fu(x)a*} > (1— e) (x log z—2) + 8 log | f(0) | - loga*, for трт. 
$— 2057 P—2 
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1 frog Íus(z)at]da > 1 | а {up(x)a*} dx 
0 n 


therefore 


т 


> га {ов r=" log „)-1(т-*=)- (т—'®)|+з1ов1/® | (1-5 ) 


+ в (log т a log r) – (1 = y 
T 


log { 1 f mòde ;> 1 f vs {us(z)a*}dar, 
0 


Since 
" T 
pu ; / 2 | 
therefore, . log | = ] пут)" ах} > &(1- ОБ log r— "е. log ra) — 3(; — 23 
+ò leg | то) |. (1 — =) (я p-o log x) -(1 -^) E 
T T T 
Henoe, log mar, f) 2350 — ут log r— gr} + O(1). 


Proceeding to limits, we get 


ro Tlogr 2 


f. 
6. Application. If A = 1 ара lim m : > 9, where n(r) denotes the ‘numrer 
Too 5 
< ma, (г, f?) «...... for all 





of zeros of f(a) m la| < т, then ms (r, f; < max(r, JP) <...... 


"> т, = T(N. 
From (5.2), we bave lim mf Я MATA ~, 5 
Tt logr 


THO 
Hence, from (2 3) with 8>1, - for rf, 
так T, f? F> log та, f) |> 1. 
msn, f) dr log r 


і.е. msan, f) < ть, llr, f?). 
Similar inequalities for higher derivatives give the result. 
S. K. Bose for suggesting the problem and his guidance in the preparation of: 


I am grateful to 


Dr. 
, this paper. 
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GEOMETRICAL SIGNIFICANCE OF THE 
ELECTRO-MAGNETIC ENERGY- 
MOMENTUM TENSOR 


By 
S. М. PaNpEY, Banaras, ОР, 
(Communicated by Prof Г.Т Narli.ar—Recewed—8September 9, 1960) 


1. Introduction. In general relativity gravitation manifests itself as a geometrical 
property of the four-dimensional space-time continuum. It is described by a Riemannian 
metric satisfying Einstein’s беја equations 


= Вау = R,,—3g,,E + Ад (1) 


which connect the distribution of matter and energy with the fundamental tensor g,, and 
its derivilives. The various forms of the energy-momentum tensor Ту, such as those 
for perfect fluid or electro-magnetic distributions impose certain restrictions on the 
possible geometry of space-time continuum, A geometrico algebruio- characterization of 
the: spuce-time continua of perfect fluids and the electromagnetic energy-momentum 
tensor has been studied in Феба] by Eisenhart (1924) and Lichnercwiez (1955) respec- 
tively. However, a global uspect of the space-time character:zed by these distributions 
is, sometimes, well-exhibited in terms of the Euclidean spaces of higher dimensions in 
which these can be immersed, If the minimum number of dimensions of a Euclidean 
space in which a given Riemannian space of n-dimensions con be immersed is n+p, the 
latter is.said to be of с] зва p. Karmarkar (1948) has studied some spherically symmetric 
metrics and by imposing class one restriction he has obtained Schwarzschild’s internal 
solution from, which the cosmological models of Einstein and de Sitter follow as particular 
cases. In a recent attempt towards а classwise study of Riemannian fourfolds, Singh 
and Рипдеу (1960) have been able to deduce a new form of non-static line.element for 
Lemuitre’s universe from class one considerations. 

While all these investigations show that a perfect fluid distribution of matter 
is compatible with class one metrics, 16 is found in the present investigation that the 
electromagnetic energy-momentum tensor ав a source term in Einstein's gravitational 
field equations is not compatible with Riemannian fourfolds of class one insofar as 
spherical symmetry is concerned. 

2. Electro-magnetic energy-momentum tensor and class one. A necessary and 
sufficient condition for a spherically symmetric metric - 


ds? = — e^dr! — (d  sin?6d93) +erdt?, А =, v = г), (2) 


to be of class one has been obtained by Karmarkar in the form 
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2 1 4 1 4 4 1 
B(1— 675 + Ват? (| — 873) (47. – Ti T) (137, тт.) =0 (3) 


where Ту“ are the components of the energy-momentum tensor defined by (1). We 
proceed to examine the compatibility of the electro-magnetic energy tensor (Eddington, 
1928). 

Ta = —F“F,, + 9 ЕВЕ а (4} 
. with (1), (2) and (8). 


1 2 2 3 
It 18 well-known that Ту, Ту, Ту, Ту, T, all vanish and 7, = T, in case of (2). When - 


these conditions are imposed on Т,” defined in (4) we get as а result the relations 


рар, + ЕР, = 0 (5) 

FUP, + ЕЕ, = 0 (6) 

FYR + Беба = 0 (7) 

Бар + ЕР, = 0 (8) 

FUP +Е Ева = (9) 

FUP +Р?*Е,, = FOP + PMPs. (10) 

Equations (5) — (10) imply that 

Р, =0, Е, = 0, Fy, =e PP, Fy, = 6 Еди (11) 


which signify that all the components of Fg beccme known if F,, and F,, are known. 
From (4) and (11) it сап be verified that 


T= Ty =D; (12) 
1 4 
Ti-T, = 0, (18) 
1 4 та 
Т.Т.-Т.Т, = 0. Е (14) 


When use is made of (12) — (14) in (3), the equation for A turns oul to be 
| ех = 1. ; (15) 
Substituting from (15) in the expressions for T," given by Tolman (1984а) and 

making use of (12) and (13) the following two equations for v are obtained : 

у’ [т+3^ = 0, (16) 
2A tv! фузјафујт = 0, (17) 
where вп overhead dash indicates a partial differentiation with respect to г. Consistency 
of (16) and (17) demands that A and v' should vanish together. Hence v is a function of 1 


alone and the line-element (2) reduces to а flat metric. Thus it follows that a spherically 
symmetric line-element of class one cannot describe an electro-magnetic situ. ation, 
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It may be pointed ont that some of the well-known solutions, obtained by the use 
of Maxwell's stress-energy tensor as source term in Einstein's gravitational field equations, 
are in agreement with this class property. For example one may mention here that 
Nordstróm's (1918) solution describing the gravitational field of a charged particle is a 
static spherically symmetric metric of class two.  Narlikar and Vaid:a (1947) have 
obtained a non-static spherically symmetric line-element describing the gravitational 
field due to a directed flow of electro-magnetic radiation.. The author has verified that 
the metric 80 obtained is of class two. 


It is important to note here that a disordered distribution of radiation which 
Tolman (1934b) defines as a special case of perfect fluid distribution characterised by the 
relation о = др is not incompatible with class one conditions. However, a confusion is 
created by the use of the terminology ‘a disordered distribution of electromagnetic 
energy’ insofar as the above established class property of electromagnetic distribution is 
concerned. That the case о = 8p may be described by a metric of class one does not 
contradict the above result, This is clear-from the fact that the stress-energy system in 
this case is not Maxwellian. 

While discussing the geometrical signifincance of the électromagnetic energy tensor 
one may also recall a result in unified field theory, In an investigation of the distant 
parallelism theory of Einstein, later modified by Levi-Civita, Tiwari (1951) has shown 
that the electromagnetic field totally disappears for a metric of the form 

48* = — Adz? — Вау? — Cd? + Ddt? 
where 4, B, C and D may be functions of the coordinates. 


The author is grateful to Prof. V. V. Narlikar for suggesting the above investigation. 
His thanks are also due to Dr, K. P. Singh for many useful discussions. 


Summary. It is shown that the electromagnetic energy-momentum tensor as a source 
term in Einstein's gravitational field equations is not compatible with spherically symme- 
trie line-elaments of class one. 
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ON SETS UNDER CERTAIN TRANSFORMATIONS—II 


Bx 
В. К. Lauri, Calcutta 


(Received —October 17, 1960) 


Kestelman (1947) has proved a number of interesting results in a bounded closed 
set of positive measure. He has proved among other results that if C is & bounded closed set 
of positive measure in E, (n dimensional Euclidean space) and p be any positive integer, 
then there exists a positive number 8 (depending on p and the measure of C) such that if 
Ay Aq -.. Ар are any p vecors in E, whose lengths |A.| satisfy the relations| А» | < ô, 
r — 1, 2, ... p; then it is always possible to find a closed set C, C C of positive measure 
having the property that if ¢ isin C,, then é+, is in C, r = 1, 2,...p. 


Now, if T,(r = 1, 2, ..., p) denote the transformation of translation 
' viz. а, = +a, t+=1,2,...n; а’ real, 


then we see from the result of Kestelman that he has, in fact, proved that if C is any bounded 
closed set of positive measure in Е, and р be any positive integer, then there exists a positive, 
number ó (depending on 2 and the measure of C) such that if 


= (af, ађ, ... а), r= 1,2, ..., p 
are vectors in E, satisfying |А | < 8, then the seb 
| СТ.СТ. С)... (T,0) = X 669] и 
is a closed set of positive uisi: E: 
(By T; О we mean the set of points эше are pointe of E, to which the points of C are trans- 


formed under T,). In a previous paper (Lahiri, 190598) we have proved a number of results 
similar to those of Kestelman by considering the transformation 7, in the following manner 


Т, 1. m = Ьа, i = 1, 2, ... п; br real 
where b, == 0 satisfy certain conditions. It has been shown. there that under some restric- 


tions on b, the set X as defined above forms a closed set of positive measure. 


1 
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In a later paper (Lahiri, 1959b) we enquire whether the set X as defined above forms 
a closed set of positive measure if T, (r = 1, 2, ... p) be any general linear transformation 
of the form’ 


nh 
24 =A аа: = 1,2,..m;af, real, 
where we assume that the determinant of the transformation, viz. 


ай; df 


аб 6 Oy 
is different from zero and aj, satisfy certain suitable conditions. 
Here also we have obtained some results parallel to those of Kestelman. 


In this paper we give a geometrical (functional) interpretation of the generalization 
of some of the results of Kestelman. Here the transformations are treated as functions 
and a neighbourhood is defined for each function making the aggregate of such functions into 
a topological space and the transformations (functions) are treated as points of this space. 


We proceed as follows : 


Let T' denote the linear nonsingular transformation, viz. 
, n . 
ax == арх Нан, += 1, 2,..п) ay real. 
=], = 


The transformation T' is completely determined by the matrix 


Qu СФ ... Qin Сал 
Og Су. Qan Сал 


Ani Ang >- Uns Ganga 


So, if we are supplied with the matrix of the transformation only, then we can fully under- 
stand the nature of the transformation and obviously different matrices give rise to 
different transformations. Now, if we write. 


Чп Фу ee Uu 01a 
f= Ga Са Oy Суза 


Чү Qng ... Qnn Ganyi 
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then by the function f we shall understand the correspondence between a set ОСЕ) and the 
set of points of E, which are the transformed points of C under the transformation T. 
Since Т is nonsingular f transforms E, onto E, in a bi-uniform and bi-continuous manner. 
We denote by f(C) the transform of CC E, by the function f (1.е. transformation T) 


Let now e be any positive number and 


Vi бя e. Vn Ving 
ЕЕ Tor Typ .. Ten Venta 

g == 
Tni Tng ee Tnn Unna 


where 2; be any real number such that 
0—6 < ty < а, tE, $—1,2,.,^; ј = 1, 2,... n+l. 


"Then we say that the set of functions g constitutes a neighbourhood of f. If the numbers 
a,j be such that a, < x, < фр +6, then we say that the set of functions g constitutes a 
right neighbourhood of f. Similar is the case for left neighbourhood. Having given є > 0, 
we say that the neighbourhood of f is determined by e. 


We present our main result of the note as follows : 


Let С be а bounded closed set of positive measure contamed in E, and let Г denote 
the function determined by tbe пх(п--1) matrix 


which leaves every point of E, invariant. 


Then corresponding to any positive integer p and positive number д (depending on 
p and the measure of С), we can determine a right neighbourhood of I such that if fy, fe 
...fp be апу p functions lying in this neighbourhood then that set 


C . f(O). (C) ... (0) 


forms a closed set of positive measure. 


In this connection, we recapitulate a well-known result (Caratheodory, 1948). If 
С denote any set in E, and T any linear transformation then 


m*(TC) = | D|m*(O) 


m.(TC) = |D|m«(O) 
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where |D] is the absolute value of the determinant of the transformation and m*(X), m.(X) 
denote the Lebesgue exterior and interior measure of X. So if C is Lebesgue measurable 
then ТО is also so and 


m(TO) = |D|m(0). 
We shall in the following understand by the term ‘measure of a set’ its Теренс measure. 


Theorem 1. Let C be a bounded closed set of positive measure contained in Ei, and 
р Бе any positive integer. Then it is possible to find a positive number М ата а number ё > 0 
depending on p and the measure of С such that tf I зз the function determined by the nx(n-+1) 
тайах 


and fı, fa, ... fy be any p functions belonging tn the right neighbourhood of I determined by the 
number 8/(M-+1)n, then the set of points belonging to 


€ .f«(0) - (б)... Јо) 
forms a closed set of posttive measure. 


Proof. Since C is bounded, there exists a sphere of radius, say M containing 0. 
Now, let U be an open set containing C and contained in the sphere of radius M satisfying 
the relation 
_2р—1 


m(U)— 0) < nO) gem a) 


Now, if 8 > 0 be the distance between C and U' (complement of U) then we may easily 
choose the set U in such & way that 


3p ` 2% | (2) 
Now, a point é of C corresponds uniquely to a point é; of f,(C) by the function fr 
(r = 1,2,... p). Since C is closed and f; is continuous, во f;(C) is closed. 
The function fr = 1, 2, ..., р) is determined by the matrix 
B bia uod. Maa 


1o owas җав 


бб Uy 6. D. Ua 
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where b}, has values such that f. lies in the right neighbourhood of J determined by the 
number 6/(M-+1)n. So, the determinant of the transformation which gives rise to the 
correspondence between C and f, (C) is 


T г 
ыы... M 


| s "T bry | 
Therefore the measure of f, (C) is given by 


m(f,(C)) =; | D, | m(C). 


Now, if X' is the corresponding point of a point X of С under ў, then it can be easily verified 
from the conditions on the elements бұ, characterizing the function fy that 


|X'—X| <8 (|x—y| denotes the distance between х and у) 


So, AY GU,rf—12/ p 


Again, it is easily obtamed that under the conditions on the elements bj, 


Р, > 1—7 bis b i, yes ba 


where the summation contains n!/2 terms. 


Again, it can be easily seen that 


а а а <(1+ ee) e 
ijr te It (M+1)n (M--1n (M+1)n 





ті 8 гы 8 
BO) Pap deg (1+ arcis! riis 
1—2^2.a1.6 


> 1— Ee from (2). 
3p 


So, m( fC) = Рт) > (= зр)" for r = 1,2, ... р. 
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Let Z denote the set f,(C).f,(C) ЈАКО). ... Л» (0), 0) = 


Then Z= Il 0) = U— X (U—f(0) 
fe тыў р 


mZ) > то) 2 (mU) т) > 40) [и (1-1) 0) 
| r-0 D 


2 


From (1) m(U)— —m(oy(1— ae z) 
m(C) 
So, | m(Z) > m(U)—(p+1) 21 > 0 


This proves: the theorem. 


Theorem. 2 Let C be a closed set of positive measure contained in an open sphere 
of radius, say М in Ел. Let {8} be any positive null sequence and f, be any function 
lying in the o neighbourhood of I (I as defined in Theorem 1) determined by the number 
8 (М- у. Then there exists a subsequence 8,3 of {8+} and a point X such that 


X ia inf, (О) f, (0) .f, (0)... 


This theorem can be easily proved. 
Theorem 3. Let C be a closed set of positive measure contained $n an open sphere 


of radius, say М in En. Then we can find a positive null sequence {ur} such that if (8,) be 
any null sequence satisfying 


0 < д, < min Џи, 2°41, nl] 


and f; 18 any function lying in the right neighbourhood of I (I as in Theorem 1) determined 
by the number 3,/(M-+1)n, then the set of points X belonging to 


C. fs, (О). fa, (0). fs, (0)... 


forms a closed set of positive measure. 
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Proof As shown in Kestelman’s paper we can obtain a sequence of open sets 
(Us) such that 


o 
U,20U,25U0U,25..; C= 0, 


fol 
and b (m(U,)—m(0)) < ©) 


If дт (20) is the distance between C and U; (complement of U r) then evidently {u} 
is a positive null sequence. We now choose the number 6, > 0 such that 


dts 1 
т = Dun |y, nA 


Then {ô} is also a positive null sequence. Let fs, by any function lying in the 
right neighbourhood of I determined by the number 6,/(M+1)n. Now, if X’ is the 
corresponding point ofa point X of C under ја, then it can be easily shown from the 
condition on the elements characterizing the function fs, that | ХХ "| < 6. 

So, №; (С) CC Ur for every ғ. 

Let 2 = fs, (0). fs. (C). fs, (0 
then if £ is in Z, £ is in fs, (C) for all n, i.e. dist (£, C) < 8, for all n. 


Since C is closed, £ is in C. 


Again, since eaoh fs, (C) is closed, so Z is closed. So, we are only to show that the 
measure of Z is positive. 


Now, Z= fs; (О). fe, (О). fs, (О)... 
= U,-3 (U,—fs, (0) 


T 2{ U,— 0,+ U,—fs,(0)} 
= C—X{U,—fs,(C)} 
+ mZ) > m(C)— E{m(0,)—m(fa,(C))} 
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As in Theorem 1, we see that 


тија, (C)) > 1—2". nl. 8) то) > m(0)— ZE 








m(Z) > m(C)—2{m(U,)—m(C)-+ О | 


mC) _ уо) —m(2)) > 0. 


This proves the theorem. 


Iam thankful to late Dr. H. M. Sengupta, for his kind help in the preparation of 
the paper. . 
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SOME ASPECTS OF SUPERPOSABILITY AND SELF-SUPER- 
POSABILITY OF FLUID MOTIONS IN GASDYNAMICS 


By 
J. N. Karur, Kanpur 


(Recewved—September 24, 1960) 


1. Introduction. To deal with the non-linear character of the hydrodynamical 
equations of motion, Ram Ballabh (1940a) introduced the concept of superposability and 
self-superposability in 1940. After that the subject was investigated by a number of authors 
Ram Ballabh [1-6] Bhatnagar (1957), Ergun (1954, 1944), Kapur (1959a, 1960a, 1960b), 
Ghildyal (1957, 1954), Prem Prakash (1954, 1955) and Strang (1942, 1948), all of whom con- 
fined themselves to the case of an mcompossible viscous Newtonian fluid. Gold and 
Krzywoblocki (1958a) in a very comprehensive study of the subject extended this concept 
to the cases of (i) incompressible non-Newtonian fluids (ii) compressible Newtonian fluids 
(iii) Compressible non-Newtonian fluids (iv) viscous conducting fluids under the action of 
magnetic fields. Independently of this Srivastava (1956a, 1956b) studied superposability 
in Non-Newtonian fluids and Kapur (1959b, 1960c), Bhatnagar (1958-59), Ram Moorty (1960), 
and Teeka Rao (1959, 1960) investigated superposability in magneto-hydrodynamics. Kapur 
(1959b, 1960c) also gave a number of applications of this concept to the case of viscous.con- 
ducting fluids. 


For steady irrotational motion of а non-viscous, non-heat-conducting gas, the 


equations of motion and continuity can be combined into the gasdynamics potential 
equation: 


(a?— p?) baa - (a? — pE) pyy - (a? ei) r 
— 29, фу фу —29:9.0:;—20,0.9 ,— 0, (1) 


where ¢ is the velocity potential and a is the velocity of sound given by 


a? = ak —}(y—1)(g2-+ 624-42), (2) 


where a, is the stagnation sound velocity and у is the ratio of specific heats. Gold and 
Krzywoblocki (19588, 1958b) referring to (1) state that “іп the given form, 16 is extremely 
difficult, if possible, to apply вирегровађшњу in the established manner". They therefore 
applied a limiting-process to solve the potential equation by reducing the solution of this non- 
linear equation to the solution of an infinite number of linear second order partial differential 
equations provided that the series of the sum of the solutions of these equations is convergent. 
The convergence is however not easy to disouss. 


2 
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In their discussion of the concept of superposability to viscous heat-conducting 
compressible fluid motions, Gold and Krzywoblocki (19582, 1958b) have been obliged to treat 
the variable Reynold number as a parameter with apparently no other justifiaction except 
that of mathematical tractability. This approach results in overspecifying the system and 
at the same time the system loses all physical meaning. 


There appears therefore the need for another approach to superposability in compres- 
sible fluids. The beginnings of such an approach have already been discussed by the author 
(Kapur 1960b) in another paper and this approach will be further discussed in the present 
paper where the results of the series of papers on steady rotational gas flows by Prim, Munk, 
Nemenyi Portisky and others (1945, 1948, 1944, 1958, 1958, 1944b, 1948a, 1948b 1948b, 1952) 
have been used. These results are stated to have been made possible by the reduction of the 
basic equations to the canonical form (1947) in the terms of the reduced velocity 
vector w = v/a’ where a’ ів the ultimate velocity along a stream-line and their approach 
was stated to be not practicable using the older formulation of the equations. It is possible 
however to obtain all these results and sometimes results of even greater generality with- 
out such a canonical reduction. Accordingly we have proved x aa from first 
principles those results which we require in our investigation. 


The main object of the present paper is to obtain three dimensional solutions from 
two-dimensional ones by using the principle of superposability. It is well-known that the 
number of exact analytical solutions for the basic equation of gasdynamics is very limited. 
Almost all of these—two dimensional vortex flow, radial flow, spiral flow, Prandtl-Meyer 
flow, flow past a wedge, Taylor-Maccoll flow past а cone and Ringlebs transonic 
flow—have been discussed in Howarth (1953), and Pai (1960). Again most of these are 
two-dimensional flows and we shall give a general method of constructing three-dimensional 
flows based on these two-dimensional flows. The number of exact solutions at our disposal 
will be considerably enlarged and these will be in terms of the velocity vector 
rather than in terms of the reduced velocity vector. 


x 2. The Basic Equations for Steady Flow. The basio equations for frictionless 
steady flow of a compressible fluid, in tensor notation, are: 





3 (ощ) = 0 Я [continuity] (2) 
pu See = SP [momentum] (3) 
Я =0 [energy] (4) 

р —'P(p)S(s) [state] (5) 


where u, are the velocity components and p, p, s denote pressure density and specific entropy 
respectively. Also-we are using the ‘product’ equation of state which includes the equation 
of state for a perfect gas as a particular case. 
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From equation (5) 


1 Op +. Pp) др S'(s) дв (6) 
p да . P(p 0% (з) Om. 





Multiplying both sides of (6) by w and substituting for 





др др дв 
Uy би. ' uar № За from (2), (3), (4) we get, 
Dus a ED) ОА 7 
да, P(p) pw да g 
or жы 8 (7) 
дњ Pr) | 
дж P(p) 
or p= TU) (8) 


where Y(p) and X are defined by (7) and (8). Substituting in (2) and making use of (3), we 


get дщ d 0X y ди 
+ (ү = —*7 <=. = 9 
kd ap | (ru o; 2. (9) 
This will be independent of the pressure if | 
P(p) 
Yi = 7 = а +b; 10) 
(P) = Ђ 0) ( 
where a and 6 are constants. 
Integrating (10), we get 
P(p) = А(ар--Б)"®, (11) 


where 4 is a constant. 


Thus we find that pressure can be eliminated from the basic equations only if the 
equation of state is of the form* 


P = A(ap+d)¥? S(s) (12) 


* For an alternative proof of this fact see Prim (1952), 
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and then the equation, satisfied by the velocity components is given by,. 


ди; : 
Mug = (13) 
ди; д Ox, 
(а—1) ww, de 1% да, дщ =0 
д 


This case, of course, includes the case when the gas is calorifically perfect, when 
the equation of state is given by | 





T 
р=Врув о (14) 


where B is a constant and P(p) has been taken as py 


3. The Definitions of Superposaibility and Self-superposability for 
compressible Fluids. For a frictionless flow which is not necessarily steady, the 
basic equations аге: 





др д EM r 
дщ ди; — др D 
08 08 i 
à 4 Qu 0 (4) 

p = P(p)8(9). (5) 


The six dependent variables w;, р, p, в characterise the flow. Two flows о Pp р, 
81; Ue, Ро, ра, ду, both of which satisfy (2’)-(5’) will be said to bo superposable or additive if 
1 is possible to find p, P, в such that и, и p, p, 8 is also a solution of (2')-(5'. A flow 
will be said to be self-superposable or self-additive if it is supersposable on itself. Thus if 
u, are the velocity components for a self-superposable flow, then both u, and 2w satisfy 
the equations (2’)—(5’), of course for different values of pressure, density and 
specific entropy. 


For gases for which the equation of state is represented by (12) i.e. for those for which 
pressure can be eliminated from the basic equations, equation (13) together with equations 
(2)-(5) shows that all steady motions for such gases are self-superposable. In fact it shows 
that if u; satisfy (13), then ku, would also satisfy it whatever positive or negative constant 
value k takes, thus showing the essentially linear character of the fluid velocity for 
steady flow of such gases and the non-linear character of the pressure, density and entropy 
fields may be responsible for the non-linearity of the basic system of equations. The principle 
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of self-superposability established here for steady flows of gases with equation of state (12) 
is also valid for the case of the more general equation of state (5). This we proceed to show 
in the next section where we shall also find the necessary adjustments in pressure, density 
and entropy fields. 


4. The substitution Principle for the Steady Motion of a Gas with 
Product Equation of State. If Kw, К'р, p,s' satisfy the basic equations (2)-(5), we 
get from (7), 





(15) 





We now choose К as such a function that 


a, 2% = y , grad К = 0, (16) 
да, 


where V is the velocity vector i.e. K is a function which is constant along the streamlines of 
the original flow. 


If both (7) and (15) are to hold 


j 1 
Ray (17) 
Substituting Ku,,1/K? p, p, s' in (2)-(5) and using (16), we find ^ 
0s’ | i 
А = 0 18 
5 (18) 
and 
; 1 
S(s') = те (Ss) (19) 


Thus we have proved that if w; p, р, s satisfy the basic equations (2)-(5) then Ku,, 
р, ШК? p, в' will also satisfy the basic equations where К ів any function which ts constant along 
the streamlines of the original (and therefore also of the new) flow and в’ is determined from (19) 
In the new flow field also the entropy is constant along the streamlines. In the new flow, the 
pressure field is the same as that of the original flow and the Mach number distribution 18 the 
game, for the local velocities of sound c, c’ and the Mach numbers M, .M' of the two flows are 
given by 








2 — др, _1 P(p), DN 1 P(p) _ Код 20 
| (5). р Po ^ ПОР) “° ae 





M Lv ultus | M'— xvi rui nul — М. (21) 
с с 
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Thus we have here an alternative proof of the substitution principle first established 
by Munk and Prim [29] for the steady flow of а non-viscous, thermally non-conducting gas 
with constant specific heats and later extended by Prim (1944) to the wider class of fluids 
having the product equation of state, viz; 

“If in any flow field, the densities are multiplied by m and the velocity by 1 where 
m is a positive function which is constant along each stream line, the new flow field will be a 
possible flow of the gas having the same pressure distribution and Mach number distribution 
88 the original flow". 

The equations (18) and (19) are important for they show that if the original flow ts 
isentropic, the new flow need nol be isentropic unless К is a constant and if the original flow із 
non-tsentropic tt will be possible to get one flow from tt which will be isentropic. The last state- 
ment is made possible, since s and K are both constants along the streamlines so that K can 
be chosen as a suitable function of s, and choosing К to be proportional to 4//S(s), we get 
S(s') = constant, во that the new flow in this case is isentropic. 


In particular by taking K to be a constant we get the result that all steady flows of 
а non-viscous, non-heat-conducting gas with a product equation of state are self-superposable, 
superposed flow has the same pressure field as the original flow and the density here is n times 
the density in the original flow where K із the ratio of the velocity in the superposed flow to that 
tn the original flow. Moreover the two flows are either both tsentropic or both non-tsentroptc. 


5. The Modified Substitution Principle for a Perfect Gas. For a gas 
for which P(p) = p!!^, (Т) gives 








д 
ищ да _ ap (22) 
С, E M 
да 


| If both w, p, p and Kw; Кур, Кур satisfy the basic equations and if К satisfies 
(16), then * 


Е, | OK ОК. 
г — a dnt. И — 
К? = = 0, 0, (28) 


да, 





80 that K, is constant. Thus for gases for which the equation of state is the product equation 
of state (5), but for which P(p) is some power of p, we have a modified substitution principle 
viz that tf for such a gas щ, p, p constitute a steady-state solution of the basic equations, then 
Ки, gi р, Кур also constitute a possible steady state solution where K, із a constant and К is 


a function which is constant along the streamlines of the two flows. The pressure field is now 
changed, though only in a constant ratio, but the Mach number distribution is still unchanged. 
The relation (19) is modified to 


x 


S(s') = = 8(8), (24) 
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so that the result that either both flows are isentropic or non-isentropic remains unchanged. 


By putting K, — 1 we derive the earlier substitution principle, as a particular case, 
though that principle was more general in another sense viz. that it held for a more general 
equation of state. The present principle, of course can be applied to a perfect gas for 

L 
which P(p) =p . 

The advantage of the present pinciple ts that it gives us а doubly infinite system of new 
solutions instead of the earlier principle which gave us only a single infinity of new solutions. 

Thus we can find new solutions in which the distribution of any one of the variables is 
unchanged or modified in а given ratio. In particular, we have that if wu, p, p, в constitute 


the steady state solution, then if K is а constant, and y is the ratio of specific heats, the 
following are also solutions, 


G) %, Kp, Kp, в'; 8(8') = EX" Sa). (25) 


In this case the velocity field is unchanged 
(ii) Ки, p, К?р, в"; 8(8") = K- Ss). (26) 


In the case the density field is unchanged. 


DE 1 ar al 

(iis) Kus pa Ps Bs 87; 88") = ру 80). (27) 
In this case the pressure field is unchanged. 

` é сй. ЈА 

(iv) Ku, К" p, K' ^ р, в. ` (28) 


In this case the entropy is unchanged. 


Alternatively we can apply the concept of self-superposability which we have so far 
applied only to the velocity field-to the other fields viz. the pressure, density and entropy 
fields also and say that for the steady frictionless flow of a perfect gas each of the fields-the 
pressure field, the density field and entropy field is self-superposable for if w, p, р, 8 is solu- 
tion, then the following are also solutions : 


1 


(v) Kp, Kp|m*, тица", S(si*) = Кі У [m*. S(s) (29) 
е EM РЕ 
(vi) Kp, K*mp, тщ, 8°; S(s*) = КЇ Ут Y S(s) (30) 


zd. gg E 2 m 
(vit) Кет "е ^o Bum Ye с p, mp (31) 
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Thus we may say that in the steady case, taken separately, each of the fields, pressure 
field, density field, entropy field is linear, but taken as a whole, the system із non-linear. 


Another class of gases for which the pressure field can be varied are those for which 





P(p) = е", (32) 
for in this саве (7) becomes 
а = ed (33) 
ш, 1 np pr 


| If both щ, р, p and Ки, Кур, Кур satisfy (2), (3) and (33), where К satisfies 
(16), we get 


1 p 
К? = KK’ Ку = constant (34) 


so that for such gases if u,, p, p constitute a steady solution, then 


Ku, K Кус" p, Кур (35) 


also satisfy the equations of steady flows where Кә is a constant and К is constant along the 
stream lines of the two flows. 


6. Motions Superposable on Steady Two-Dimentional Flows. For steady 
isentropic flow of a perfect gas, the equations for determining the velocity components are : 


—(y—1) Uy Uys] =( cc ) ЈЕ (86) 


where а comma followed by а suffix denotes differentiation with respect to the талаша 
corresponding to that suffix. 


ЛЕ [u (X1, 24), ито, 25), 0] satisfies these equations we get 








—(y —Y)[uy, tu, иу, а] = [wit кш Ri atta, VI, (37) 
21 28 
—(у—1)[ шуш, зима] = | Portas etsi tal от (88) 


If [0, 0, u(z,, 2, z4)] satisfies (36), then we get 


шз=0 (39) 
so that и; has to be a function of z,, 2, alone. 





SOME ASPEOTS OF FLUID MOTIONS IN GASDYNAMIOS 103 


Tf [0, 0, изба, 25)] is superposable on [и (21, £e), (жү, 25), 0] then. [u;(z,, 25), Ua(1, Ta) 
чз (21, x,)] must satisfy (30). This gives 





2 2 
э йн ee RM [ы Real абы ttg] . (40) 





2 2, 
—(у—1)[ ш us usi] = Б: каба aa ee жыйы) „ 6) 


—(у—1)[иу из, у иди 4] = 0 (42) 


If (42) is satisfied, then (40), (41) reduce to (87) and (38) . 


Г 


Thus the condition that the flow [0, 0, u4(z;, z,)] may be superposable on the two- 
dimensional flow [и (21, £a), us(a,, %а), 0] is 


WW + Mg, = 0 (43) 


1.0. му is constant along the stream lines of the two-dimensional flow so that wu, is 
a function of у, where у is the stream function of the two-dimensional flow. Thus if 


> > mid 
Ча = Urti + Ug ty | (44) 


gives a possible two-dimensional flow, then 


P. шу сту => 
Qs = Uy ty и + P(t, (45) 
is also a possible three dimensional flow of the same. gas. 


This, of course, includes the case when ¢(7r) is a constant, a case which Portisky 
(1945) obtained from Newtonian relativity considerations. 


On making use of (43) and the fact that u, is not a function of z,, we find from equa- 
tions (2) and (3) that pressure and density fields‘ are the same for both the flows (44) and 
(45). Also from (5), the new flow is isentropic. 


We could have also expeoted this result, on intuitive grounds from Newtonian rela- 
tivity considerations for-1f we bave a steady two dimensional flow and if the fluid moves in 
a direction perpendicular to the plane of motion in such a way that the velocity 1s the same 
for points on the same stream line, though it may vary from streamline to streamline, the 
flow pattern m any plane parallel to the original plane of motion will remain unchanged. 
This consideration also shows that the principle will also be applicable in the case of 
incompressible flows and for both incompressible and compresible flows it gives a means for 

| generating three-dimensional flows from two-dimensional flows. 


3 ^ 
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The vorticity components of the three-dimensional motion are 
Ма — Uggs 13 — Uprup Una 
and these will vanish if 
па = 0; 210; Uy — Upp = 0 (46) 


i.e. if из is constant and the original flow is irrotational. Thus in general, the new velocity 
field will be rotational even if the original flow 18 irrotational. 


For seeing whether the principle is applicable for non-isentropic flows, we obtain the 
conditions that both the flows I and II satisfy the basic equations (2)-(5) where those flows 
are characterised by : 


I: (21,2), wu 25), 0) (ть. та); рај, 25); 8(ш,х,). (47) 
п : az, Xs), (т, 23), из(®\, 24) Н p(x, Ta) ; р(®ү®у) ; a(x, у). (48) 
If flow П satisfies the basic equations, then 


Ри ива) Кир, Чар == 0 (49) 


РА 


Ру у у ита) = — ) рии) = — Ра) рии H Ugga) == 0 (50) 


Uy 811-4183 = 0 (51) 


If flow І is known to satisfy the basic equationa, then flow Н will also satisfy these 
equations if the additional condition 


UUs Ustiga = 0 (52) 


is satisfied. (52) is, of course, the same as (43). 


Thus we find that for both isentropic and non-tsentroptc two dimensional flows, a third 
velocity component wus(z,, %q)= ф(у) normal to the original plane of flow can be super. 
posed on the two dimensional flow field and the superposed flow will give a possible flow with 
the same pressure, density and entropy as in the original flow. 


It is easily seen that the above principle will also hold for diabatic flow with О, the rate 
of heat addition, as a function of жу, жү. 
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7. Construction of Three-Demensional Exact Analytical Solutions from 
Two Dimensional Solutions and Properties of these Solutions. We now combine 
the principles of superposition and substitution to construct three-diemensional solutions 
from two-dimensional solutions, Thus let 


> = Э 
ду = (Фу, 25) Нат, 23) % i^ (53) 


be a two-dimensional motion with stream function y with pressure p and density p, then by 
the substitution principle, a more general two-dimensional flow is 


> > > 
Ча = Ју ји «Думы , (54) 


and for this flow the pressure will be Kp, the density will be Kp/{f(y)}* and. the entropy 
&' will be determined from 


S(s') = SED 


FOT) S(s) | (55) 
where K is any constant. 

Now by using the superposition principle, the three dimensional flow, 
& 


> > > 
да = f(V)uyfs ЕД Ua tet oly) (58) 


will also be a possible flow with same pressure, density and entropy as for flow (54). 
Moreover from Croccos’ vortex theorem 
> > | 
qxCurlg = grad (+ 4) — T grad s 


or 


> ә 
Охо = grad h, —Т grad s, 


— 
where o is the vorticity vector, Љ ів the enthalpy, № is the stagnation enthalpy, and. T is 
the absolute temperature. 


For flow (53) 





> > 
à = Y.P | 14 2 
axo, = gad( 2 7 Tie) 28 grad s (58) 
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and for flow (56) 


> > , p 1 ; 
@ X иң = 269) $Y) god - (25 Pest nf grad y 
ОЛОР) grad (etar) (59) 
Р К” tg'(s) grad a— 2- S(s)f(J) grad ү]. 
Pg USC) rod e— уру SIX) grad V] 


We now wish to investigate the conditions under which the new flow will be (i) irrotational 
(ii) Beltrami (11) isoenergio (iv) isoentropic. р 


(а) The new flow will be irrotational if 


oe > > > 
©з = ourl qs = Ф (И), Ф (И) ба 


HY ша за) иза шаа) & = 0 (80) 


If the original flow is irrotational, the new flow will be irrotational only tf both f(y) and 
p(y) are constants. If the originl flow is rotational, then the new flow will be irrotational 
if P(r) is constant and f(y) is determined from 


РО) Тл ta — и +h) Гиз, 1— 3] 0 (61) 


For this to be possible it is necessary that 
0,10 — Pat] (gs 1 а) 


should be a function of у alone. Thus tf the original flowis rotational, itis not in general possible 
to find а flow from it, which will be trrotational. In one particular case, however, it will always 
be possible viz. in the case when the velocity magnitude із constant along the stream lines for in this 
саве density, pressure, temperature would also be constant along the streamlines and in the 
two-dimensional case the vorticity is also constant along streamlines. [Kapur A]. 


(b) The new flow will be tsentropic if (3) either the original flow із tsentropic and 
Јој) = constant or | 


(ш) +f ДИ) is chosen proportional to S(s) which itself is a function of y since s is 
constant along the stream lines and is thus a function of у. 


(c) The new flow will be tsoenergic if 
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either ($) the original flow $8 saoenergic and Фр) = constant or (53) Lf) Pho t+ ley)? = 
constant — A (say), 
where №, ta the original stagnation m Since A, is а function of y, it is possible 
to choose f(y) and ф(уг) to satisfy (62), but the constant stagnation enthalpy А has to be 
such that A—A, [f(V)]? is positive throughout the flow region. Thus given A, the function 
f(y) cannot be completely arbitrary or given f(y), A cannot be arbitrary, but A must be 
greater than the maximum value of # [У in the flow region. 


(d) The new flow will be both isentropic and isoenergic of [f(J^)]* is proportional to 
S(s) and (62) is also satisfied. Thus both Ду) and (y) will be determined. Thus, given 
a two-dimensional flow which may be neither isentropic nor tsoenergic $ ів possible from tt to 
generate а three-dimenstonal flow which will be both isentropic as well as ssoenergic. From 
Crocco’s vortex theorem it then follows that we can always generate a three-dimensional Beltrami 
flow from a given two-dimensional flow. | 

(e) Fora flow to be Beltrami, it is sufficient but not necessary that it should be both 
isentropic and isoenergic. Thus if the new flow is Beltrami, the necessary condition is 


grad HE (Y) hti] ТОЈ) grad s' = 0, | = 4 (63) 


where к 


[ЛР 


Since Т is, in general, not a function of у, it follows that in all cases where Т 18 not 
a function of у, the new flow will be Beltrami only tf it із both isentropic and isoenergic. If, 
however, the velocity magnitude is constant along streamlines tn the original flow $.e. tf in the ori- 
ginal flow the streamlines are concentric circles or parallel straight lines, [Kapur A] then 
T will also be constant along the streamlines and this will be a Jincan of y. In this case (03) 
gives the only restriction on f(y) and olp). 


(f) The density in the new flow is e/[f(yr)]* and f(y) can be chosen to make it constant 
provided p is constant along the stream lines of the original flow. Similarly the square of 
the velocity in the new flow is [f(V.)*g?--[G(V) * where g is the velocity in the original flow. 
This can also be made constant in-the new flow provided g is constant along the stream 
lines of the original flow. Similarly since the temperature in the new flow is ТД), it is 
possible to find a flow in which the temperature is throughout uniform provided in the original 
flow the temperature is constant along the stream lines. Thus we find that 1f in a two or three- 
dimensional flow the velocity magnitude 13 constant along streamlines, then it is possible to 
generate from them three dimensional flows for which either (+) the velocity magnitude is constant 
throughout the region of flow or (ti) the density ts constant throughout this region of flow or 
(itt) the temperature із constant throughout the region. The second case is of special interest 
as it provides a means of connecting compressible and incompressible flows. 

Thus we find that flow fields in which velocity magnitude is constant along stream- 
lines are of special significance. Prim [38] proved that for two-dimensional flows of this 
type, streamlines must be either concentric circles or straight lines and axially-symmetric 
flows of this type must be purely axial. Alternative proofs of these facts together with a 
discussion of the three-dimensional case have been given by Kapur [A]. 


Se!) = E S(s) (64) 
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8. Some Particular Examples. (a) Generalised vortex flow: 


For a two-dimensional vortex of straight K, 
z PANE ers a 
g = 0.4 5 ig 0.4. (65) 
Since the streamlines are circles, for the generalised flow 
> > > 
vu fe ene (66) 


The streamlines are the intersections of the surfaces, 
т = const, f(r)y-32—4ó(r)Ü = const. (67) 


i.e. they are right circular helices. 


If p, p, s denote the pressure, density, entropy of the original field, then for the new 
flow the pressure is proportional-to p, the density is proportional to p/[f(r)? and the entropy 
8'— 8, is proportional to —log[f(r)?. Since for the original flow the density vanishes inside 
a certain circle, for the new flow also it will be vanish inside a circular cylinder except in one 
case i.e. when f(r) is chosen as 4/p, so that the density is throughout uniform. The pressure 
will increase with r. 


(b) Generalised radial flow : 
For two-dimensional radial flow. 


= с Saat на 
= — | s 68 
CEN TUER tr+0. %-0 . tp | (68) 


where с is a measure of the strength of the source and p/p, is given by 





2 
Lx gi р 1 69 
n dng? as и [1—(o]p9) 3] 58 


The generalised flow is 





ro dieu s ae rd 
9 = АРХ 1-0 . 292-9 (n), . | (70) 
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We can choose f(r), (г) so that in the new field either ‘the velocity magnitude or 
the density or the temperature are constant throughout. 


(c) Generalised Prandil-Meyer Flow. 


Here for two-dimensional flow 





> > > > 
9 = Фаз Bin ЛО +,-|-ЛЧтаз соз AD 4,0. +, (71) 
Pol cos A0, D — сов? 1240, @ дов Аб, Eri cos? AG ^ (72) 
Ра Ps a, T, 
у =— — Ада f 008 28+110 |. (73) 
where 
—1 1 1 
Д? = am 3 = —— ,; 2 == — ^ у 
Ут, ведрина (74) 
The generalised flow is 
ә Р | > > | 
9 = Ў) Imas sin АӨ Е ДУРА qua 008 AD %--$(У/) is | 08) 


The stream lines are given by 


= —Qmaz Г 00828 A0 = const (76) 
2 = const [seh A0 ад | | (77) 


The properties of the above three flows and other three-dimensional flows that can be 
generated from known two-dimensional flows by using the methods of the present paper 
will be discussed separately ın a paper which will also deal with the axially-symmetric case. 


МАТНЕМАТТОВ DHPARTMENNT, 
INDIAN INSTITUTE OF THOHNOLOGY 
KANPUR 
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LAMINAR FLOW OF NON-NEWTONIAN FLUID IN 
CHANNELS WITH POROUS WALLS 


Ву 
S. К. Datta, Calcutta 
(Received—October 26, 1960) 


In this paper laminar flow of a certain type of non-Newtonian fluid between two 
porous plane parallel boundaries has been discussed assuming constant suction at both the 
walls. The type of fluid taken satisfies the stress-strain relation postulated by Rivlin (1948). 


Introduction. Laminar flow of a Newtonian viscous fluid near the entrance end of 
channels of plane parallel porous walls has been considered by Berman and Yuan for small 
values of the cross flow Reynold’s number R (Berman, 1963) and also for large values of R 
(Yuan, 1966). In both the cases a suitable form of the stream function has been chosen 
and then solution has been obtained as a series in В (if R is small) or as а series in R~ (if 
R is large). In our case we have taken the same form of the stream function and have shown 
that for such a choice the flow pattern does not depend on the cross-viscous coefficient, though 
the pressure drops in the axial direction as alsoin the transverse direction depend on the cross- 
viscous coefficient. This means that the normal stresses will also depend on the same, In 
& previous communication (Dutta, 1960) we studied the same problem for a fully developed 
laminar flow. ` 


Equations. The stress, rate of strain relation for a certain non-Newtonian fluid 
as postulated by Rivlin (1948) is 


ту ==— рду 206, + А itive ess 
ay == io а), 1) 
ду = 1,%=Ј;=0,% 2), 


where д is the viscous parameter and д, the cross-viscous parameter. This equation oan also 
be obtained from the Oldroyd's formulation (Oldroyd, 1958) for an elastico-viscous non- 
Newtonian fluid by the choice of the suitable values of the constants. 


We consider a channel of rectangular section, the distance between the porous planes 
being taken to be very small in comparison with the longitudinal dimensions of the walls, 
This condition enables us to treat the problem as a two-dimensional one. Both channel walls 
are taken to be of equal permeability, so that the velocity of the fluid leaving the walls may 
be taken as independent of the distance from the entrance end. We choose Cartesian co- 
ordinates with origin midway between the walls y-axis perp. to the walls and z-axis parallol 
to flow. Taking A = y/h, 2h being the distance between the walls. We have 


4 
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where д/р = у, Kelp == у, and the equation of continuity is 


But by virtue of our foregoing assumptions we may take 


до __ 
z7 


The boundary conditions are 


v(z, 0) = 0, 


v(z, +1) = v, = constant. 


(8) 


Since the flow is symmetrical about & plane midway between the porous walls, we consider 
the flow over half the channel, i.e. from the middle to one wall. We take the stream function 


as 


h(x, A) = Du (0) — ват ДА), 


(9) 
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where 10) is the average velocity (axial) at = = 0. So 
ча, A) = 10) tn] АЈ FA). | (10) 


А) = v, ДА). (11) 


Substituting (10) and (11) in in (2) and using (4) we obtain 


— ie --[wo- || па (а 591) rn] (12) 
120 um gr) ре [и Инда“ | (13) 

and the boundary conditions (5)—(8) give 
10) =f = Го) =0, f1-1 - (14) 


From (12) and (13) by a suitable reduction, we get 
Rf? —ff" I+ w? == C, ‚ (15) 


where R = v,h/v and C is a constant. This is the same equation as obtained by Berman and 
Yuan. 


Case I. Во ife, is small i.e. В is small, we assume following Berman 
fefeBÁÉGOEW Se, (16) 
O= ORO аб: 4 (17) 
with the appropriate boundary conditions 
= fill) =fx0) =0, n>0, 
A021; f0-29 21 


The first order perturbation solutions for f(A) and C are then (Berman, 1953) 


JOA) =F 2828) + gu (AS 2A—M), 


OW =—8+ = R. (18) 
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The expression for the axial pressure drop can be obtained from the equation (13) in. the form 





700, А) (=, А) | = [8-18 (2)] 
30v (0) AN, Му, \ № 
als. > 2/3 _ at 
x | 3+ ŠI во твЕ 18KRA 2 то) |, (19) 
where 
— Yet _ УВ 
hy BC 


Thus if v, is finite and => 0(у), we may neglect the product KR in (19) and write 





pO, А) —р(а, А) ____1 (2)[s- 16R 


io ONE x. (+) L-9* з5 8-189] (20) 


where ЈУ у, is the entrance Reynold's number = 4^u(0)/v. This expression contains К and 
so depends on the cross-viscous coefficient д,. Thus, although the flow pattern does not 
depend on y,, the pressure drop is a function of the same. We can also see that the stress 
components Tey, Tas, Туу at the walls are given by, 

a(0)? = (#)] [= Вад 5] 





_ Lu MBR[3 R | 
[т фр (О) у „+1 =— Pen = ye (p+ ggg ON TAB} af 


EENE и (0A? — 729-2) | + 


NR, 


и 


+8] 





у _ hes аа R m 
р 0) ал = уз (20 A8) as, (29 —1A1—2) 


+1 


+8K| =n 


п (ot 7as—2)} i 


а-я, (1) 


x[-9- 48А-195 | _, 


= EL 
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where 


Бери РА 2-00 oe; [IEF ) || -3+ 8] -y a-2y- 


where C is an additive constant. Here we have neglected powers of В higher than 1 and 

the products KR eto. 

This shows that 7,, is independent of К, and та and ty, depend on the same. 
Case II, If В ів large, then following Yuan we write the equation (16) as 


Л ГА 1 и’ __ 
ИЕ" = а, 


(21) 
and assume as a solution 
1 1 
= Ки Лт. , (22) 
1 1 
я = NUS род ра 04 +... (23) 


Substituting (22) and (23) in (21) we obtain as the first order perturbation solution of (21) 
as (Yuan, 1956) 





FMA) = sin 5 MI {1438+ л | ed \ TÀ 


sin TE 


: ла 
+! 0:002--7 log 





tan ТА | 2 sin 18 — А oos a} 1824 вила | 


а RU en 
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Thus the expression for the axial pressure drop is 


ЕЕ 


neglecting 1/R* and terms of the same order. 


This also depends on К i.e. оп v, We can obtain the expressions for the stress components 


as in the previous case and show that та; and ry, are dependent on К. 


І express my thanks to Prof. Dr. B. Sen, Р.Во., E.N.I. for his kind help in preparing 


this work. 
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Abstract, In the present paper, the steady, incompressible flow of conducting viscous fluid 
сд an annular channel with fixed boundaries, in the presence of radial magnetic field has been disoussed. 
Changes in the velocity profile are shown with the variations in the (1) strength of the magnetic Неја 
(proportiona! to В) (ii) distance between the two cylinders (proportional to 9). 


A pair of the points of inflexion in the flow region may exist and the condition for their 
existence has been obtained, velocity profiles are asymmetrical and 1 has been shown that the interaction 
between viscous and induction drags causes the shift of the point of maximum velocity towards the outer 
boundary. Induction drag is dommant over the viscous drag in the proximity of the inner cylinder, 
while the оазе is reversed near the outer cylinder. This gives rise to a sharp velocity drop near the 
outer boundary, 


Introduction. Recently Pai (1954) has discussed the flow of a conducting fluid 
in a tube in the presence of a radial magnetic field and pointed out the similarity in pressure 
distribution with the corresponding tubulent flow in hydrodynamics. Globe (1959) has 
disoussed similar problem in a cylindrical annulus. Magnetic field has been taken to be radial 
but is different in the form from Pai’s field. The Globe’s problem does not appear to have 
been discussed in full. He has obtained the solution of the system and the limiting form of 
the solution for parallel plates. The behaviour of the solution is quite interesting which 
has been discussed in the present paper. 


`. It has феей shown that for a given value оЁ 6, there éxist, in general, a critical value 
of f(—,) such that for all values of В > f, there occurs a pair of poirits of inflexion in the 
flow region. Similar result is true if 2 is fixed. The velocity profiles are asymmetric and 
the point of maximum velocity occurs nearer to the outer boundary. This fact is contrary 
to the results in the absence of magnetic field. A physical interpretation has been provided 
for it. There also occurs a boundary layer near outer cylinder for large values of f and 0. 


The solution of the problem with moving boundaries and their limiting cases have 
been considered in а, separate.paper. [Kapur and Jain (1982) Part II of the present paper] 
The behaviour of the solutions has also been discussed. 


Notation and the Solution. An infinitely long cylindrical channel of a and b as 
its inner and outer radii is considered. An external magnetic field (w/r, 0, 0) is present. 
Here о i8 a constant and characterises the magnitude of the magnetic field. Cylindrical 
coordinates (ғ, 9, 2) of a point have been used throughout the present paper. The common 
axis of the oylinders has been taken as z-axis. ` | , 
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On the assumptions that (i) fluid is viscous and electrically conducting (2) с is во 
large as to neglect displacement currents and convection currents (3) free charge density 
is zero (4) д (permeability) and с (electrical conductivity) are constant (5) Lorentz force is 
the only body force (8) the flow is steady and incompressible (7) axial symmetry is prevalent; 
(8) flow is parallel to the axis: 5 = (o, 0, 0,), and the applied magnetic field ој“ fixes the 
normal component of the magnetic field at r = а and at r = b for all values of z, and this 
is the only field impressed, it has been shown that (Globe, 1959, 1960, Kapur and Jain, 1960) 


(i) о, = 0,(r) 
Gi) E= |9,0, Нр) | 


on the assumption that H,=f(z)+¢(r) or Н, = /@)ф(). The solution of the problem 
is (Globe, 1959) 








Ра? _ 6? sinh (f In y)—sinh (£ In y/0) 
= кру |” sinh (Bn 8) | ~ 
H.=— Pæ о { yi — (0 д cosh (A In у/Ө)—@? cosh (В In y)— 1 4-6? cosh (В In 0) 2 
: (fi—4) А 2 f sinh (f In 6) } (2) 


where f? + 4, and 


Ра? gi 





Ly 2 = 
(£x ш y= a ш б— 2 + о (7) }. (4) 


The equations (1) and (2) can alternatively be written as 


~ 











Раа (s. (0—6-5)y8-- (08 —65)ју 
= Dy џи 08—0-8 j (6) 
Рай о [90 | yP(9-8—82)--y-F (085—602) —1 4 0409 40-8) 
H, их | 2 + | 09 


F(98—879) 


where —P is a constant pressure gradient in the radial direotion, 


у == "ја, 0-—bja, ү = p, Ё = Porh and А = l/(&muo) 
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Existence of the Points of Inflexion. The existence of the points of inflexion 
in the velocity profile in the flow region is generally connected with the problem of 
the stability of the flow. Therefore it is desirable to study their existence before discussing 
the stability problem. The necessary and sufficient condition for the existence of a 
point of inflexion at y = £ is 


Gas = 0 and ( ©) 52 0, where 1< £ < б. (7) 


The second derivative of the velocity is 








Po k BB —1)y83(8—6-8)-4- B(B-+1)(09 — 68) y 

dy? 4 | = 08—0-В | (8) 
where | Е = Pa?[n, ff? 55 4, во that 
4, _ _ ke @—l 1 

OE al In 0 at | (9) 

and 
dui _ k "lw 
(3i), xm [20 (1 я) (6: 108+2 у)]. (10) 





First, we shall determine the zeros of T As б> 1, ot from (9) is always 
negative and therefore velocity profile, in the absence of tho magnetio field, does not have 


any point of inflexion. 





(S), one has 


2AyB32— By28 C = 0 (11) 
where А = 08—8-8, В = f(f —1)(0?—0-8), and C = f(f--1)(08 —0*). 
If д> 2, then 

А> 0, В> 0, О;> 0. 


In addition, if р is a positive integer greater than 2, left hand side of the equation (11) will 
represent a polynomial z 


Ду) = —Вув -24yP13—0 
in descending powers of y. Clearly it has two changes of sign and therefore can have at the 


most two positive real zeros. Even if f = s , Where p and » are any positive integers 


and p > 2, Ку) can be made to represent a polynomial : 


f(x) == — Beta) 1 9 Aging - - he 2 ag 
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1 
by the transformation y" = к and the above result is still true. In fect f(y) corresponds 
to the flow which has to behave regularly . Therefore by the continuity of argument, f(y) 
can have at the most two real positive zeros for all real values of f > 2. 


If 2 = 1/n where n is any positive integer, one has 
A>0, Bg0,C<0, 


equality for B being true when n = 1. Using the same transformation, the polynomial in 
descending powers of x becomes 


—B43--2A4x192— C, 


which does not have any change of sign. Therefore for 2 < 1, there is no zero of f(y) for 
all values of 0. 


Similarly, if £ = 1+1/n, n being any positive integer other than 1, it can be shown that 
Ду) can be represented as a polynomial and shall have two changes of sign. Thus one sums 
up that for all values of f > 1, f(y) has atmost two positive real zeros for all values of 0 and 
for 8 x; 1, f(y) has no zero. Obviously the flow patterns do not have any point of inflextion 
for all Ө and f <1. 


Next question to be decided is about the position of the zeros. We consider 


/ Я E k 
(ара јуна 7 90-9) 





[B,(9-8 —08) + 28(—68)} el (12) 








d, k e T 
and ( са) y- =- ggg ну. Di OP) +2066 —6-9)] ~. (13) 
у 
where В, = f*—f—2 and [8 5 4. 
Also B,S0 when #22. 
div, ` div, | 

One observes that for all values of 0 and f? + 4, ( dy? ) y=1 and ( dyi bn are both nega- 
tive. Moreover when f* = 4 | 

Py k 

(St) ye =—g 0! in 0430-1) ... (14) 
dy k 
d Е ) M 22 —1)]. 

an ( d) узо =— с-ту 1 6190110) (5) 
Here again ( A ER and (=>) ЕГ are negative for all values of 0. Thus both the 


zeros of f(y) will be either within the annulus er outside it. 
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To decide finally that the two zeros of Eo are points of inflexion, one has to 











ду? 
consider the third order derivative of v, - 
dv, = Е #8+1(8+2) валу 2. 6- 
(E) эы = agre e (6° 0) = p15 1и—2)7 0 0-5] a9 
and 
до, _ __156%1п 0 
PNE ux a 
The zeros of a are given by = 
_ (#+1(8+2) 08—60? (18) 


(B—1)0—2) ®—0-8 


when 254 and 


ween f? = 4. One observes that right hand sides of an (18) [for all values of f except 
for В=1, for which there, obviously exist no zero of 5 and (19) are positive. Thus 


in the case of our interest [6 > 1], there exist at most one real positive zero of El 
7 





This zero of dv, cannot coincide with the zero of div, It is because du, has two 
Й Фуз ау? ау? 


zeros and since function is continuous there must occur an extremum of ae between 


its two zeros where necessarily it has to vanish. To conclude, we have proved that : 


(i) for P < 1, there is no point of infloxion, 


and (ii) for Ё > 1, either there is no point of inflexion or there are two pointa of 
5 inflexion within the flow region, . 
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FIG.A. rr SHOWS THAT FOR 6-2 
APPEAFRSONLY FOR Р Io 


‚ POINT OF INFLEXION 


Fig,2 IT SHOWS THO EXISTENCE or 
POINTS СЕ IMELEXION FOR OFF RRENT 
VALUES OF Р WHEN 9 =з 





Fig. 1. shows that for 0 = 2, the point of inflexion &ppears only for values of Ё 10. 
But as 0 increases the limit of f after which point of inflexion appears first (/,) goes on de- 
creasing. It becomes clear when figures 1 and 2 are compared. When 6 = 5, Ё, has been 
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shown (fig. 3) to be 4.16. For.all values of f > f, there will essentially be two points of in- 
flexion as the values of "s for y = 1 and у = 6 are always negative for all values of fand 

Y 3 . 


E vanishes as 2—0. The distance between the points of inflexion goes on increasing 
as increases. But it appears that the two points of inflexion will always occur on one side 
of the point of maximum velocity. 





Е 10, 25 IT GIVES Pc FOR $925 


Ё, for a given 0 or 0, for a given f can be obtained by the condition that the curve 
dy touches the y-axis. This will be given if the equations dy = 0 and dy =0_ 


have a common root. ‘The later of these two equations gives the equation (18). If that is 





also & root of the equation it = 0, the condition will be 


р ү 1 
_{ BUB+1) 08—0° 18-2  [(f--1f--2) 08—0* 28 — 
r =| 020 в) хя 2) тва) —° v 








which gives the critioal value of the other when either of the two 0 or f is specified. 


The importance of the existence of the points of inflexion in the velcoity profile 
lies in the study of the problem of stability of the flow. Since dv,/dy is non-zero except 
at the point of maximum velocity, therefore, there exist oylindrical vortex sheets within 
the flow, each element of which has the same vorticity, vortioity is decreasing from inner to 
the outer boundary. Appearance of the points of inflexion indicates the extremum of the 
vorticity. Though in the presence of dissipative terms, there is no theoretical result which 
shows the effect of a point of inflexion on the stability problem of the system. But there 
is an evidence of it (Berman (1958) in the experiments of the fluid flow through tubes with 
porouswalls. One more point is to be observed here. The induction current has only one 
component in Ó—direction and vanishes on both the boundaries. The point of maximum 
velocity and maximum induction current coincide. This point can also be said as a point 
of inflexion of the magnetic field vector, 
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- Asymmetry ‘of the Velocity Profiles. In the absetice of the шан field, the 
yelocity profiles are given by 





kp 681 | 
Ne | i In y— (у) 5 - (21) 
Since the relation 
vy, 6) == v,(14-0—y, 0) (22) 


is not independent of y, therefore the velocity profile given by (21) is asymmetrical. Not 
only this, but if f is also present, from the equations (1) and (3), one observes that the relation 


Ру, 9, В) = v, (1--9—у, 9, 8) (23) 


eorresponding to the relation (22) is not independent of ‘у. "Therefore again the velocity 
profiles are asymmetrical. We shall first study the asymmetry of the solution (21) and then 
see how this asymmetry differs from the asymmetry present in the solutions in the presence 
of magnetic field, 


The point of maximum velocity y, of the solution (21) is given by 


А (24) 


Table 1 gives the values of ул for different values of 0. It is clear from there that the point 
of maximum velocity always occurs nearer the inner boundary. Moreover as 9 increases 
there is an actual shift of the point towards the boundary. This asymmetry is due to the 
geometrical asymmetry of the problem. 


The asymmetry in the velocity profiles may be intuitively perceived by the following 
consideration. The radius of the inner cylinder being smaller than that of the outer one, 
the ratio of the viscous to the inertial forces is expected to be larger near the inner cylinder 
„во that the velocity gradient will also be larger near the inner cylinder with the result that the 
maximum velocity will be reached in a shorter distance from the inner cylinder as compared 
-with the distance of this point from the outer cylinder. By the same argument the larger 
the value of, we expect relatively closer to the inner cylinder will be the point of maximum 
velocity. If the geometrical asymmetry of the problem is removed that is if one considers 
limiting form of the solution (21) [Flow between two parallel plates], the velocity profile 
becomes symmetrical. 
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CABLE 1. It gives ym the pomt of maximum velooity for diffrent values of 98 ==0) shift d - T1 —ym 
from the middle of the annulus [It 18 positivo if the shift 18 towards inner cylinder and vega- 


tive, 1f the shift is towards outer cylinder.] and the percentage shift у = ns X100. Every qu. 


я 2 
antity has been evaluated upto three places of decimals. 

















5.8 


One can easily show that the maximum velocity can never occur at the middle 
of the annulus. If so, Ó is to be given by 
8-1 
«(8— 
‚ ЛӨ) =e бы) 
Consideration of the sign of the first and second derivative of f(0) gives it to be continuously 


decreasing function as 0 varies from one to infinity. Thus @ = 1 is the only possible zero 
of КО) and so the point of maximum velocity сап never be at the middle point of the annulus. 


—@ = 0. (25) 


Similarly the maximum velocity occurs at the inner oylinder if and only if 
Е(0) = 5 (9—1)—n 8-0 (26) 


whose first derivative F'(0) == (62—1)/6 is positive for all values of 0. Therefore F is a conti- 
nuously increasing function and 0 = 1 is the only possible zero. One concludes that in the 
absence of the magnetio field, the point of maximum velocity always lies between the inner 
cylinder and the middle point of the annulus. "Rc 


This asymmetry in the velocity profile continues to be present even in the presence 
of magnetic field, but surprisingly the point of maximum velocity is shifted towards outer 
boundary. This point is given by 


2y (0? —8-5)— B((8 — 0-9) у8—(08—6%уу-8] = 0 (ет) 
when 2824 "and 


Jing) m6 
Ту 6 —1 


when f? = 4. Had there been no #, for a given value of 0, the point of maximuni velocity 
would have occurred, near the inner boundary. The magnetic feld shifts the point towards 
outer boundary. The shift increases with the increase in 2. Fig. 4 shaws how the point 
of maximum velocity shifts towards outer cylinder with the increase in f. For 0 = b, В = 10, 
the ratio of the distances to the point ул, from outer-aud inner boundaries is approximately 
3/13 while for 0 = 5, В = 4 this ratio is 5/11. We expect, naturally, that for given 0 there 


(28) 
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will be some value of Д, for which the point of maximum velocity lies at the centre of the an- 
nulus. The values of such f’s are given by putting у = (14-0)/2 in the relation (27). We 
Bhall obtain 


28-1925 —1)1-4-0)92— удана — 1)(1--0)95---276 (998 — 0842) = 0. (29) 


Taking the analogy for f = 0, we except that for a specified f 0, the point of maximum 
velocity would move towards inner cylinder with the increase in 0. But, in fact, the result 
is just opposite; it still shifts towards outer cylinder. lt becomes clear when figures 4 and 5 
are compared. For 0 — 4, f — 10 the ratio of the distances to the point y,, from outer 
and inner cylinders is approximately 1/4. Since this ratio for 0 = б, В = 10 is 3/13, there- 
fore the point of maximum velocity for Ө = 5 is nearer to the outer boundary than that for 
0 = 4. The reason for it lies in the interaction of magnetic and viscous forces. 





{0 15 20 25 зо 75 40 45 5 
v ——Ó—— 
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FIG.8 VELOCITY PROFILES 
FOR Ф = 


The magnetic force is proportional to g?r АУ where JA is the characteristic magnetic 
field (w/r) perpendicular to the velocity field and the viscous force is proportional to ру +12. 
The ratio of the magnetio force to the viscous force is given by Hartmann’s number If 

dre а ОЛО ИА 
у? py rn? 





(30) 
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where 0 has been taken as the characteristic length of the flow. Magnetic forces and viscous 
forces, both are dissipative forces; velocity falls in their presence. But while magnetic 
force reduces the velocities considerably throughout the flow region, viscous force is more 
effective near the boundaries. The comparison of the values of M from (30) clearly indicates 
that it is greater near inner boundaries than outer one. Consequently, while the reduction 
of velocity near inner boundary is more due to induction drag, it is due to viscous drag near 
outer boundary. And we e&pect that the point of maximum velocity to occur near the outer 
cylinder. Expression (30) shows that for fixed w/v, the shift is again towards outer cylinder. 


Effect on the Thickness of Boundary Layer. As has been seen already, the 
interaction between viscous drag and induction drag shifts the point of maximum velocity near 
the outer boundary. Moreover, since shift increases with £ and 6, therefore for large values 
of р and 0, while іп most of the region of the annulus velocity increases slowly, it falls down 
abruptly from maximum to zero in a narrow region near the outer boundary. From Fig. 
4 one observes that for # = 10, 0 = 5, in less than one-fourth of the annulus the velocity 
falls from maximum to zero, while in more than three-fourths of the annulus it increases 
from zero to maximum. For f = 4, 0 = 5, the sha?p decrease is in about 3/8 of the total 
region. Thus, though the volocity as a whole decreases markedly for large values of A, 
but simultaneously it provides a narrow region near outer boundary where the velocity gradient 
is sufficiently large. It can be argued that the thickness of the boundary layer near outer 


oylinder is increased markedly when large magnetic fields are present while the inner boundary 
layer is almost unaffected. 


Since, if 2 is kept constant and @ allowed to vary, the point of asymmetry again 
shifts towards the outer cylinder. So we expect some further increase in the thickness of 
the boundary layer. In conclusion, when 6 and £ both are large, this region is more significant 
as could easily be seen by comparing the figures 4 and 5. 
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ON MUKHOTI'S THEOREMS 


By 
D. B. Goopnur, Tallahassee, U.S.A. 


(Communicated by the Seoratary Received—October 25, 1960) 


Generalizations of Rolle’s theorem, the law of the mean, Cauchy’s generalized law 
of the mean, and Darboux’s mean value theorem for ‘derivatives have been given by 8. A. 
Mukhoti (1958). These generalizations are of particular interest since they do not require 
the existence of the derivative and, hence, may serve us when there is no derivative. Тһе 
purpose of this paper is to strengthen and to extend Mukhoti's results. which are contained 
in the first parts of our first, second, fourth, and fifth theorems. Theorems of Karamata 
(1951) and Vuékovic’ (1952) as corrected by Utz (1955) are special cases of our results, 


Let f be a function defined on a finite closed intervala < x < b. Throughout this paper 
we will use D,f, Dif, D,f, and D,f to denote the upper right, upper left, lowe. right, and lower 


left Dini derivates [Karamata(1951), p. 188; McShane (1944), p. 87], respectively, of f. We 
use the name derivate to avoid confusion with one-sided derivatives [McShane (1944), p. 189; 
Vuékovic’ (1952), p. 159]. Та addition, if p and g are real numbers, we will use 


pD,f+qDif 


to denote any one of the expressions 


pau, тра рај ўар 


Also; we will use pD,f+-qDif to denote all four of the above expressions. Finally, we will 
use braces to indicate alternate statements. | 
Our first theorem is a generalization of Rolle's theorem. Since we require only that 


Ка) > Ка) Сао) < Га}, the first part of the theorem is slightly stronger than Mukhoti's 
Theorem 1 [Mukhoti (1958), p. 87]. This enables us to strengthen later results. 


Theorem 1. (i) If f is upper {lower} semi-continuous at each point of the finite 


closed interval a < x < b, if f(a) = f(b), and if there is an zy, а <a, < b, such that f(a) > f(a) 
(f(x) < fla}, then there exists а ё, а < ё < b such that 


DPS) <O < р (f) «9« DAE. 
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(и) If, in addition, f has а right derivate D,f and a left derivate Гу} each of which їз 
finite at Ё, then there are numbers p and q with p > 0, а > 0, 2+д = 1 for which 


Рр јео +ар је) = 9. 
Proof. ($) A proof of this part may be obtained by making minor modifications in 
Mukhoti's proof [McShane (1944), p. 88]. The details will be left to the reader. 


(it) If D,f(é) and D;f(£) are both finite, we see that D,f(é).  D;f(f) <0. If, in 
addition, D,f(£) = 0, we may choose р = 1 апа 4 = 0. If D,f(£) 5 0, we may choose 


— DO md g= DID 
PU DAO-D4O ~ * 7 DHO- DIG 
In either case gD,f(£)--qD,f(£) = 0 which was to be shown. 


Theorem 2 is a generalization of the law of the mean. Its proof, as well as the proofs 
of our remaining theorems, is based on Theorem 1. 


Theorem 2. (i) If f is upper {lower} semi-continuous at each point of the finite 
closed interval a <a < b and if there із an zs, а < zy < b, for which 





then there exists а ё, а < Ё < 0, such that 


pipe) 2:109 > раја (Бука < ФЛ < p gy). 


b—a _ 





(ii) If, in addition, f has a right derivate D,f and a left derivate Dif each of which із 
finite at ё, then there exist numbers p and q with р > 0,9 > 0 p+q = 1 for which 


SOAD орка) © ШЕ 


Proof. (i) A proof of this part may be obtained by making miner modifications in 
Mukhoti’s proof [Mukhoti (1958), p. 89-90]. The details will-be left to the reader. 


(ii) If D,f(f) and D,f(é) are both finits, the auxilliary function 


Fie) = f) O10 „ 
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satisfies the conditions of Theorem 1, part (ii. Hence, there are numbers p and q with 
20,9420, ptq = 1 such that pD,F(£)--qD,F(£) = 0. 
If follows that [McShane (1944), p. 191] 


1010) — pp,fe)+ Dill) 


a 
which was to be shown. 
Theorem 3 is an extension of Taylor’s formula with remainder. 
Theorem 3. ($) Let the function f, defined on the finite closed interval a < = < b, 


be such that tts first n—1 derivatives exist and are upper {lower} semi-continuous on a < x 
< b (only right and left derivatives are required at a and b, respectively) and let 


n-1 


4=0 





If v is a number such that п— > 0 and tf there із an x, а < ту < b, for which 


nol 
(b— zo) f (Xp) ДЕ о К > f(b) 


i! (0а)? 











(90 


n-1 


ГУ Ва ea) | 0—29*" к c np, 


ч иди 
2 il (b—a) 





then there exists а ё, а < ё < b, such that 


(n —1)l(n—v) 


АША K > D fU 
Eras 5? PU 


ру") > 


[ D, < ead К< реј). 


(8) If, in addition, f*- has a right derivate D, f ^? and a left derivate Df 7 each of 
which is finite at ё, then there exist numbers p and q with p > 0, q > 0, рід = 1 for which 


f(b) = Ја) 4- (b —a) (а)... —— ш 


b— n—t b— © "и = 
ар Гер, fo) + ai f" "Q1 
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‚ Proof. А proof of this theorem can be obtained by noting that the auxilliary function 
1 Toe s HI voes *. 5 ^ o us 
F(x) = fb) +f) + (0 —2)/ (2) +... 
(b—2)"3 pinay (b—zy* 
ш Л" ај» 


satisfies the conditions of Theorem 1. The details of the proof are left to the reader. - 


+ 








Our next theorem is an extension of Darboux’s intermediate value theorem for deri- 
vatives. 


| Theorem 4. (i) Jf f is upper {lower} semi-continuous at each point ој. the m 
closéd interval а < x <b and if k {з а number such that 


Dj) < k< Ва {DO > k > Ра), 


then there exists af, а ТЕ < b, for tohich 
DAE) < k Бле) Xf) < k < РАЛЕ). 


(ii) If, in addition, f has a right derivate D, f and a left derivate D,f each of which ia 
finite at ё, then there are numbers p and q with р > 0, q > 0, p+q = 1 for which 


»D,f(E)--qDif() =k (D. f(E)--aDv ЈЕ) = k). 


Proof. ($) Let f be upper semi-continuous on a < x < b and define the function 
_ F by F(z) = (а) — Ке. Then РР) < 0 and D,F(a) <0. Hence there is a point ё, a < ё 
<b, such that F has a maximum at ё. Then D,F(t) <0 < БР). It follows that 
D, f(E) < k < D, f(é). A similar proof may be given when f is lower semi-continuous. 


(и) НО, КЕ) and Di f(é £) are finite, then D,F(£) and D,F(é) are finite and by Theorem 
1, part (i$), there exist numbers, p > 0, q > 0, р + 9 = 1 such that TEST с = 0. 
Hence 


ИЛЛЕ = 


which was to-be shown: A similar proof may be given for f lower semi-continuous, 


"m Let f and g'be functions which are defined at each point of the closed interval 
а с ж < 5. We shall say that the derivate Df is partially bounded above a the DE 
Dg provided there exists a number £, a <§< b, such that [g(b)—g(a)]Df(£) 

[/(b) —f(a3)]1Dg(£). In addition, we shall say that D,f+D,f is partially а above = 
D,9-2- Dig mir there exist numbers р, q, and £ with p > 0, q > 0, ptq=l,anda<é 
<b such that 


[g(b) — XD, f дару < 1700) —/(9)12Р()+4ри()1. 
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Moreover, we shall say that the four expressions D, {+D, f are partially bounded above by 


Dg+Dyg provided each of the expressions is partially bounded above by Dg+Dy. The 
notion of partially bounded below is defined by replacing "above" by “below” and “<” 
by “>” in the preceding definitions. 


Our final results, which are extensions of Cauchy's generalized law of the mean. 
wil be stated as three related theorems. We shall prove only one relation from the first 
theorem. The other relations can besimilarly proved. The reader will observe that additional 
relations may be obtained by interchanging upper and lower and reversing the inequalities 
in the theorems. 


Theorem 5a. (+) If the function f зз upper semi-continuous and the function 
g tslower semi-continuous on the closed interval a < x < b, of f(b) —f(a) > 0, if g(b)—g(a) > 0, 
and sf there is an жу, а < хо < b, for which 


[g(5) —g(a)]Lf (xo) —/(«)] 22 1700) — Да) 1900) —9(9)) 
then D, f and D, f are partially bounded above by D,g und Dg, respectively, and D, f and Dif 
are partially bounded below by Dig and Dig, respectively 


(й) 1 f, in addition, the one sided wpper and lower derivates of f and g are finite «t Ё, 
then D,f+Dif is рата у bounded above by D.g+Dy and зз partially bounded below by Dig 
+Dy. Moreover, Df+Dif and Dif +Dif are partially bounded above and below by Dg4- 
Dig and Dj 4-Dig, respectively. 

Proof. (3) Let 


P(x) = [9(6) —9(«)]Lf C9) —Ја)] — L0) — (9) (2) gla). 
Since F satisfies the conditions of Theorem 1, part (2), there exists a ё, а < Е < b, such that 
D,F(é) «0x БЕ). It follows that the relations of part (+) hold. 


(и) Тан derivatives of f and у are finite at ё, then the function F satisfies the condi- 
tions of Theorem 1, part(i:). Hence there exist numbers p and q with p > 0,g > 0, 2+9 = 1 


such that pD,F(£)--gD,F(£) = 0. But [McShane (1944), p. 190] 


D,F(«) > 1906) —4(9)1D,f (x) ГЛ) Ле) 
and 


D,F(z) > [g(0) —9(2)1 Df (x) —L/(0) f(a) Dyle). 


Tt follows that 
[9(6)—9(a) еб.) DFE] < L/06)—f())L»Dig (6) -aDig(£)] 


which was to be shown. 
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Theorem 5b. (i) If the function f is lower semi-continuous and the de g i8 
upper semi-continuous on the closed interval a « x < b if f(b)—f(a) < 0, У g(b) —9(a) < 0, 
and if there їз an zy, а < zy < b, for which 


[g(b) —g(a) о) Ла) > 1700) — Ла) Јо) —9(9)1, 


then Df and D,f are partially bounded. above by Dg and Пу, respectively, and Dif and Dif 
are partially bounded below by Dig and Dig respectively. 

(8) If, in addstion, the one-sided upper and lower derivates of f and g are finite at rE, 
then D,f+Dif ts partially bounded above by Dg+Dig and is partially bounded below by Dg 
Dg. Moreover, Бури and D,f4- Dif are partially bounded above and below by Dig+Dy 
and Dy4-Dg, respectively. | 

Theorem 5c. x If the b adu) f and A are upper semi-continuous on the finite 


closed interval a < x < b, if f(b) fla) < 0, if g(b) —g(n) > 0, and if there ts an ара <a, < b, 
Jor which 


[9(b) —9(@) Lf (%o) —f(a)] 22 [Л6)—ЛДе 9%) ga), 


then D,f and D,f are partially bounded above by Dry and D,g, respectively, and Dif and ру ате 
partially bounded below by Dig and Dy, respectsvely. 

(8) If, in addition, the one-sided upper and lower derivates of f and g are finite at Ё, 
then D,f+Dif and D,g-- Dig are partially bounded above by Dyg+Dyg and D,f-- Dif, respec- 
pectively, and are partially bounded below by Dg+-Dyg and D,f-- Dif, respectively. Moreover 
D,f--Dif and D,f4-D,f are partially bounded above and below by Dg+Dy and DgiDg, 
respectively. | | 
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